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SECTION 1 


INTRODUCTION 


The primary objective of the Langley Atmospheric Information 
Retrieval Ssytem (LAIRS) is to make it possible to obtain 
accurate estimates of atmospheric pressure, density, tem- 
perature, and winds along the Shuttle reentry trajectory 
for use in postflight data reduction. One possible appli- 
cation of this information will be in the Aerodynamic Co- 
efficient Measurement Experiment (ACME) , which is part of 
the Shuttle Flight Test Program. The objective of this 
experiment is to improve estimates of the aerodynamic 
coefficients characterizing the Shuttle. 

The LAIRS Program operates in one of two basic modes , de- 
fault or adjusted. In the default mode, pressure, density, 
temperature, and wind values are interpolated for specified 
points along a Shuttle reentry trajectory using standard 
(default) atmospheric models stored in the computer. In 
the adjusted mode, atmospheric modeling parameters are 
adjusted to fit in situ measurements of atmospheric quan- 
tities from various sources, such as sounding rockets, 
balloons, or ground instrumentation. Several modeling 
approaches are available when LAIRS is operating in the 
adjusted mode, including the Jacchia-Roberts model for the 
upper atmosphere and both linear and nonlinear least squares 
polynomial fitting models for the lower atmosphere. These 
models are discussed in detail in this document. 

The purpose of this document is to describe the LAIRS Pro- 
gram and its operation. Section 2 presents the User's 
Guide, including descriptions of all input and output. 
Section 3 describes the overall system and includes subrou- 
tine descriptions for all of the LAIRS subroutines. Section 
4 contains descriptions of the computational models used 
to determine the required atmospheric parameters , and 



Section 5 describes the development of an optimal procedure 
for using the LAIRS Program with data from the first Shuttle 
flight (STS-1) . Appendix A contains descriptions of the 
COMMON blocks , Appendix B specifies the formats of the 
files that are created and used by the LAIRS Program, and 
Appendix C contains benchmark decks used for testing the 
LAIRS Program. File maintenance routines are described in 
Appendix D, and the results of noise analysis experiments 
with the LAIRS Program are discussed in Appendix E. 

The authors wish to acknowledge the assistance and guidance 
of Joseph M. Price, Robert C. Blanchard, and Harold R. 
Compton of the NASA Langley Research Center, who contributed 
much to the successful development and implementation of 
the LAIRS Program, They also wish to acknowledge the exten- 
sive contributions to this document by Ajay K. Kapoor, 
Richard A. Kuseski, and James 0. Cappellari, Jr., of CSC. 
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SECTION 2 


USER'S GUIDE 


The Langley Atmospheric Information Retrieval System (LAIRS) 
Program has been implemented on the Langley Research Center 
(LaRC) CDC Cyber 170 computer with the NOS 1.4 operating 
system. This section contains a description of the input 
required (both user input and data ) , the output options and 
reports, and the output files for the LAIRS Program. In 
thelse descriptions, it is assumed that the reader has a 
working knowledge of the computer and the operating system. 


2.1 USER INPUT 


The user input for the LAIRS Program is composed of the 
following three types of input: 


• ENTREE Trajectory File - This file contains Shuttle 
positions as a function of time along with the 
uncertainties associated with those positions. 

• Meteorological Data Files - These files contain 
observations of temperature, pressure, density, 
and winds made at up to three stations along the 
Shuttle reentry path. 

• Formatted Keyword Cards - The information on these 
cards will include the following: 

- Atmospheric model to be used 

- Type of profile to be calculated (ENTREE 
Program trajectory. Local Atmospheric Profile, 
or both) 


In addition, a previously created permanent file (KP) 
containing the 3-hour geomagnetic indices (K ) and the 

It 

daily Ottawa 10.7-centimeter solar flux values is needed 
for use with the Jacchia-Roberts density model, which rep- 
resents the upper portion of the atmosphere in the LAIRS 
Program. 
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2.1.1 ENTREE TRAJECTORY 


The LAIRS input trajectory data is a standard LaRC Best 
Estimate Trajectory (BET) file, output by the ENTREE Pro- 
gram. The BET file gives Shuttle altitude, latitude, and 
longitude, as well as the uncertainties in these parameters, 
as a function of time. 

2.1.2 METEOROLOGICAL DATA FILES 

The LAIRS Program can accept up to three meteorological data 
files, each containing a vertical profile of atmospheric 
parameters measured at a given station. These files may be 
created by a preprocessor (PREMET) from a meteorological 
data tape supplied by the National Oceanic and Atmospheric 
Administration (NOAA) . Table 2-1 shows the format of this 
tape; more information on the preprocessor and the files it 
creates is given in Appendix B. These files may also be 
created from raw data supplied to LaRC by NOAA, the Air 
Force, or other sources. 

2.1.3 KEYWORD CARD DESCRIPTIONS 

All of the LAIRS Program input, with the exception of the 
permanent, working, and meteorological files, is provided via 
keyword cards. Each keyword card includes an alphanumeric 
identifier in card columns 1-10. This identifier is called 
the keyword. This section describes each keyword card, its 
use, and the information that must be supplied on it. 

The keyword card listing is alphabetical, and the name of 
the keyword appears in the upper right-hand or left-hand 
corner of the page for quick reference. The keywords con- 
form to a format that is described on each keyword page. 

In all cases, this format consists of one alphanumeric field 
(AlO) and seven real fields (7F10.3). Data supplied in the 
seven real fields will be accepted by LAIRS as single- 
precision numerical data and must be expressed as a valid 

(text continues on page 2-23) 
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Table 2-1. 


Format of LAIRS Meteorological Data Tape 


(1 of 2)^ 


WORD 

SYMBOL 

DESCRIPTION 

UNITS 

1 

LAT 

Latitude 

Degrees, + N 

2 

LON 

Longitude 

Degrees, + E 
to 360 

3 

FLAG 

Data flag: 

= 0, measured data 
= 1, modeled data 
= 2 , combined measured 
and modeled data 


4 

— 

Spare 


5 

ALT 

Geometric altitude 

Feet 

6 

WS 

Horizontal wind speed 

Feet/second 

7 

WD 

Direction horizontal 
wind is coming from 
relative to true 
north, north being 
0 degrees, increas- 
ing positively 
clockwise 

Degrees 

8 

TE 

Ambient temperature 

Degrees C 

9 

PR 

Ambient pressure 

Millibars 

10 

D 

Ambient density 

2 

Grams/meter 

11 

DW 

Dew point (not used) 

Degrees C 

12 

TEU 

Ambient temperature 
systematic uncer- 
tainty 

Degrees C 

13 

PRU 

Ambient pressure 
systematic uncer- 
tainty 

Millibars 

14 

DU 

Ambient density 

systematic uncer- 
tainty 

Grams/meter 

15 

HWSUS 

Horizontal wind speed 
systematic uncer- 
tainty 

Feet/second 

16 

HWSUN 

Horizontal wind speed 
noise or fluctua- 
tion uncertainty 

Feet/second 


^All quantities are real variables that are one word in length. 
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Table 2-1. Format of LAIRS Meteorological Data Tape (2 of 2) 


— 

WORD 

SYMBOL 

DESCRIPTION 

UNITS 

17 

VWSUN 

Vertical wind speed 
noise or fluctua- 
tion uncertainty 

Feet/second 

18 

HWDUS 

Horizontal wind direc- 
tion systematic 
uncertainty 

Degrees 

19 

HWDUN 

Horizontal wind direc- 
tion noise or fluc- 
tuation uncertainty 

Degrees 

20 


Spare 

— 









CREATE 



Description 


CREATE - Keyword to create a 
new file in standard LAIRS 
meteorological file format from 
the output USE file. Default 
(no CREATE card) is to omit 
creation of this file. 

Blank 
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DC 


DC 


• Card format: AlO, 7F10.3 


• Detailed format; 
Columns Format 


1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 


AlO 

F10.3 
FIO. 3 

F10.3 

FIO. 3 

FlO. 3 
FIO. 3 


71-80 FIO. 3 


Description 


DC - Keyword to perform differ- 
ential corrections on atmos- 
pheric data 

Blank 

Use temperature data in differ- 
ential correction: 

= 0., yes (default) 

= 1. , no 

Use pressure data in differ- 
ential correction: 

= 0., yes (default) 

= 1. , no 

Use density data in differential 
correction : 

= 0., yes (default) 

= 1. , no 

Blank 

Polynomial segment number for 
this DC card:* 

= 0., this card applies to all 
segments (default) 

= 1. , top segment (65-90 kilo- 
meters ) 

= 2 ,, middle segment (25-65 kil- 
ometers) 

= 3., bottom segment (0-25 kil- 
ometers) 

Continuity with next segment 
imposed at upper segment bound- 
ary : 

= 0., yes (default) 

= 1. , no 


Up to three DC keyword cards (one for each segment) may be 
included in the deck. 
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DCEDIT 


DCEDIT 


• 

Card format: AlO, 

7F10.3 

• 

Detailed format: 


Columns 

Format 

Description 

1-10 

AlO 

DCEDIT - Keyword to set edit 
criteria for defferential 
corrections 

11-20 

F10.3 

Blank 

21-30 

F10.3 

RMS multiplier 
(default = 3.0) 

31-40 

F10.3 

Initial value of RMS 
(default = lO.ElO) 

41-50 

FIO . 3 

RMS adder parameter 
(default = 0.0) 

51-60 

FIO. 3 

A priori temperature uncertainty 
(default = 5 percent) 

61-70 

FIO. 3 

A priori pressure uncertainty 
(default = 5 percent) 

71-80 

FIO. 3 

A priori density uncertainty 


(default = 5 percent) 


NOTE: Any observation for which the absolute value of 

the weighted (0-C) is larger than the quantity 
(weighted RMS x RMS multiplier + RMS adder) is 
edited out. The weighted RMS used is obtained 
from the previous iteration (the initial value 
of RMS is used for the first iteration) . 
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DC ITER 


DCITER 

• Card format: AlO , 


• 

Detailed format 

Columns 

Format 

1-10 

AlO 

11-20 

FIO. 3 

21-30 

FIO. 3 

31-40 

FIO. 3 

41-80 

4F10. 3 


7F10.3 

Description 

DCITER - Keyword to set DC ite 
ration flags 

Blank 

Maximum number of iterations 
Convergence criterion 
Blank 
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DEBUG 


DEBUG 


• Card format: AlO, 7F10.3 

• Detailed format: 

Columns Format Description 

1“10 AlO DEBUG - Keyword to obtain debug 

printout 


11-20 

F10.3 

Blank 

21-30 

FIO. 3 

Debug 
= 1., 
= 2., 

= 3., 


= 4., 
= 5. , 


= 8 ., 

= 9., 

= 10. 
= 15. 


type : 

overall flow debug 
working file listing 
parameter calculation 
debug 

BET file listing 
created meteorological 
file listing (to be used 
with CREATE card) 
translational coefficient 
adjustment debug 
meteorological file listing 
(to be used in a polyno- 
mial model run) 

, differential correction 
debug 

, onboard measurement 
calculation debug 


31-80 5F10.3 Blank 
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DERIVE. 


DERIVE 

• Card format; AlO , 7F10.3 

• Detailed format: 

Columns Format Description 

1-10 AlO DERIVE - Keyword to override 

pressure and density calcula- 
tion with derivation of pres- 
sure and density from 
temperature via physical rela- 
tionships 


11-20 

FIO. 3 

Blank 

21-30 

FlO.3 

Integration step size (kilomete: 

31-40 

FlO.3 

Source of pressure boundary 
condition: 

= 0 . , this card 

= 1., default file 

= 2., meteorological profile 

41-50 

FlO.3 

Reference altitude (kilometers) 

51-60 

FlO.3 

Pressure at the reference al- 

2 

titude (newtons/meter ) 

61-80 

2F10 . 3 

Blank 
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END 


END 

• 

Columns 

1-10 

11-80 


Card format: AlO, 7F10.3 

Detailed format: 

Format Description 

AlO END - Keyword to specify end of 

keyword card deck 

7F10 . 3 Blank 
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ENTREE 


ENTREE 


• Card format: AlO , 7F10.3 

• Detailed format: 


Columns 

1-10 

11-20 

21-30 

31-40 

41-50 

51-70 


Format 

AlO 

FIO . 3 

Flo . 3 

FIO . 3 
FIO. 3 
2F10 .3 


Description 

ENTREE - Keyword to set ENTREE 
profile type (along Shuttle 
reentry path) 

YYMMDD. of first point on input 
BET file (date of Shuttle re- 
entry) 

Interval (in number of points) 
at which BET file will be 
sampled 

(default = 1., i.e., read each 
point) 

Lower altitude limit for BET 
file processing (kilometers) 
(default = bottom of BET file) 

Upper altitude limit for BET 
file processing (kilometers) 
(default == top of BET file) 

Blank 


71-80 


FIO .3 


Run identification number 



FIN 


FIN 


• Card format: AlO, 7F10.3 

• Detailed format: 


Columns Format 

I- 10 AlO 

II- 80 7F10.3 


Description 

FIN - Keyword following last 
stacked deck. This keyword 
terminates the run. 

Blank 
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INTJIRPOL 


INTERPOL 


• Card format: AlO, 7F10.3 

• Detailed format: 


Columns Format 


Description 


1-10 

AlO 

11-20 

FIO. 3 

21-30 

FIO. 3 

31-40 

FIO. 3 


41-80 4F10.3 


INTERPOL - Keyword to set in- 
terpolator type 

Blank 

Vertical interpolator type: 

= 1., first-order Lagrange 
= 2., second-order Lagrange 

Horizontal interpolator type: 
= 0., first-order interpola^ 
tion in horizontal distance 
from profile. 

= 1., first-order bivariate 
interpolation in latitude 
and local solar time 
= 2. f second-order bivariate 
interpolation in latitude 
and solar time 

Blank 
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LAP 


LAP 


• 

Card format: AlO, 

7F10.3 

• 

Detailed format: 


Columns 

Format 

Description 

1-10 

AlO 

LAP - Keyword to set Local 
Atmospheric Profile type 

11-20 

F10.3 

YYMMDD. for LAP 

21-30 

FIO .3 

HHMMSS.SS for LAP (GMT) 

31-40 

FIO . 3 

Latitude for LAP (decimal de- 
grees) 

41-50 

FIO. 3 

Longitude for LAP (decimal 
degrees) 

51-60 

FIO. 3 

Lower altitude for LAP (kilo- 
meters) 

61-70 

FIO. 3 

Upper altitude for LAP (kilo- 
meters) 

71-80 

FIO. 3 

Altitude interval for LAP 


(kilometers) 
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LOWMOD 


LOWMOD 


• 

Card format: AlO, 

7F10.3 

• 

Detailed format: 


Columns 

Format 

Description 

o 

I — 1 

1 

1 — 1 

AlO 

LOWMOD - Keyword to set model 
options for lower atmosphere 

11-20 

F10.3 

Blank 

21-30 

FIO .3 

Atmospheric model type; 


= 1., interpolate T, P, and D 
from default files 
(default) 


= 2., interpolate T, P, and D 
from meteorological pro- 
files 

= 3., calculate T, P, and D 

from default polynomial 
model 

=4,, calculate T, P, and D 

from adjusted polynomial model 

31-40 FlO.3 Wind model type;* 

= 0., none (default) 

= 1. , interpolate from default 
file 

= 2., interpolate from meteoro- 
logical profiles 
= 3 . , calculate from analytical 
model 


41-50 Flo. 3 


51-60 F10.3 


61-70 F10.3 


Include latitude variations; 

= 0 . , no (default) 

= 1. , yes 

Include diurnal and semidiurnal 
variations : 

= 0., no (default) 

= 1., yes 

Include weighting in adjusted 
polynomial model; 

= 0. , no (default) 

= 1 . , yes 


'ic 

wind model choice is applicable only to atmospheric model 
types 3 and 4 . Types 1 and 2 automatically return inter- 
polated wind values (unless wind type 0 is chosen) . 
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LOWMOD 


LOWMOD 

Columns 

71-80 


(Cont ' d) 


Format Description 

Flo. 3 Recalculate final temperature 

to agree with gas law: 

== 0., yes (default) 

= 1. , no 
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MODBQUND 


MODBOUND 

• Card format: AlO, 7F10.3 

• Detailed format: 


Columns 

Format 

Description 

1-10 

AlO 

MODBOUND - Keyword to set 
boundaries of atmospheric model 
segments 

11-20 

FIO. 3 

Blank 

21-30 

F10.3 

Boundary between first and sec- 
ond segments (kilometers) 
(default = 90 kilometers) 

31-40 

Flo. 3 

Boundary between second and 
third segments (kilometers) 
(default = 65 kilometers) 

41-50 

FIO. 3 

Boundary between third and 
fourth segments (kilometers) 
(default = 25 kilometers) 

51-80 

3F10. 3 

Blank 
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MODREF 


MOD REF 


• Card format: AlO, 7F10.3 


• 

Detailed format: 


Coliimns 

Format 

Description 

1-10 

AlO 

MODREF - Keyword to set model 
reference values 

11-20 

FlO. 3 

Blank 

21-30 

Flo. 3 

Model reference latitude (de- 
grees) 

31-40 

FlO. 3 

Model reference local solar 
time (decimal hours) 

41-80 

4F10. 3 

Blank 
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PRINTALL 


PRINTALL 

• Card format: AlO, 

• Detailed format: 

Columns Format 

1-10 AlO 


11-80 7F10.3 



Description 

PRINTALL - Keyword to print all 
polynomial fitting residual 
reports in addition to the 
standard reports. Default 
(no PRINTALL card) is to print 
only the standard reports . 

Blank 
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TRANADJ 


TRAN AD J 


• Card format: AlO , 7F10.3 

• Detailed format: 


Columns 

1-10 

11-20 


Format 

AlO 

FIO . 3 


Description 


TRANADJ - Keyword to invoke 
adjustment of translational 
coefficients 

Blank 


21-30 Flo, 3 


31-40 FIO. 3 


41-80 410.3 


Diurnal/semidiurnal coefficient 
adj ustment 
= 0., no (default) 

= 1 . , yes 

Latitude gradient adjustment 
= 0., no (default) 

= 1. , yes 

Blank 


2-21 



UPMOD 


UPMOD 


• 

Card format: AlO, 

7F10.3 


• 

Detailed format: 



Columns 

Format 


Description 

1-10 

AlO 

UPMOD - 
options 

Keyword to set model 
for upper atmosphere 

11-20 

FIO. 3 

Blank 


21-30 

FIO. 3 

Atmospheric model type: 


= 1., interpolate T, P, and D 
from default files 
(default) 

=2., interpolate T, P, and D 
from meteorological pro- 
files 

= 3., calculate T, P, and D 
from default Jacchia- 
Roberts model 

= 4., calculate T, P, and D 
from adjusted Jacchia- 
Roberts model 

31-40 FIO . 3 Wind model type* 

= 0., none (default) 

= 1. , interpolate from default 
file 

= 2 . , interpolate from meteoro- 
logical profiles 

41-50 FIO. 3 Include latitude variations** 

= 0., no (default) 

= 1 . , yes 

51-60 FIO, 3 Include diurnal and semidiurnal** 

variations 
= 0., no (default) 

= 1 . , yes 

61-80 2F10.3 Blank 

*Wind model choice is applicable to atmospheric model types 
3 and 4 only. Types 1 and 2 automatically return inter- 
polated wind. 

**Latitude and diurnal/semidiurnal variations are auto- 
matically included in the Jacchia-Roberts atmospheric 
models, types 2 and 4. 
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floating point representation (i.e., either a number punched 
within the field containing a decimal point or a number 
punched with an E decimal exponent right-justified in the 
field) . Any data value supplied for the real fields (even 
though the data may be a whole number) must be a floating 
point representation. 

The keyword descriptions state the default data values where 
appropriate. A default value will be used if the data field 
is left blank or if a value of zero is supplied in the field. 
Any value in the field other than zero or blank will override 
a specified default. 

The data units used (unless otherwise specified) are meters, 
seconds, meters/second, and degrees. Units for most quanti- 
ties are standard M.K.S. units. Dates supplied on LAIRS 
input cards are usually words in the packed form: 

yymmdd 

hhmmss . ss 

2.2 OUTPUT REPORTS 

The output reports for the LAIRS Program consist of several 
separate reports. These include the LAIRS Initialization 
Reports, the File Statistics Report, the Polynomial Fitting 
Reports, the Atmospheric Parameters Output Report, and the 
Summary Report. Each report will be explained in more 
detail later. Briefly, the output reports perform the 
following functions; 

1. Initialization Reports - The first of these reports 
shows the input keyword card images. This report 
is followed by the Initial Conditions Report for 
the lower and upper atmosphere models, which in- 
cludes a description of the output parameters. 
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2. File Statistics Report - This report summarizes 
the features of the working file or the meteorolog- 
ical files that are the source of data for this 
LAIRS run. 

3. Polynomial Fitting Reports - These reports detail 
aspects of the polynomial fitting (both linear and 
nonlinear) such as observations, residuals, weight- 
ing factors, and calculated coefficients. 

4. Atmospheric Parameters Output Report - This report 
consists of a block printout in tabular form of 
all the output parameters for each trajectory point 
as a function of altitude. 

5. Error Summary Report - This report lists any errors 
that were encountered during the LAIRS run. 

2.2.1 INITIALIZATION REPORTS 

The first page of the printed output lists the input data 
keyword cards for verification (see Figure 2-1) . 

The second page of output is the Initial Conditions Report 
(Figure 2-2) . This report contains a description of the 
run as specified by the control parameters. The type of pro- 
file is reported along with the epoch of the run. The lower 
and upper atmosphere parameters are specified, including the 
type of atmosphere and wind models, the variations desired, 
and the interpolation type. 

The last pages of the report contain a description of the 
output parameters (see Figures 2-3 and 2-4) . This includes 
the title, the units, and a description of each parameter 
that will be output in the Atmospheric Parameters Output 
Report . 
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LAIRS KEYWORD CARD IMAGES 


ENTREE 

810414*0000 

10.0000 

60.0000 

95.0000 

0.0000 

0 * 000 < 

UPMOD 

0*0000 

3*0000 

2.0000 

1.0000 

1.0000 

0*000< 

LOWMOD 

0*0000 

4.0000 

2.0000 

1 .0000 

1.0000 

0*000( 

MODBOUND 

0*0000 

90*0000 

68.0000 

14.0000 

4.0000 

0*000( 

INTERPOL 

0*0000 

2.0000 

0.0000 

0.0000 

0.0000 

0.000< 

DC 

0*0000 

0.0000 

0.0000 

0.0000 

0.0000 

0 . 000 < 

PRINTALL 

0*0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.000( 

STATION 

0*0000 

17.0830 

0.0000 

0.0000 

0.0000 

0*Q00< 

STATION 

0*0000 

17.2167 

0.0000 

0.0000 

0.0000 

0*000< 

STATION 

0.0000 

18.3670 

0.0000 

0.0000 

0.0000 

0 . 000 < 

END 

0*0000 

0.0000 

0.0000 

0.0000 

0.0000 

0*000< 


Figure 2-1. LAIRS Keyword Card Images 



LANGLEY ATMOSPHERIC INFORMATION RETRIEVAL SYSTEM 


INITIAL CONDITIONS REPORT 

EPOCH =810414 174230»000 PROFILE IS 

LOWER ATMOSPHERE MODEL 

TYPE = ADJUSTED POLYNOMIAL 

WIND MODEL = METEOROLOGICAL INTERPOL ♦ 

LAT. VARIATIONS INCLUDED ~ YES 

D/SD VARIATIONS INCLUDED - YES 


ENTREE ID NUMBER = 0 

UPPER ATMOSPHERE MODEL 
TYPE = JACCHIA-ROBERTS 
WIND MODEL = METEOROLOGICAL 
LAT. VARIATIONS INCLUDED “ 
D/SD VARIATIONS INCLUDED - 

2ND ORDER LAGRANGE 

- 1ST ORDER UNIVARIATE 
YES 

- NO 


ALTITUDE INTERPOLATION USED - 
METEOROLOGICAL INTERPOLATION USED 
DIFFERENTIAL CORRECTION APPLIED - 


DERIVE PRESSURE AND DENSITY FROM TEMPERATURE 
RECOMPUTE TEMPERATURE FROM GAS LAW “ NO 


INTERPOL* 

YES 

YES 


Figure 2-2. Initial Conditions Report 
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DESCRIPTION OF OUTPUT PARAMETERS - PART I 


TITLE 

UNITS 

DESCRIPTION 



ALT 

KM 

ALTITUDE 



LAT 

DECIMAL DEG 

LATITUDE 



LONG 

DECIMAL DEG 

LONGITUDE 



GMTSEC 

DECIMAL SEC (TIME) 

GMT IN SECONDS 


LST 

DECIMAL HOURS 

LOCAL SOLAR 

TIME 

TEMP 

DEG KELVIN 

ATMOSPHERIC 

TEMPERATURE 

PRESSURE 

N/M**2 

ATMOSPHERIC 

PRESSURE 

DENSITY 

KG/M**3' 

ATMOSPHERIC 

DENSITY 

WNDSP 

M/SEC 

WIND SPEED 



UNDDR 

DEG FROM N 

WIND DIRECTION 


TEMPU 

DEG KELVIN 

UNCERTAINTY 

IN 

TEMPERATURE 

PRESU 

N/M**2 

UNCERTAINTY 

IN 

PRESSURE 

DENU 

KG/M*))t3 

UNCERTAINTY 

IN 

DENSITY 

USU 

M/SEC 

UNCERTAINTY 

IN 

WIND SPEED 

WDU 

DECIMAL DEG 

UNCERTAINTY 

IN 

WIND DIRECTION 


Figure 2-3. Description of Output Parameters (Part 1) 
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DESCRIPTION OF OUTPUT PARAMETERS 


PART II 


TITLE UNITS DESCRIPTION 

ALT KM ALTITUDE 

6MTSEC DECIMAL SEC* (TIME) 6MT IN SECONDS 


M*WT 

GM/MOLE 


SC.HT 

KM 


MACH NO 

SPEED OF 

SOUND 

PN 

N/M**2 


PCDEM 

PERCENT 


E-W WD 

M/SEC 


N-S WD 

M/SEC 


ALTU 

KM 


LATU 

DECIMAL 

DEG 

LONU 

DECIMAL 

DEG 

E-W U 

M/SEC 


N-S U 

M/SEC 



MEAN MOLECULAR WEIGHT 
PRESSURE SCALE HEIGHT 
MELOCITY(MACH NUMBER) 

ONBOARD PRESSURE 
PERCENT DEVIATION FROM GAS LAW 
EAST-WEST WIND COMPONENT 
NORTH-SOUTH WIND COMPONENT 
UNCERTAINTY IN ALIITUDE 
UNCERTAINTY IN LATITUDE 
UNCERTAINLY IN LONGITUDE 
EAST-WEST WIND UNCERTAINTIES 
NORTH-SOUTH WIND UNCERTAINTIES 


Figure 2-4. Description of Output Parameters {Part 2) 



2.2.2 FILE STATISTICS REPORT 


This report shows the latitude and longitude of each meteor- 
ological station supplying data for the run, as well as the 
number of observations and the altitude limits on the files 
provided by each station (see Figure 2-5) . Also shown are 
the atmospheric segment boundaries and the temperature and 
density at the top boundary, which are used by the Jacchia- 
Roberts model and by the differential correction polynomial 
model. The value ZREF, which is the reference altitude sup- 
plied by the user for the (optional) pressure and density 
override integration from temperature, is given, together 
with the pressure at that altitude. 

If the default model is used for atmospheric parameter calcu- 
lation, the above information will be taken from the default 
file . 

2.2.3 POLYNOMIAL FITTING REPORTS 

The exact nature of the polynomial fitting reports will de- 
pend on whether the user has selected a simple, linear fit 
to the meteorological data or has selected the nonlinear 
differential correction process. As the latter possibility 
includes all the reports that would appear in the former, 
the reports illustrated here are those that a user would see 
when selecting the differential correction option. 

The first report (see Figure 2-6) shows the result of the 
initial fit to temperature data which establishes the default 
coefficients for the differential correction process (this 
report is printed segment by segment; only the report for 
the upper segment is shown here) . The temperature observa- 
tions for each altitude point on each file are tabulated, 
along with the computed value at that point and the residual. 
Because no uncertainties were provided for this sample data. 
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FILE STATISTICS REPORT 

SOURCE OF data: meteorological files 

STATION NO. 1 LONGITUDE 200.20 LATITUDE 22.00 

31 POINTS BETWEEN 58.00 KM AND 92.00 KM 
STATION NO. 2 LONGITUDE 240.90 LATITUDE 34.10 

35 POINIS BETWEEN 23.32 KM AND 67.25 KM 
STATION NO. 3 LONGITUDE 242.27 LATITUDE 34.98 

100 POINTS BETWEEN .72 KM AND 30.78 KM 

ATMOSPHERIC SEGMENT BOUNDARIES (KM) 1 BOUNDl = 90.00 B0UND2 = 68.00 BOUNDS = 14.00 

T(BOUNDl) = 191.000 DEGREES 

D(BOUNDl) = .326E-05 K6/Mit:)K3 

ZREF = 30.000 KM 

P(ZREF) = .113E+04 N/M**2 


Figure 2-5. File Statistics Report 



LAIRS POLYNOMIAL FITTING REPORT 


ITERATION NUMBER 0 


UPPER SEGMENT 


ALTITUDE 

GBS 

OBSERVED 

COMPUTED 

RESIDUAL 

UNCERTAINTY 

OBSERVED 

COMPUTED 

WEIGHTED 

WEIGHT 

EDIT 

(KM) 

TYPE 

0 

C 

□ -C 

DEL 0 

LOG 0 

LOG C 

RESIDUAL 


FLAG 

92.00 

TEMP 

189.0 

189.6 

-.6060 

0. 



-.6413E-01 

, 1058 

1 

91.00 

TEMP 

191,0 

189.1 

1,893 

0. 



,1982 

.1047 

1 

90.00 

TEMP 

191.0 

188.8 

2.236 

0. 



.2342 

. 1047 

1 

85.00 

TEMP 

189.2 

189.5 

- .2895 

0. 



-.3063E-01 

. 1058 

1 

84,00 

TEMP 

185.1 

190.1 

-5.040 

0 . 



-.5479 

.1087 

1 

83.00 

TEMP 

185.0 

190.9 

-5.917 

0. 



— . 6466 

,1093 

1 

82.00 

TEMP 

188,0 

191,9 

-3,934 

0, 



-.4252 

, 1081 

1 

81,00 

TEMP 

193.0 

193 . 1 

-.1067 

0, 



-. 1129E-01 

.1058 

1 

80.00 

TEMP 

197 . 0 

194.4 

2.547 

0. 



.2653 

.1042 

1 

79,00 

TEMP 

199.0 

196.0 

3.023 

0. 



,3133 

, 1036 

1 

78.00 

TEMP 

201 . 0 

197,7 

3.317 

0. 



.3419 

.1031 

1 

77.00 

TEMP 

205.0 

199,6 

5.424 

0. 



.5506 

. 1015 

1 

76.00 

TEMP 

206.0 

201.7 

4.340 

0 . 



,4406 

, 1015 

1 

75,00 

TEMP 

208,0 

203.9 

4.061 

0. 



.4102 

.1010 

1 

74,00 

TEMP 

209,3 

206.4 

2,931 

0, 



.2945 

,1005 

1 

73.00 

TEMP 

210,7 

209.1 

1,601 

0. 



.1601 

. lOOOE+00 

1 

72.00 

TEMP 

210.0 

212,0 

-1 .930 

0. 



-.1940 

. 1005 

1 

71 ,00 

TEMP 

212.4 

215.0 

-2.668 

0. 



-.2654 

.9950E-01 

1 

70.00 

TEMP 

213.7 

218.3 

-4,614 

0. 



-.4568 

.9901E-01 

1 

69.00 

TEMP 

215.0 

221,8 

-6.852 

0. 



-.6751 

.9852E-01 

1 

68.00 

TEMP 

216,2 

225.5 

-9.306 

0. 



- .9123 

.9804E-01 

1 

67.00 

TEMP 

221.4 

229.4 

-7.979 

0 . 



-.7635 

.9569E-01 

1 

66.00 

TEMP 

232,7 

233.5 

-.8732 

0. 



-.7938E-01 

.9091E-01 

1 

65 . 00 

TEMP 

244r9 

237.9 

7.007 

0. 



.6040 

,8621E-01 

1 

67.25 

TEMP 

242.8 

228.4 

14,34 

0. 



1 .199 

.8359E-01 

1 

64.53 

TEMP 

245.0 

240.0 

4,954 

0, 



,4117 

.8310E-01 

1 


WEIGHTED RMS FOR THIS ITERATION= .49909279 
TOTAL NUMBER OF POINTS = 26 


Figure 2-6 . LAIRS Polynomial Fitting Report 



an a priori uncertainty of 5 percent (which can be changed 
by the user) was used for all parameters. It is necessary 
to provide uncertainties in order to balance the contribu- 
tions that each of the parameters makes to the fit. As no 
logarithms are used in the initial temperature fit, the 
colximns for LOG O and LOG C are blank. The edit flag indi- 
cates whether an observation has been edited out. An edit 
flag of 1 means that the observation has been used; a flag 
of zero means that it has been edited. At the bottom of the 
report are the weighted root mean square of the residuals 
and the total number of points used in the fit. 

The second residual report (Figure 2-7) is output during 
the differential correction process. It is identical in 
format to the first report, except that, for pressure and 
density, the logarithms of the observed and computed values 
(which are the actual subjects of the fitting process) are 
displayed. The weights and weighted residuals are also ex- 
pressed in log space. 

The report shown in Figure 2-8 follows these residual reports 
for each segment and for each iteration. The number of data 
points available in this segment and the number of points 
used (i.e. , not edited) are given, as are the values of the 
current, previous, and predicted root mean squares of the 
residuals. A summary of the current status of the solve- 
for parameters is also presented. 

In a differential correction run, the variance-covariance 
matrix of the solve-for parameters is computed in each 
iteration and printed out at the end of the iteration (see 
Figure 2-9) . The diagonal elements of this matrix define 
the variances of the solve-for parameters, and the off- 
diagonal elements are used to compute correlation coefficients 
between two different solve-for parameters. The square root 
of the variance of a solve-for parameter gives the standard 
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LAIRS NON LINEAR LEAST SQUARES ESTIMATION REPORT 


ITERATION NUMBER 3 


UPPER SEGMENT 


ALTITUDE 

DBS 

OBSERVED 

COMPUTED 

RESIDUAL 

UNCERTAINTY 

OBSERVED 

COMPUTED 

WEIGHTED 

WEIGHT 

EDIT 

(KM) 

TYPE 

0 

C 

0-C 

DEL 0 

LOG 0 

LOG C 

RESIDUAL 


FLAG 

92.00 

PRES 

.1270 

.1315 

-.4512E-02 

0. 

-2.064 

-2,029 

-.7155 

20.50 

1 

91.00 

PRES 

.1510 

.1567 

-.5739E-02 

0. 

-1.890 

-1.853 

-.7646 

20.50 

1 

90.00 

PRES 

. 1790 

,1872 

-.8219E-02 

0. 

-1.720 

-1.675 

-.9201 

20.50 

1 

85.00 

PRES 

.4783 

.4638 

. 1449E-01 

0. 

-.7375 

-.7683 

.6307 

20,50 

1 

84,00 

PRES 

.5802 

.5566 

.2357E-01 

0. 

-.5444 

-.5859 

.8502 

20,50 

1 

83,00 

PRES 

.6993 

,6676 

.3168E-01 

0. 

-.3577 

-.4041 

.9501 

20.50 

1 

82.00 

PRES 

.8383 

.8000 

.3830E-01 

0. 

-.1764 

-.2231 

.9585 

20.50 

1 

81,00 

PRES 

.9990 

.9576 

.4144E-01 

0. 

-.1008E-02 

-.4337E-01 

.8683 

20.50 

1 

80.00 

PRES 

1.183 

1.145 

.3853E-01 

0, 

.1681 

.1350 

,6787 

20.50 

1 

79.00 

PRES 

1.406 

1.366 

.4020E-01 

0. 

.3407 

.3117 

.5946 

20.50 

1 

78.00 

PRES 

1 .662 

1 .627 

.3539E-01 

0, 

.5080 

.4865 

.4411 

20.50 

1 

77.00 

PRES 

1.958 

1.934 

.2479E-01 

0, 

.6721 

.6594 

.2612 

20.50 

1 

76.00 

PRES 

2.300 

2.293 

.6337E-02 

0. 

.8328 

.8301 

.5655E-01 

20.50 

1 

75.00 

PRES 

2.697 

2.714 

- . 1696E-01 

0. 

.9923 

.9985 

-.1285 

20.50 

1 

74.00 

PRES 

3.168 

3.205 

-.3723E-01 

0. 

1 . 153 

1.165 

-.2395 

20.50 

1 

73,00 

PRES 

3.716 

3.775 

-.5921E-01 

0. 

1.313 

1.328 

-.3240 

20.50 

1 

72.00 

PRES 

4.351 

4.436 

-.8509E-01 

0, 

1.470 

1.490 

-.3969 

20,50 

1 

71.00 

PRES 

5.097 

5.201 

-.1044 

0, 

1.629 

1.649 

-.4154 

20,50 

1 

70.00 

PRES 

5.939 

6.083 

-.1440 

0. 

1.781 

1,805 

-.4910 

20.50 

1 

69,00 

PRES 

6.998 

7.097 

-.9958E-01 

0. 

1.946 

1.960 

-.2896 

20.50 

1 

68.00 

PRES 

8.190 

8.262 

-.7255E-01 

0. 

2.103 

2.112 

-.1808 

20.50 

1 

67.00 

PRES 

9.547 

9.597 

-.4942E-01 

0, 

2.256 

2.261 

-.1058 

20,50 

1 

66.00 

PRES 

11,09 

11.12 

-.2803E-01 

0, 

2.406 

2.409 

-.5173E-01 

20.50 

1 

65,00 

PRES 

12.79 

12.86 

-.7782E-01 

0, 

2.548 

2.554 

-.1244 

20.50 

1 

67.25 

PRES 

9,179 

9,243 

-.6419E-01 

0. 

2.217 

2.224 

-.1428 

20.50 

1 

64.53 

PRES 

13.30 

13.77 

-.4683 

0, 

2,588 

2,622 

-.7092 

20.50 

1 


Figure 2-7. LAIRS Nonlinear Least Squares Estimation Report (Part 1) 
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LAIRB NON LINEAR LEAST SQUARES ESTIMATION REPORT 


ITERATION NUMBER 3 

UPPER SEGMENT 


DC HAS CONVERGED 





AVAILABLE 

USED 

NUMBER OF 

TEMPERATURE DATA 

26 

26 

NUMBER OF 

DENSITY DATA 

26 

26 

NUMBER OF 

PRESSURE 

DATA 

26 

26 

CURRENT 

WEIGHTED 

RMS 

.69619249 


PREVIOUS 

WEIGHTED 

RMS 

.69753787 


PREDICTED 

WEIGHTED 

RMS 

.69619348 



TEMP WRMS = .72058210 PRES WRMS = ,55948992 DENS WRMS = ,78853309 


SOLVE-FOR PARAMETER REPORT 

APRIORI VALUE CURRENT VALUE CORRECTION 

CURRENT-PREVIOUS 

TEMP COEFF AO 
TEMP COEFF A1 
TEMP COEFF A2 
TEMP COEFF A3 
LOG PRES 


188,76371 185.91303 -.20062635 

-.26482847 -1.8133701 , 36539087E-01 

.78974132E-01 .23398251 - . 85495997E-03 

.40869166E-03 - . 35267601E-02 - . 9454 1438E-05 

-1.5906200 -1.6754745 - . 23142881E-02 


Figure 2-8. LAIRS Nonlinear Least Squares Estimation Report (Part 2) 



i 


COVARIANCE HATRIX AND CORRELATION COEFFICIENTS 


A1 


A2 


A3 


AO LOG PRES 


Isj 

I 

LJ 

U1 


Ai 

.168583E-I^01 

-.135638E+00 

♦313935E-02 

-.539155E+01 

A2 

-.961177E+00 

.118125E-01 

-.286776E-03 

. 36561 lElOO 

A3 

♦902397E+00 

-.984773E+00 

.717912E-05 

-.762422E-02 

AO 

-.849936E+00 

.688537E+00 

-.582424E+00 

.238694E+02 

LOG PRES 

-.249686E+00 

.739109E-01 

-.513650E-03 

.667027E+00 


-.A83610E-02 
.169390E-03 
-.290209E-07 
.A87183E-01 
♦ 4>»4648E“03 


TEMPERATURE PRESSURE 


DENSITY 


UEIGHTED ROOT MEAN SQUARE 
AVERAGE UEIGHTED RESIDUALS 
STANDARD DEV* OF UEIGHTED RESIDUALS 


.62743467 

.12560247 

.61473432 


1.0350223 

-.10205341 

1.0299788 


.98876449 

.10205172 

.98348394 


! 

Figure 2-9. Covariance Matrix and Correlation Coefficients Report 


deviation of the solve- for parameter. The upper triangular 
part (including the diagonal elements) of the matrix printed 
out represents the variance-covariance matrix and the lower 
triangular part represents the correlation coefficients. 

When the differential correction process is completed, an 
Estimation Summary Report (see Figure 2-10) is produced. 

This report displays the calculated temperature coefficients 
for each segment and the boundary value of pressure at the 
top of each segment. 

2.2.4 ATMOSPHERIC PARAMETERS OUTPUT REPORT 

This report, which is divided into two parts, provides tables 
of all the output parameters for each trajectory point (see 
Figures 2-11 and 2-12) . In addition to the given (ENTREE) 
information regarding altitude, latitude, longitude, and 
time (GMT seconds) , Part I of the report prints the value 
of local solar time (solar hour angle, measured from zero 
hours at noon) and the calculated values for temperature, 
pressure, density, wind speed, and wind direction. Also 
printed are the uncertainties in all the calculated param- 
eters except local solar time. Part II of the report lists 
the mean molecular weight used in calculating the parameters 
and the calculated values of pressure scale height, Mach 
number, onboard pressure measurement, and the east-west and 
north-south wind components. The percent deviation from the 
gas law of the calculated values of temperature, pressure, 
and density is reported, although this deviation is removed 
in the final reported parameters (the percent deviation 
serves as a lower estimate on the error inherent in the 
LAIRS processing) . Part II of the report also gives the 
uncertainty in the input trajectory point and in the Cartesian 
wind components and repeats the altitude and time values as 
an index . 
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LAIRS NON LINEAR LEAST SQUARES 
ESTIMATION SUMMARY REPORT 
UPPER MIDDLE 






SEGMENT 

SEGMENT 


TEMP 

COEFF 

AO 

105.90976 

221 .52518 


TEMP 

COEFF 

A1 

-1.8151987 

6.6014731 

to 

1 

TEMP 

COEFF 

A2 

.23430978 

-.26004440 

■37 

TEMP 

COEFF 

A3 

-.35552341E-02 

.24965560E 


PRESSURE AT TOP BOUNDARY -1.6754178 2,1125525 


JACCHIA-ROBERTS BOUNDARY VALUES READJUSTED TO MATCH POLYNOMIAL MODEL? 
NEW TEMPERATURE = 105.910 DEGREES KELVIN 
NEW DENSITY = .304E-05 KG/M**3 


Figure 2-10. Estimation Summary Report 
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ATMOSPHERIC PARAMETERS OUTPUT REPORT J PART II 


ALT 

6HTSEC 

H.UT 

SC.HT 

MACH NO 

PN 

PCDEi .^ 

E-U UD 

N-S WD 

ALTU 

LATU 

L0N6U 

E-W U 

N-S U 

93.47 

64337.00 

28.4480 

5,735 

27.43 

. 79828E+02 

1.82 

- 9999.00 

- 9999.00 

.09 

.00 

.00 

0.00 

0.00 

92.07 

64347.00 

28.6261 

5.676 

27.48 

. 10151E+03 

1.18 

- 9999.00 

- 9999.00 

.09 

.00 

.00 

0.00 

0.00 

90.70 

64357.00 

28.7695 

5.637 

27.50 

, 13156E - f03 

.68 

44.21 

- 5.55 

.09 

.00 

.00 

0.00 

0.00 

89.34 

64367.00 

28.9644 

5.503 

27.73 

. 16447E - f03 

-.95 

20.19 

3.25 

,09 

.00 

.00 

0.00 

0.00 

88.01 

64377,00 

28.9644 

5.660 

27.34 

. 22341E+03 

3.09 

1.20 

11.81 

.09 

.00 

.00 

0.00 

0.00 

86.71 

64387.00 

28.9644 

5.826 

26.93 

. 29101E+03 

6.96 

- 18.22 

24.51 

.09 

.00 

.00 

0.00 

0.00 

85 . 4S 

64397.00 

28.9644 

5.500 

27.71 

. 32917E+03 

1.42 

- 26.69 

29,31 

.09 

.00 

.00 

0.00 

0.00 

84.23 

64407.00 

28.9644 

5.288 

28.24 

. 39982E+03 

- 2.27 

- 29.14 

31.43 

.08 

.00 

.00 

0.00 

0.00 

83.07 

64417.00 

28.9644 

5.297 

28.19 

, 46877E+03 

- 2.05 

- 29.19 

29.39 

.08 

.00 

.00 

0.00 

0.00 

81.99 

64427,00 

28.9644 

5.332 

28.07 

. 56940E+03 

- 1.60 

- 28,14 

24.39 

.08 

.00 

.00 

0,00 

0.00 

80.99 

64437.00 

28.9644 

5,382 

27.90 

. 67577E+03 

- 1.03 

- 27,35 

16.98 

.08 

.00 

.00 

0.00 

0.00 

80.08 

64447.00 

28.9644 

5.438 

27.72 

. 76711E+03 

-.45 

- 26.59 

8.86 

.08 

.00 

.00 

0.00 

0,00 

79,28 

64457.00 

28.9644 

5.505 

27.50 

. 90438E403 

.14 

- 25.41 

1.04 

.08 

.00 

.00 

0.00 

0.00 

78.61 

64467.00 

28.9644 

5.568 

27,30 

. 98323E+03 

.65 

- 24.08 

- 5,59 

.08 
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Figure 2-12, Atmospheric Parameters Output Report (Part 2) 



nil 


2.2.5 ERROR SUMMARY REPORT 

The Error Summary Report is the last report generated. If 
any errors occurred during the execution of LAIRS, a message 
identifying these errors is printed in the Error Summary 
Report . 
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SECTION 3 


SYSTEM DESCRIPTION 


This section describes the requirements for LAIRS and the 
software developed to fulfill these requirements. 

3 . 1 SOFTWARE REQUIREMENTS 

The LAIRS software must perform the following functions 
(Reference 1) : 

• Respond to user input specifying program options 

• Receive as input an ENTREE Program trajectory 
file containing Shuttle position vectors as a 
function of time 

• Receive as input meteorological observation files 
containing vertical profiles of temperature, 
pressure, density, and winds taken at up to three 
locations along or near the Shuttle reentry path 

• Produce values of the temperature, pressure, 
density, and winds for selected points along 
the Shuttle trajectory using the input meteoro- 
logical data or default data 

• Produce output reports and an output user file 
containing the values of the calculated parameters 

3 . 2 DESCRIPTION OF LAIRS MODELING CAPABILITIES 

LAIRS is designed for flexibility. A variety of parameter 
calculation modes are available, and the user may choose 
the one or ones best suited to his needs. All LAIRS model 
types are designed to operate in two basic modes: default 

and adjusted. In the default mode, the source of data used 
by the model is the Edwards Reference Atmosphere (0-25 kilo- 
meters) (Reference 2) , updated to remove discrepancies, and 
the Cospar International Reference Atmosphere (CIRA) 1972 
profiles (25-110 kilometers) (Reference 3) . This data is 
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stored on a permanent file (see Appendix B) , which also 
contains latitude gradients and diurnal and semidiurnal vari- 
ation coefficients for temperature, pressure, and density. 

These values were obtained from Reference 3 except for the 
latitude gradients below 25 kilometers, which were obtained 
from the Handbook of Geophysics (1965) (Reference 4) . In 
the adjusted mode, the meteorological observations described 
in Appendix B are used as input to the model. 

The models used in LAIRS are : 

• Interpolation of atmospheric parameters at points 

on the Shuttle trajectory. This is the most straight- 
forward model available in LAIRS. 

When more than one meteorological profile is avail- 
able, the user may choose between two interpola- 
tion schemes. The first consists of univariate 

* 

interpolation in altitude on each profile, followed 
by univariate interpolation in distance between 
profiles. The second, and more sophisticated, 
scheme involves univariate interpolation in altitude, 
followed by bivariate interpolation in latitude 
and local solar time. This method requires the 
use of latitude gradients . 

When only one profile is available, as in the de- 
fault case or when only one meteorological observa- 
tions file is available, the interpolation model 
consists merely of a univariate interpolation in 
altitude followed by (at the user's request) trans- 
lation of the interpolated values in latitude and 
local solar time to the Shuttle location. This 
translation makes use of the diurnal and semidiurnal 


* 

All altitude interpolation is Lagrangian for the temperature 
and exponential for the pressure and density. 
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coefficients and the latitude gradients stored 
on the permanent file. 

• The Jacchia-Roberts model (for use only above 90 
kilometers) . This model returns values of tempera- 
ture, pressure, and density for any input point 
above 90 kilometers (wind components are not pro- 
duced by the Jacchia-Roberts model) . The only data 
that must be input to the model are the geomagnetic 

index (K ) and solar flux Ft,. - values (stored on a 
p 10.7 

permanent file) , together with boundary values of 
temperature and density at 90 kilometers. These 
boundary values are chosen to match the model being 
used below 90 kilometers and are thus taken from 
either the default CIRA profiles or the meteoro- 
logical profiles. 

Diurnal and latitude effects are incorporated in the 
Jacchia-Roberts model. The coefficients and grad- 
ients stored on the permanent file are not required. 
No uncertainties are estimated by the model. 

• The fitted polynomial model (for use only below 90 
kilometers). In this scheme, the atmosphere is 
broken into three segments (0-25 kilometers, 25-65 
kilometers, and 65-90 kilometers) , and the data in 
each segment are fitted separately to polynomials. 
Winds are not treated in this model, but may be in- 
terpolated from the meteorological data. The poly- 
nomial coefficients produced by the model correspond 
to a single vertical profile at a reference latitude 
and reference local solar time. At each altitude on 
the Shuttle trajectory, the atmospheric parameters 
are calculated from the polynomial coefficients and 
translated to the Shuttle latitude and local solar 
time by means of the diurnal and semidiurnal coeffi- 
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cients and the latitude gradients. Uncertainty 
estimates are derived from the residual sum of 
squares of the polynomial fit. 

When the default profile is used, or when only one 
meteorological file is available, the model reference 
latitude and local solar time are taken to be those 
of the profile. Then the profile is sent to the 
polynomial fitting routines for calculation of the 
polynomial coefficients. If more than one meteoro- 
logical file is to be used, the data must first be 
translated to a single reference latitude and local 
solar time by utilizing diurnal/semidiurnal coeffi- 
cients and latitude gradients. The single resulting 
reference profile is then sent to the polynomial 
fitting routines. 

The polynomial fitting itself can take one of two 
forms at the user's request. The first is simply 
linear least squares fitting of the temperature and 
the logarithms of pressure and density. The second 
is the differential correction (DC) process dis- 
cussed in Section 4 . The linear least squares pro- 
cess requires as input a reference profile containing 
temperature, pressure, and density measurements, 
while the DC process can use as input a reference 
profile containing any one of these measurement types 
alone or any combination of tnem. 

A related aspect of the polynomial model is the 
integration of temperature to yield pressure and 
density after the polynomial coefficients have been 
calculated. This is an option available when the 
linear least squares process is used, but occurs 
automatically when the DC process is used, as the 
DC process produces only temperature coefficients. 
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• Adjustment of diurnal/semidiurnal coefficients and 

latitude gradients. This is not actually a complete 
model, but is an option the user may wish to select 
in addition to one of the preceding models. It 
consists of the utilization of meteorological data 
to update the default coefficients and gradients 
on the permanent file. At least two meteorological 
profiles must exist for this option to be used. 

The user should note, however, that this process is 
highly dependent on the quantity and quality of 
the meteorological data and may not always yield 
good results . 

3 . 3 LAIRS STRUCTURE AND FLOW 

The diagrams in Figure 3-1 show the hierarchy and modular 
structure of the LAIRS Program. Routine LAIRS serves as 
the executive driver of the system. After READIN has pro- 
cessed the input keyword card deck, WFILE builds the working 
file from the permanent file and, if required, calls TRNADJ 
to update the latitude gradients and diurnal/semidiurnal 
coefficients. MODADJ controls the model adjustment process; 
after it is called, the polynomial coefficients exist in 
COMMON. PARAMS is responsible for calling the proper I/O 
routines to input a trajectory point and output the calculated 
parameters and for calling MODELS, an executive that selects 
the proper atmospheric modeling routines for each trajectory 
point. 
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Figure 3-1. LAIRS Program Hierarchy (2 of 14) 



























Figure 3-1. LAIRS Program Hierarchy (3 of 14) 






U) 

I 



Figure 3-1. LAIRS Program Hierarchy (4 of 14) 








Figure 3-1 


LAIRS Program Hierarchy (5 of 14) 
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Figure 3-1. LAIRS Program Hierarchy (6 of 14) 
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Figure 3-1. LAIRS Program Hierarchy (7 of 14 
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Figure 3-1. LAIRS Program Hierarchy (9 of 14) 
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Figure 3-1. LAIRS Program Hierarchy (10 of 14) 
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Figure 3-1. LAIRS Program Hierarchy (11 of 14) 
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•DUMMY STUB FOR FUTURE IMPLEMENTATION. 


Figure 3-1. LAIRS Program Hierarchy (14 of 14) 








3 . 4 SUBROUTINE DESCRIPTIONS 


Subroutine descriptions for all the subroutines that comprise 
the LAIRS Program are given in this section. These descrip- 
tions are also included in the prologue of each respective 
subroutine. The subroutine descriptions are presented in 
alphabetical order and are preceded by a list that gives the 
name and a brief description of the function of each sub- 
routine . 
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LAIRS MODULE DESCRIPTIONS 


BARODE 

COVUP 

DCERR 

DELTIM 

DFLINT 

DIFFDE 

DIGRAD 

DJUL 

DPAR 

DSDADJ 

ENTRD 

GAS LAW 

GRAND 


Solves the barometric equation for the 
region between 90 and 125 kilometers for 
the Jacchia-Roberts model 

Deletes or adds a row or column in the 
normal matrix 

Controls differential correction error 
estimation process 

Turns packed date (YYMMDD.) into Julian 
date 

Interpolates temperature, pressure, den- 
sity, and winds from working (default) file 

Solves the diffusive equilibrium equation 
for the region between 20 and 125 kilo- 
meters for the Jacchia-Roberts model 

Produces longitudinal gradients from 
diurnal and semidiurnal coefficients 

Computes the Julian date (called by DELTIM) 

Calculates density partial derivatives for 
the differential correction process 

Adjusts the diurnal and semidiurnal coeffi- 
cients 

Reads a block of Shuttle trajectory points 
from the ENTREE file 

Recalculates temperature from pressure and 
density via gas law and performs related 
gas law computations 

Calculates integrand for the differential 
correction process 
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GWIND 
HI ALT 

HMSCON 

IBI 

lUNI 

JACCWF 

JACROB 

JR 

LAIRS 

LAP 

LATADJ 

LOWALT 

LSQPOL 

METINT 

METINl 


Calculates the geostrophic wind 

Calculates the temperature, pressure, and 
density in the region above 125 kilometers 
for the Jacchia- Roberts model 

Converts hours, minutes, and seconds into 
decimal seconds 

Performs bivariate interpolation (FORTRAN 
library routine) 

Interpolates using Lagrange's univariate 
method (FORTRAN library routine) 

Retrieves values from K file for Jacchia- 
Roberts routines 

Controls Jacchia-Roberts temperature, 
pressure, and density calculations 

Controls flow of Jacchia-Roberts model 

Controls flow of Langley Atmospheric Infor- 
mation Retrieval System (LAIRS) 

Generates altitude points for a Local 
Atmospheric Profile (LAP) 

Adjusts latitude gradient 

Computes the temperature, pressure, and 
density in the region between 90 and 125 
kilometers for the Jacchia-Roberts model 

Fits linear least-squares polynomial 
(FORTRAN library routine) 

Interpolates temperature, pressure, density, 
and winds from meteorological observation 
files 

Interpolates raw meteorological profiles 



METREF 

MODADJ 

MODELS 


NLSPOL 

OUTDO 

OUTPRF 

OUTPT 

OUTPTl 

OUTRES 

OUTSUM 

PARAMS 

PHYSCS 

PNCAL 

POLADJ 

POLCAL 

POLERR 


Produces reference meteorological profile 
(unfitted) 

Controls model adjustment process 

Retrieves temperature, pressure, density, 
and wind components for specified trajec- 
tory point via chosen model 

Controls fitting of nonlinear least squares 
polynomial 

Outputs differential correction iteration 
report 

Outputs atmospheric parameter report and 
user file 

Outputs LAIRS report header and summaries 

Outputs second half of atmospheric para- 
meter report 

Outputs residual report 

Outputs summary report of the differential 
correction coefficients 

Controls flow of the atmospheric parameter 
calculations for Shuttle trajectory or LAP 

Calculates pressure and density from temp- 
erature via physical relationships 

Calculates estimated onboard pressure meas- 
urement and Mach number 

Controls polynomial model adjustment 

Calculates temperature, pressure, and den- 
sity from polynomial coefficients 

Controls error estimation for the polynomial 
model 
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POLFIT 


POLY 

PPAR 

PREUSE 

PS GALE 
READ IN 
REST 
ROOTS 

SETREF 

SIMP 

SOLVE 

SYMINV 

TPAR 

TRNADJ 

TRNDSD 

TRNLAT 

VCMTRX 


Fits reference meteorological profile to 
polynomial model 

Controls calculation of parameters via 
polynomial model 

Calculates pressure partial derivatives 
for differential correction process 

Prepares a new meteorological file from a 
standard LAIRS output user file (USE) 

Computes pressure scale height 

Reads input keyword card deck 

Calculates local solar time 

Calculates desired complex roots of a given 
polynomial 

Sets reference values of parameters 

Performs integration by Simpson’s rule 
(FORTRAN library routine) 

Solves the normal equation and updates the 
solve-for parameters 

Inverts a symmetric matrix 

Calculates temperature partial derivatives 
for differential correction process 

Controls adjustment of translational coeffi- 
cients 

Translates atmospheric parameters in local 
solar time 

Translates atmospheric parameters in latitude 

Completes calculation of the 5x5 variance- 
covariance matrix for differential correc- 
tion error estimation 
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WCNVRT 

WFILE 

WIND 

XALFA 


Converts wind components (U and V) to 
wind speed and direction, or vice versa 

Builds working file from permanent file 

Controls the calculation of the analytic 
wind components 

Computes integrals for differential correc- 
tion variance-covariance matrix 
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SUBROUTINE BARODE ( RHO r T » P ) 


PURPOSE 

SUBROUTINE BARODE IS CALLED BY LOUALT TO COMPUTE THE ATMOSPHERIC 
DENSITY BETUEEH 90 AND 100 KM. 

METHOD! li EVALUATE THE SERIES BI 

2. COMPUTE S(R) 

3. E'v'ALUATE THE EXPONENT PI 

4. COMPUTE THE MOLECULAR MASS AT THIS ALTITUDE 

5. COMPUTE THE F -FACTORS 

6. COMPUTE THE UN CORRECTED DENSITY 

7. COMPUTE THE TEMPERATURE AND PRESSURE FROM DENSITY 


ARGUMENT 

i. 1ST: 




ARGUMENT TYPE 

I ./ 0 

DESCRIP'f ION 

Rl-IG 

R 

0 

DENSITY 

T 

R 

0 

TEMPERATURE 

P 

R 

0 

PRE 

SSURE 

C A L L- 1 N G 

SUBROUTINE 

: LOWAl. 

T 


SUBROUTINES CALLED 

: NONE 



COMMON 

BLOCK PARAMETERS 



NAME 

DIMEN 

COMMON X 

/O 

DESCRIPTION 

AC 

7 

JRCOM 

I 

MOLECULAR POWER SERIES COEFFICIENTS 

A D 

6 

JRCOM 

I 

CONSTANT COEFFICIENTS FOR p. ( N ) BE RTF 

BD 

6 

JRCOM 

I 

LINEAR COEFFICIENTS FOR B(N> SERIES 

FKL 

1 

JRCOM 

I 

FACTOR INUOLVED IN RHO COMPUTATION 

FLC4 

1 

JRCOM 

I 

MODIFYING FACTOR 

GLO 

. 1 . 

JRCOM 

I 

MEAN SURFACE GRAUITY 

HGT 

1 

JRCDM 

I 

HEIGHT OF SPACECRAFT 

RC 

1 

JRCOM 

;t; 

UNIUERSAL. GAS CONSTANT 

RCM 

]. 

JRCOM 

I • 

AMERAGE EARTH RADIUS 

RHGZ 

1 . 

^.! R (!; 0 M 

I 

ATMOSPHERIC DENSITY AT 90KM 

RLl 

1 

JRC OM 

I 

ROOT OF POLYNOMIAL IN INTEGRAND 

RL2 

J. 

JRCOM 

I 

ROOT OF POLFNOMIAL IN INTEGRAND 

TX 

1 

JRCOM 

T 

INFLECTION POINT TEMPERATURE 

TZ 

1 

JRCOM 

I 

TEMPERATURE AT HEIGHT Z 

TO 

1 

JRCDM 

I 

TEMPERATURE AT MINIMUM HEIGHT 

UC 

•n 

JRCOM 

I 

FUNCTIONAL M ALLIES AT RLl AND RL2 

UCDI 

1 

JRCOM 

I 

FACTOR INCLUDED IN RHO COMPUTATION 

WC 

'‘j 

JRCOM 

I 

FUNCTIONAL U A LUES AT RLl AND RL2 

XCDI 

1 

JRCOM 

I 

FACTOR TNCLUDEDIN RHO COMPUTATION 

X L P S 

J . 

R c; 0 M 

I 

ROOT OF polynomial IN INTEGRAND 

Y L F' 

1 

JRCOM 

Y 

R i.i 1 f I..1 F P I..I !., Y N '...1 n .1. A 1... I N I N T E G R A N D 

7 10 

j 

..) R C 0 M 

I 

M I N ]. M 1...! n H E I G H T 


F' F'! U 1.7 F-! A M M !:. 

J. P. M(]L.!;NEAUX ? COnPUTER SCIENCES CORPORATION 


Subroutine BARODE 
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SUBROUTINE COUUP ( MODE » NN r MN r IJ r ND f S » lERR ) 

purpose: 

COVUP DELETES THE ROWS AND COLUMNS OF NORMAL MATRIX 
WHICH ARE ZEROS OR ADD THE ZEROED ROWS AND COLUMNS 
TO THE INVERSE 
ARGUMENT LISTt 

ARGUMENT TYPE I/O DESCRIPTION 

MODE I I =1 DELETE ROW/COLUMN 

=2 ADD IJ ZERO ROW/COLUMNS SPECIFIED BY MN 
NN I I SI INDEX OF ROW/COLUMN INDEX TO BE REMOVED 

MN II ROW NUMBERS THAT HAVE BEEN DELETED FROM 

NORMAL MATRIX. MAXIMUM IS 20. 

IJ 10 TOTAL NUMBER OF ROW/COLUMNS DELETED 

ND II ORIGINAL DIMENSION OF NORMAL MATRIX 

S R I NORMAL MATRIX IN UPPER TRIANGULAR FORM 

CALLING subroutine: SYMINV 

SUBROUTINES CALLED: NONE 

programmer: ann welker> esc 

MODIFIED FOR LAIRS BY D. E, BOLAND MARCH 1981 


Subroutine COVUP 
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SUBROUTINE DCERR ( I SE6M » ZBOUND r DZ » T » P f D » UT j UP » UD f I ERR ) 

purpose: 

nCERR COMPUTES 1-SIGMA UNCERTAINTY FOR ATMOSPHERIC 
PARAMETERS CALCULATED THROUGH THE DC PROCESS 

method: 

1. COMPUTE PARTIAL DERIVATIVES OF TEMPERATURE» PRESSURE^ 
AND DENSITY WITH RESPECT TO THE SOLVE-PARAMETERS 

2. COMBINE THE PARTIAL DERIVATIVES WITH THE VARIANCE- 
COVARIANCE MATRIX TO FIND UNCERTAINTIES 

ARGUMENT LIST: 

ARGUMENT TYPE I/O DESCRIPTION 
ISEGM I I SEGMENT NUMBER MINUS ONE 

ZBOUND R I SEGMENT UPPER BOUNDARY ALTITUDE 

DZ R I ALTITUDE WITH RESPECT TO THE UPPER 

BOUNDARY 

T R I TEMPERATURE 

P R I PRESSURE 

D R I DENSITY 

UT R 0 COMPUTED TEMPERATURE UNCERTAINTY 

UP R 0 COMPUTED PRESSURE UNCERTAINTY 

UD R 0 COMPUTED DENSITY UNCERTAINTY 

lERR I 0 ERROR FLAG 

CALLING subroutines: POLERR 

SUBROUTINES CALLED: XALFA 

COMMON BLOCK PARAMETERS USED: 

COMMON NAME PARAMETERS I/O DESCRIPTION 


/COFCOM/ 

ct 

I 

ARRAY OF TEMPERATURE COEFFICIENTS 


CP 

I 

ARRAY OF PRESSURE COEFFICIENTS 


COVAT 

I 

VARIANCE-COVARIANCE 

MATRIX 

/SOLAND/ 

KINT 

0 

NUMBER OF INTEGRALS 
BY XALFA 

TO BE COMPUTED 


ZREF 

0 

REFERENCE ALTITUDE 

(SAME AS ZBOUND) 


A 

0 

SOLVE-FOR PARAMETER 

ARRAY 


EXTERNAL DATA SETS USED: NONE 

programmer: t. lee» esc September i981 


Subroutine DCERR 
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SUBROUTINE BELT I M ( MODE » YHD » HMS » TO JUl. » REF JUL t HT > 

PURF'OSE 

TO COMPUTE THE TIME IN SECONDS RELATIVE TO A REFERENCE DATEt 
OR THE JULIAN DATE GIUEN A PACKED CALENDAR DATE 
CALLING SEQUENCE 

CALL DELTIM( MODEr YMD . HMS f TO JUL » REF JUL ? DT ) 

CALLING SEQUENCE VARIABLES 
INPUT 

MODE = l=COMPUTE JULIAN DATE AND RETURN IN TOJUL 

2=CQMPUTE ELAPSED TIME REFERENCED TO REJULy FROM YMDj HMS 
YMD = PACKED YEAR MONTH DAY 
HMS = PACKED HOUR MINUTE SECOND 
REFJUL = REFERENCE JULIAN DATE 
OUTPUT 

TOJUL = JULIAN DATE OF YMDr HMS 

DT = ELEPSED EPHEMERIS SECONDS PAST REFJUL 

REFERENCES 

GTDS TASK NAME - TIME RELATIUE TO REFERENCE DATE (DELTIM) 
PROGRAMMER 

ALLEN L. COHEN - COMPUTER SCIENCES CORP. 


Subroutine DELTIM 
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SUBROUTINE DFLINT < N » POS » TPDW f COEF » lERR ) 


purpose; 

DFLINT CONTROLS THE DEFAULT INTERPOLATION FOR TEMPERATURE^ 
PRESSUREfDENSITY»E“W AND N-S WIND COMPONENTS* IT ALSO 
CONVERTS E~W AND N“S WINDS INTO WIND SPEED AND DIRECTION, 

method: 

1, READS IN DATA POINTS FROM WORKING FILE AND STORES VALUES IN 
ARRAYS. 

2, PASSES ARRAYS TO SUBROUTINE lUNI FOR INTERPOLATIONS, 

3, USES LOGARITHMIC INTERPOLATION FOR PRESSURE AND DENSITY, 

4, TRANSLATES TEMPERATURE? PRESSURE? AND DENSITY FOR LATITUDINAL? 
DIURNAL? AND SEMIDIURNAL VARIATIONS IF DESIRED, 

5, CONVERTS E-W AND N-S WIND COMPONENTS INTO WIND SPEED AND 
DIRECTION . 


ARGUMENT LISTt 

ARGUMENT TYPE 
N I 


POS R 


TPDW R 


COEF R 


lERR I 


I/O 

I 


I 


0 


0 


0 


DESCRIPTION 
ATMOSPHERIC SEGMENT 

=0 ? NO INTERPOLATION NEEDED 
(FLAG FOR MET, DATA) 

=1 ? UPPER 
=2 ? LOWER 
3-WORD array: 

POS(l) = ALTITUDE (KM) 

P0S(2) = LATITUDE (DEGREES) 

P0S(3) = LOCAL SOLAR TIME (HOURS) 
P0S(4) = LONGITUDE 

7-WORD array: 

TEMPERATURE 
PRESSURE 
DENSITY 
E-W WIND(U) 

N-S WIND(V) 


TPDW( 1 ) 
TPDW(2> 
TPDW(3) 
TPDW(4) 
TPDU(5) 
TPDW(6) 
TPDW(7) 


(DEG, K) 

(N/M2) 

(KG/M3) 

(M/SEC) 

(M/SEC) 

WIND SPEED(WS) (M/SEC) 

WIND DIR,(WD) (DEG. FROM N) 
15-WORD array: 

LATITUDE GRADIENTS OF TEMPERATURE? 
PRESSURE? AND DENSITY? DIURNAL 
AMPLITUDE? PHASE? SEMIDIURNAL AMPLI- 
TUDE? PHASE? FOR T? P? D 
ERROR FLAG 

=0 ? NO ERROR 


CALLING subroutines: PARAMS 
SUBROUTINES CALLED: WCNVRT?IUNI 


Subroutine DFLINT (1 of 2) 
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FUNCTIONS CALLED 


TRNLATrTRNDSD 


COHNON BLOCK PARAMETERS USED: 


COMMON NAME 

PARAMETER 

I/O 

DESCRIPTION 

/USECOM/ 

ILATL 

I 

INCLUDE LAT» VARIATIONS IN L.A. 
0=N0 
1 = YES 


IDSDL 

I 

INCLUDE D/SD VARIATIONS IN L»A» 
0=N0 
1=YES 


ILATU 

I 

INCLUDE LAT* VARIATIONS IN U.A» 
0=N0 
1 = YES 


IDSDU 

I 

INCLUDE D/SD VARIATIONS IN U*A. 
0=>=N0 
1=YES 


INTP 

I 

INTERPOLATER TYPE 

1«1ST ORDER LAGRANGE 
2=2ND ORDER LAGRANGE 


IDB6 

I 

DEBUG FLAG 

/FILCOM/ 

NPRINT 

I 

PRINTED OUTPUT LUN 


NWRK 

I 

INPUT LUN 

/INTCOM/ 

IPT 

I/O 

POINTER USED IN lUNI 
INITIALLY SET TO -1 


IPTFLG 

I 

FLAG FOR FIRST CALL TO DFLINT 


ALT 

I 

ARRAY FOR ALTITUDES PASSED TO lUNI 


TPDWAR 

I 

ARRAY FOR TfP»Dr AND WIND COMP* 
PASSED TO lUNI 

EXTERNAL DATA 

ARR 

SETS used: 

I 

ARRAY FOR COEFF. PASSED TO 
TRANSLATING FUNCTIONS 

NAME 

LUN 

I/O 

OPERATIONS PERFORMED 

WRK 

NWRK 

READ 

OUT 

NPRINT WRITE 

programmer: r« 

A* KUSESKI* 

CSC 



Subroutine DFLINT (2 of 2) 
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SUBROUTINE DIFFDE ( RHO » T » P ) 

PURPOSE 

DIFFDE IS CALLED BY LOWALT TO COMPUTE THE ATMOSPHERIC DENSITY 
VALUES FROM 100 KM TO 125 KM» 


method: 

1. COMPUTE THE EXPONENTS QI 
2» COMPUTE THE F FACTORS 
3» COMPUTE THE CONSTITUENT DENSITIES 
A* COMPUTE THE TOTAL DENSIITY 

5» COMPUTE TEMPERATURE AND PRESSURE FROM DENSITY 


ARGUMENT LIST! 


ARGUMENT TYPE 
RHO R 

T R 

P R 


CALLING subroutine: 
SUBROUTINES CALLED: 


I/O DESCRIPTION 

0 DENSITY 

0 TEMPERATURE 

0 PRESSURE 

LOUALT 

NONE 


COMMON 

NAME 

BLOCK 

DIMEN 

PARAMETERS 

COMMON 

I/O 

DESCRIPTION 


ADT 

6 

JRCOM 

1 

THERMAL DIFFUSION COEFFICIENTS 


CM 

6 

JRCOM 

I 

MASSES OF ATMOSPHERIC CONSTITUENTS 

IN 

FD125 

5 

JRCOM 

0 

GM/MOLE 

INDIVIDUAL NUMBER DENSITIES 


FKL 

1 

JRCOM 

I 

FACTOR INVOLVED IN RHO COMPUTATION 


FLC4 

1 

JRCOM 

I 

MODIFYING FACTOR 


6L0 

1 

JRCOM 

I 

MEAN SURFACE GRAVITY 


HGT 

1 

JRCOM 

I 

HEIGHT OF SPACECRAFT 


RC 

1 

JRCOM 

I 

UNIVERSAL GAS CONSTANT 


RCM 

1 

JRCOM 

I 

AVERAGE EARTH RADIUS 


RLl 

1 

JRCOM 

I 

ROOT OF POLYNOMIAL IN INTEGRAND 


RL2 

1 

JRCOM 

I 

ROOT OF POLFNOMIAL IN INTEGRAND 


SD 

5 

JRCOM 

I 

COEFFICIENTS FOR NUMBER DENSITY AT 

lOOKM 

TCIL 

1 

JRCOM 

I 

TEMPERATURE AT 100 KM* 


TINF 

1 

JRCOM 

I 

EXOSPHERIC TEMPERATURE 


TZ 

1 

JRCOM 

I 

TEMPERATURE AT HEIGHT Z 


UC 

2 

JRCOM 

I 

FUNCTIONAL VALUES AT RLl AND RL2 


VCDI 

1 

JRCOM 

I 

FACTOR INCLUDEDIN RHO COMPUTATION 


UC 

2 

JRCOM 

I 

FUNCTIONAL VALUES AT RLl AND RL2 


XCDl 

1 

JRCOM 

I 

FACTOR INCLUDEDIN RHO COMPUTATION 


XLPS 

1 

JRCOM 

1 

ROOT OF POLYNOMIAL IN INTEGRAND 


YLPS 

1 

JRCOM 

I 

ROOT OF POLYNOMIAL IN INTEGRAND 


ZD 

7 

JRCOM 

I 

POUER SERIES COEFFICIENTS FOR RHO(l0O>/M 


PROGRAMMER 

P. MOLINEAUXr COMPUTER SCIENCES CORPORATION 


Subroutine DIFFDE 
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SUBROUTINE D I GRAD ( DA » DPH » SDA r SDPH t ST r DLONG ) 

PURPOSE} 

DIGRAD PRODUCES THE LONGITUDINAL GRADIENT CORRESPONDING TO 

THE INPUT DIURNAL AND SEMIDIURNAL COEFFICIENTS (FOR ANY QUANTITY) 

METHOD} 

1. EVALUATE THE TIME DERIVATIVE OF THE DIURNAL VARIATION, 

2, EVALUATE THE TIME DERIVATIVE OF THE SEMIDIURNAL VARIATION, 

3* CONVERT THE TOTAL DER I V AT I VE . TO UNITS OF DEGREES LONGITUDE, 

4, RETURN, 

ARGUME:NT LIST 

ARGUMENT TYPE I/O 

D A R I 

DPH R I 

SDA R I 

SDPH R I 

ST R I 

D I... 0 N G R 0 


C A L 1. 1 N G S U B R 0 U T I M E } U I M D 

SilBROUTINES CALLED} NONE 

C 0 M M 0 N B L.. 0 C K P i^i R A M E T E- R S U S P. D } 
COMMON NAME PARAMETER I/O 

/CONCOM/ TUG PI I 

HTD I 

EXTERNAL DATA SETS USED! NONE 
PROGRAMMER} D, E, BOLAND, CSC 


Subroutine DIGRAD 


DESCRIPTION 
TWO TIMES PI 

HOURS TO DEGREES CONVERSION FACTOR 


DESCRI PTION 

A M !"■ L. I T U D E 0 F D I U R N A L.. V A R I A T I 0 N 
PHASE OF DIURNAL VARIATION 
A M F' 1.. I T LI D F: 0 F S E M I D I U R N A 1. V A R I A T I 0 N 
PHASE OF SEMIDIURNAL VARIATION 
LOCAL SOLAR TIME AT WHICH DERIVATIVE, 
IS TO BE CALCULATED (HOURS) 

L. 0 N G I T U D I N A L. D E R I V A T I V E C 0 R R E S F^' 0 N D I N G 
TO INPUT D&SD COEEFICIENTS 


3-33 



REAL FUNCTION r.UUL ( Y ^ XM , H » HR » TM r SEC ) 


FORTRAN FUNCTION SUBPROGRAM 
PURPOSE 

TO COMPUTE THE MODIFIED JULIAN DATE OF A GIVEN GREGORIAN 
DATE AFTER 1950*0 (JULIAN DATE - 2430000 > 

USAGE 

REAL FUNCTION D JUL ( Y i XM » D f HR » TM r SEC ) 

INPUT 

Y - YEAR 
Xft - MONTH 
D - DAY 
HR - HOURS 
TM - MINUTES 
SEC " SECONDS 

OUTPUT 

DJUL - THE MODIFIED JULIAN DATE 


REMARKS 

INPUT DATE MUST OCCUR AFTER 1950*0 

SUBROUTINES REQUIRED 
NONE 

PROGRAMMER 

M* MC GARRY 


Real Function DJUL 



SUBROUTINE DPAR ( T » DTDS f RI NT i DDDS > 

PURPOSE? 

DPAR COMPUTES PARTIAL DERIVATIVES OF DENSITY WITH RESPECT TO 
SOLVE-FOR PARAMETERS 

METHOD? 

COMPUTE THE PARTIAL DERIVATIVES OF DENSITY WITH RESPECT TO 
SOLVE-FOR PARAMETERS USING TEMPERATURE r TEMPERATURE PARTIAL 
DERIVATIVES AND INTEGRALS COMPUTED BY SIMP 


ARGUMENT 

ARGUMENT 

LIST? 

TYPE 

I/O 

DESCRIPTION 

T 

R 

I 

COMPUTED TEMPERATURE 

DTDS 

R 

I 

COMPUTED TEMPERATURE PARTIAL DERIVATIVES 

RINT 

R 

I 

COMPUTED INTEGRALS 

DDDS 

R 

0 

LOG DENSITY PARTIAL DERIVATIVES COMPUTED 


CALLING SUBROUTINE? NLSPOL 
SUBROUTINES CALLED? NONE 


COMMON BLOCK PARAMETERS USED? 
COMMON NAME PARAMETER I/O 

/SOLANO/ IMAX I 

KS0LA4 I 

NPOL I 

A 1 


DESCRIPTION 

MAX* NUMBER OF INTEGRALS 

=O.DOES NOT SOLVE FOR BOUNDARY TEMP* 

=1 SOLVES FOR BOUNDARY TEMP. 

ORDER OF TEMP, POLYNOMIAL 
SOLVE-FOR parameter ARRAY 


EXTERNAL DATA SETS USED? NONE 


PROGRAMMER? T, LEEf CSC MARCH 1981 


Subroutine DPAR 
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SUBROUTINE DSD/)Dj<NFf Ff Al fPl »A2*P2»0FLATf TAUyRLAT) 

purpose: 

DSDADJ COMPUTES NEW DIURNAL AND SEMIDIURNAL COEFFICIENTS 
FOR ANY INPUT PARAMETER F <T# Pt OR D> FROM METEOROLOGICAL 
DATA. 

method: 

1. TRANSLATE EACH INPUT PARAMETER F (UP TO THREE POINTS — 

ONE FROM EACH MET FILE) FROM ITS LATITUDE RLAT TO THE 
MODEL REFERENCE LATITUDE REFLAT TO REMOVE THE LATITUDE 
COMPONENT OF THE VARIATION OF F. 

2. LOAD THE *C* ARRAY t CONSISTING OF THE COSINE TERMS 
TO PE USED IN SUPSEOUENT CALCULATIONS. 

3. IF ONLY THE DIURNAL COMPONENT OR ONLY THE SEMIDIURNAL 
COMPONENT HAS KNOWN DEFAULT VALUES* DO A LINEAR LEAST- 
SQUARES FIT FOR THE AMPLITUDE OF THE KNOWN COMPONENT. 

4 . IF DOTH COMPONENTS HAVE KNOWN DEFAULT VALUES* THEN 
SOLVE FOR THE AMPLITUDES OF THE DIUPNAL AND SEMIDIURNAL 
VARIATIONS (DEFAULT VALUES OF THE PHASES ARE USED). 

IF THIRE ARE ONLY TWO POINTS <NF»2)* CALCULATE THE 
AVERAGE VALUE OF F (FAV) USING THE DEFAULT VALUES OF 
THE COEFFICIENTS AND USE IT TO CALCULATE EXACTLY THE 
THE values of THE DIURNAL AND SEMIDIURNAL AMPLITUDES, 

IF NF»3, DOTH FAV AND THE AMPLITUDES MAY BE CALCULATED 
exactly* THOUGH ONLY THE AMPLITUDES ARE CALCULATED. 

3, RETURN. 

ARGUMENT list: 

ARGUMLnT ( YPE I/O DESCRIPTION 

NF R 1 NUMBER OF POINTS OF F 

F R I parameter to be TRANSLATED 

Al R I/O DIURNAL amplitude OF F 

PI R I/O DIURNAL phase OF F 

A2 R I/O SEMIDIURNAL AMPLITUDE OF F 

P2 R I/O SEMIDIURNAL PHASE OF F 

TAU R I LOCAL SOI AR TIME OF MET PROFILE 

DFLAT R 0 LATITUDE GRADIENT OF F 

RLAT R I LATITUDE OF MET PROFILE 

CALLING SUBROUTINE! TRNADJ 
SUBROUTINES called: TRNLAT 

COMMON BLOCK PARAMETERS USED: 

COMMON NAME PARAMETER I/O DESCRIPTION 
/USECOM/ REFLAT I MODEL REFERENCE LATITUDE 

PEFTAU I MODEL REFERENCE LOCAL SOLAR TIME 

EXTERNAL DATA SETS USED! NONE 

programmer: D. E, BOLAND* CSC* march 1981 


Subroutine DSDADJ 
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SUBROUTINE ENTRD < NZ f I END » I ERR ) 


purpose: 

ENTRD READS SHUTTLE TRAJECTORY POINTS FROH THE ENTREE FILE AND 
STORES THEN IN COMMON /CALCOM/ » 


method: 

1. READ A BLOCK OF TRAJECTORY POINTS FROM THE ENTREE FILE AND 
STORE THEM IN /CALCOM/ . 

2» CALL RLST TO CONVERT LONGITUDE AND GMT INTO LOCAL SOLAR TIME 
(ALSO STORED IN COMMON)* 

3* RETURN* 


ARGUMENT LIST: 


ARGUMENT TYPE 
lERR I 

NZ I 

lEND I 


I/O DESCRIPTION 

0 ERROR flag: 

=0» NO ERROR 

1 NUMBER OF POINTS IN BLOCK 

0 END-OF-FILE FLAG (NO, PTS* IN LAST BLOCK) 


CALLING subroutine: LAIRS 


SUBROUTINES CALLED** RLST 


COMMON BLOCK PARAMETERS USED** 
COMMON NAME PARAMETER I/O 
/CALCOM/ TRJ(5»10) 0 


UTRJ(5» 10) 
V( 10) 
ALPHA( 10) 


/FILCOM/ NPRINT I 

NENT I 

/USECOM/ IDBG I 

NSAMP I 


DESCRIPTION 

BLOCK OF 10 TRAJECTORY POINTS** 

ALT I TUDE » L AT I TUBE » LONGITUDE? GMT ? 

LOCAL SOLAR TIME 

TRAJECTORY POINT UNCERTAINTIES 

BLOCK OF SHUTTLE VELOCITIES 

BLOCK OF SHUTTLE ANGLES-OF- ATTACK 

PRINTED OUTPUT LUN 

ENTREE FILE LUN 

DEBUG FLAG 

SAMPLING INTERVAL 


EXTERNAL DATA 
NAME 
OUT 
ENT 


SETS used: 
LUN 

NPRINT 

NENT 


I/O OPERATIONS PERFORMED 
WRITE(DEBUG ONLY) 

READ 


programmer: D, E, BOLAND? CSC» march 1981 


Subroutine ENTRD 
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SUBROUTINE GASLAW ( Z t TPDW f WTM y SCH » PCDEV ) 

purpose: 

GASLAW CALCULATES CERTAIN PARAMETERS ASSOCIATED WITH THE IDEAL 
GAS LAW AND RECOMPUTES THE TEMPERATURE TO AGREE WITH THIS LAW* 

method: 

1* CONVERT MOLECULAR WEIGHT TO MKS AND CALCULATE THE RATIO 
OF PRESSURE TO DENSITY. 

2. COMPUTE THE PERCENT DEVIATION FROM THE GAS LAW BY 
COMPUTING THE GAS CONSTANT FROM THE INCOMING PARAMETERS 
AND COMPARING IT TO THE KNOWN VALUE. 

3. RECOMPUTE TEMPERATURE IN ACCORDANCE WITH THE GAS LAW. 

4. COMPUTE THE MOLECULAR WEIGHT. 

5. COMPUTE THE PRESSURE SCALE HEIGHT. 

6. RETURN. 


ARGUMENT LIST 
ARGUMENT 

• 

♦ 

TYPE 

I/O 

DESCRIPTION 


Z 

R 

I 

SHUTTLE ALTITUDE <KM) 


TPDW(7) 

R 

I 

CALCULATED ATMOSPHERIC 

parameters: 

WTM 

R 

0 

TrPyDyUyVfWSf WD (ONLY 
MOLECULAR WEIGHT 

TypyD USED) 

SCHT 

R 

0 

PRESSURE SCALE HEIGHT 


PCDEV 

R 

0 

PERCENT DEVIATION FROM 

GAS LAW 


CALLING subroutine: PARAMS 

SUBROUTINES CALLED: NONE 

COMMON BLOCK PARAMETERS USED: 

COMMON NAME PARAMETER I/O DESCRIPTION 


/USECOM/ 

BOUNDl 

I 

BOUNDARY BETWEEN UPPER AND LOWER 
ATMOSPHERIC MODELS <DEFAULT=90 KM) 


IGAS 

I 

FLAG FOR DERIVATION OF T FROM GAS LAW 

/JRCOM/ 

RAV 

I 

AVERAGE EARTH RADIUS (KM) 


RC 

I 

GAS CONSTANT (MKS UNITS) 


GLO 

I 

SEA-LEVEL ACCELERATION OF 
GRAVITY (M/SEC**2) 


CM 

I 

SEA-LEVEL MEAN MOLECULAR WEIGHT 


EXTERNAL DATA SETS USED: NONE 

programmer: D. E. BOLANDr CSC JULYf 1981 


Subroutine GASLAW 
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SUBROUTINE GRAND(X»WK) 

purpose: 

GRAND PREPARES THE INTEGRAND FOR SIMPSON INTEGRATOR 


method: 

USING TEMPERATUREf MOLECULAR WEIGHTf GAS CONSTANT» AND GRAVITATIONAL 
ACCELERATION^ COMPUTE THE INTEGRAND OF INTEGRALS NEEDED FOR PRESSURE 
AND DENSITY COMPUTATION AT A GIVEN ALTITUDE 


ARGUMENT LIST: 
ARGUMENT TYPE I/O 
X I 

WK 0 


DESCRIPTION 

ALTITUDE 

COMPUTED INTEGRAND 


CALLING subroutine: SIMP 
SUBROUTINES CALLED: NONE 


COMMON BLOCK PARAMETERS USED: 


COMMON NAME 

PARAMETERS 

I/O 

DESCRIPTION 

/JRCOM/ 

RE 

I 

MEAN RADIUS OF THE EARTH 


rstar 

I 

GAS CONSTANT 


GS 

I 

SEA-LEVEL GRAV, ACCELERATION 


UM 

I 

MEAN MOLECULAR WEIGHT 

/SOLAND/ 

NPOL 

I 

ORE^ER OF POLYNOMIAL 


IMAX 

I 

MAX. NUMBER OF INTEGRANDS 


ZREF 

I 

REFERENCE ALTITUDE: SEGMENT BOUNDARY 


A 

I 

SOLVE-FOR PARAMETERS 

EXTERNAL DATA 

SETS used: 

NONE 


programmer: t» 

LEEf CSC 


MARCH 1981 


Subroutine GRAND 
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SUBROUTINE GUINO< THETA FOfOPLATFDPLONGfUfV) 

purpose; 

GUIND CALCULATES THE GEOSTROPHIC UIND FROH A GIVEN PRESSURE 
DISTRIBUTION. 

method; 

1. CONVERT THE LATITUDINAL AND LONGITUDINAL PRESSURE GRADIENTS 
FROM UNITS OF 1/DEBREES TO UNITS OF 1/METERS. 

2. CALCULATE THE WINDS FROM THE GEOSTROPHIC WIND EQUATON. 

3* RETURN* 

ARGUMENT LIST; 

ARGUMENT TYPE I/O DESCRIPTION 

THETA R I LATITUDE 

DPLAT R I LATITUDE PRESSURE GRADIENT 

DPLONG R I LONGITUDE PRESSURE GRADIENT 

D R I DENSITY 

U R 0 EAST-WEST WIND COMPONENT (M/SEC) 

V R 0 NORTH-SOUTH WIND COMPONENT <M/SEC> 

CALLING SUBROUTINE; UIND 

SUBROUTINES CALLED; NONE 

COMMON BLOCK PARAMETERS USED; 

COMMON NAME PARAMETER I/O DESCRIPTION 

/CONCOM/ R I EARTH EQUATORIAL RADIUS 

ANG I EARTH ANGULAR VELOCITY 

RTD I RADIANS TO DEGREES CONVERSION 

FACTOR 

EXTERNAL DATA SETS USED*. NONE 
programmer; D. E. BOLANDr CSC 


Subroutine GWIND 
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SUBROUTINE H I ALT < RHO t T » P ) 


PURPOSE 

HIALT IS CALLED BY JACROB TO FURNISH DENSITY VALUES ABOVE 125 KM. 

method: 1. EVALUATE CONSTANTS THAT RECUR FREQUENTLY 

2. COMPUTE NUMBER DENSITIES OF N2» AR»HE»02f 0 

3. COMPUTE THE SEASONAL LATITUDINAL DENSITY OF HELIUM AND 
MULTIPLY IT BY HELIUM'S CONTRIBUTION 

4. CONVERT THE NUMBER DENSITIES TO MASS DENSITIES 

5. COMPUTE HYDROGEN DENSITY 

6. SUM THE SONSTITUENT DENSITIES TO GIVE THE TOTAL» 
UNCORRECTED DENSITY 

7. COMPUTE TEMPERATURE AND PRESSURE FROM DENSITY 


ARGUMENT LIST 
ARGUMENT 

• 

TYPE 

I/O 

DESCRIPTION 

RHO 

R 

0 

ATMOSPHERIC DENSITY 

T 

T 

0 

TEMPERATURE 

P 

R 

0 

PRESSURE 


CALLING subroutine: JACROB 


CALLING SEQUENCE 
CALL HIALT (RHO) 

WHERE RHO IS THE OUTPUT ATMOSPHERIC DENSITY (PARTIAL» DUE TO 
THE MODELLING OF EFFECTS ABOVE 125 KM.) 

COMMON BLOCK VARIABLES 


NAME 

DIMEN 

COMMON 

I/O 

DESCRIPTION 

ADT 

6 

JRCOM 

I 

THERMAL DIFFUSION COEFFICIENTS 

AVG 

1 

JRCOM 

I 

AVOGADRO'S NUMBER 

CF 

1 

JRCOM 

I 

CORRECTION FACTOR 

CFL 

5 

JRCOM 

I 

POWER SERIES COEFFICIENTS FOR CF 

CM 

6 

JRCOM 

I 

MASSES OF ATMOSPHERIC CONSTITUENTS 

D125 

5 

JRCOM 

I 

LOG INDIVIDUAL NO* DENSITIES AT 12 

DELO 

1 

JRCOM 

I 

SUN'S LATITUDE 

GLO 

1 

JRCOM 

I 

MEAN SURFACE GRAVITY 

HGT 

1 

JRCOM 

I 

SPACECRAFT HEIGHT 

PHJ 

1 

JRCOM 

I 

S/C LATITUDE 

RC 

1 

JRCOM 

I 

UNIVERSAL GAS CONSTANT 

RCM 

1 

JRCOM 

I 

AVERAGE EARTH RADIUS 

TINF 

1 

JRCOM 

I 

EXOSPHERIC TEMPERATURE 

TUOPI 

1 

CONCOM 

I 

TWO TIMES PI 

TX 

1 

JRCOM 

I 

INFLECTION POINT TEMPERATURE 

TO 

1 

JRCOM 

I 

TEMPERATURE AT MINIMUM HEIGHT 

ZJX 

1 

JRCOM 

I 

INFLECTION POINT HEIGHT 

ZJO 

1 

JRCOM 

I 

MINIMUM HEIGHT 


PROGRAMMER 

A.K. KAPOORr COMPUTER SCIENCES CORPORATION 


Subroutine HIALT 


3-41 



FUNCTION HMSCON(HMS) 

purpose: 

HMSCON CONOERTS HOURS» MINUTES, SECONDS TO SECONDS. 

method: 

1. DECOMPOSE HOURS, MINUTES, SECONDS AND COMPUTE SECONDS. 

2. RETURN. 

ARGUMENT LIST*. 

ARGUMENT TYPE I/O DESCRIPTION 

HMS R I HOUR, MINUTES, SECONDS (HHMMSS.SSS) 

HMSCON R 0 SECONDS CORRESPONDING TO HMS 

CALLING subroutines: READIN 

SUBROUTINES CALLED: NONE 

COMMON BLOCK PARAMETERS USED: NONE 

EXTERNAL DATA SETS USED: NONE 

programmer: D. E. BOLAND, CSC, march, 1981 


Function HMSCON 


3-42 



SUBROUTINE JACCUF(T) 


UERBION 2.3 
PURPOSE 

JACCUF IS CALLED TO FILL COMMON BLOCK TKPTC WITH A TIME ITr 
MAGNETIC ACTIVITY INDICES AND EXOSPHERIC TEMPERATURES. 


METHOD! 1. 

2 « 
3. 


4 , 


6 . 


7. 


8 . 
9 , 

10 . 
11 . 


CONVERT THE STARTING AND ENDING TIMES TO MODIFIED JULIAN 

JULIAN DATES» USING SUBROUTINE DELTIM 

READ THE JACCHIA PERMANENT FILE HEADER FILE KPILE 

IF THE ARC WILL FIT IN COMMON/TKPTC/ ( END-START<20 DAYS>r 

GO TO STEP 5 

SINCE THE ARC IS MORE THAN 20 DAYS» SIMPLY DETERMINE THE 
RECORD NUMBER CONTAINING THE REQUEST TIME» READ IT INTO 
COMMON/TKPTC/» AND PROCEED TO STEP 10 IF IT IS ON THE 
PERMANENT FILEt OR WRITE AN ERROR MESSAGE AND RETURN IF 
IT IS NOT. 

THE ARC IS 20 DAYS OR LESS. READ THE RECORD CONTAINING 
THE START OF THE ARC INTO THE COMMON ARRAYS, IF IT IS 
PRESENT IN THE PERMANENT FILE, OTHERWISE, WRITE AN ERROR 
MESSAGE AND RETURN. 

THE<21-DAY-ARC TIME SPAN MAY RESIDE ON TWO RECORDS. CHECK 
IF THE WHOLE ARC IS ON THE RECORD. IF SO, PROCEED TO 
STEP 10. 

IF THE RECORD CONTAINING THE REST OF THE ARC IS NOT ON 
THE FILE, WRITE AN ERROR MESSAGE AND RETURN. IF IT IS 
PREPARE TO READ IT BY SLIDING THE DATA PERTAINING TO 
THE FIRST PART OF THE ARC INTO THE STARTING LOCATIONS 
OF ITS ARRAY IN COMMON/TKPTC/. 

READ THE REST OF THE ARC INTO THE LAST LOCATIONS OF THE 
ARRAYS. 

RESET THE STARTING DAY OF THE COMMON BLOCK TO ACCOUNT 
FOR THE SHIFTING. 

EXPAND THE PACKED KP VALUES SO THAT EACH LOCATION IN 
ARRAY KP WILL HAVE ONLY ONE KP NUMBER, 

RETURN TO THE CALLING PROGRAM. 


CALLING SEQUENCE 

CALL JACCWF(T) 

WHERE T IS A TIME AT WHICH A KP AND TC ARE REQUIRED, IN MJD. 
CALLING SUBROUTINE! JACRDB 
SUBROUTINES CALLED! OPENMS , READMS , CLOSMS 


COMMON BLOCK VARIABLES! 


NAME 

IT 

KP 

TC 


DIMEN COMMON 
1 TKPTC 
21 ,8TKPTC 
20 TKPTC 


I/O DESCRIPTION 

TIME OF FIRST DAY OF TC DATA IN COMMON 
MAGNETIC ACTIVITY 3-HGUR INDICES 
EXOSPHERIC TEMPERATURE 


EXTERNAL DATA SETS USED! 
NAME LUN I/O 

KPILE 7S READ 


OPERATION PERFORMED 

READS THE KP VALUES FROM THE FILE AND 
STORES THEM IN COMMON TKPTC 


PROGRAMMER 
A.K. KAPOOR 


COMPUTER SCIENCES CORPORATION 


Subroutine JACCV7F 
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SUBROUTINE JACROB ( RHO r T » P ) 


PURPOSE 

JACROB IS CALLED BY JR TO FURNISH ATMOSPHERIC DENSITY USING 
JACCHIA'S 1971 ATMOSPHERIC MODEL AS MODIFIED BY ROBERTS IN 
1971 » COMPUTATIONS PERFORMED ARE TIME DEFINITIONf EXOSPHERIC 
TEMPERATURE^ CORRECTIONS TO THE ATMOSPHERIC DENSITY^ AND 
DETERMINATION OF WHICH HEIGHT RANGE TO USE IN COMPUTING THE 
DENSITY » 

method; 

1» IF THE SPACECRAFT HEIGHT IS LESS THAN 90KM» SET DENSITY 
TO ITS VALUE AT 90KM AND RETURN 
2» CHECK IF COMMON/TKPTC/ CONTAINS DATA FOR THE TIME IN 

QUESTION, IF NOT» CALL SUBROUTINE JACCWF TO FILL IT WITH 
21 DAYS OF KP AT 3-HOUR INTERVALS < KP < I ) j J » 1=1 » 21 r J=1 ^ 8 ) 
AND 20 DAYS OF TC < TC< I ) » 1=1 » . , , > 20 ) , 

3, COMPUTE THE LOCAL TEMPERATURE 

4, COMPUTE CORRECTIONS TO THE LOCAL TEMPERATURE DELTA TG) AND 
TO DENSITY OCCASIONED BY GEOMAGNETIC FLUX. 

5, CORRECT THE LOCAL TEMPERATURE FOR GEOMAGNETIC EFFECTSr 
GIVING THE EXOSPHERIC TEMPERATURE 

6, COMPUTE THE INFLECTION POINT TEMPERATURE (TEMPERATURE 
AT A HEIGHT OF 125KM ABOVE THE SUBSATELLITE POINT), 

7, COMPUTE AN ADDITIONAL CORRECTION TO THE DENSITY DUE 

TO SEMI-ANNUAL VARIATIONS AND ADD IT TO THE CORRECTION 
TO DENSITY ALREADY COMPUTED 

8, COMPUTE THE CORRECTION TO THE DENSITY DUE TO SEASONAL 
LATITUDINAL VARIATIONS IN THE LOWER THERMOSPHERE » ADDING 
IT TO THE CORRECTION TO DENSITY ALREADY COMPUTED, 

9, CALL THE SUBROUTINES TO COMPUTE UNCORRECTED DENSITY BY 
ROBERTS formulas; 

IF Z<125KM CALL LOWALT 
IF Z>125KM CALL HIALT 

10, CORRECT THE DENSITY 

11, CONVERT DENSITY TO MKS UNITS 


ARGUMENT LIST; 




ARGUMENT 

TYPE 

I/O 

DESCRIPTION 

RHO 

R 

0 

THE OUTPUT ATMOSPHERIC DENSITY 

T 

R 

0 

TEMPERATURE 

P 

R 

0 

PRESSURE 

CALLING subroutine; 

JR 


SUBROUTINES CALLED; 

JACCWFf 

HIALTf LOWALT 


Subroutine JACROB (1 of 2) 
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REFERENCES 

JACCHIA^ L. REVISED STATIC MODELS OF THE THERMOSPHERE AND 

EXOSPHERE WITH EMPIRICAL TEMPERATURE PROFILES^ S«A»0*» SPECIAL 
REPORT NO. 332f MAY 5»1971. 

ROBERTSf C. E.» JR.» AN ANALYTICAL MODEL FOR UPPER ATMOSPHERE 

DENSITIES BASED UPON JACCHIA'S 1970 MODELS. CELESTIAL MECHANICS 
4. (1971 ) . 

COMMON BLOCK VARIABLES 


NAME 

DIMEN 

COMMON 

I/O 

DESCRIPTION 

DELO 

1 

JRCOM 

0 

SUN'S LATITUDE 

DTR 

1 

CONCOM 

I 

DEGREES TO RADIANS CONVERSION FACTOR 

EBASE 

1 

CONCOM 

I 

BASE OF NATURAL LOGARITHM 

HCJ 

1 

JRCOM 

0 

HOUR ANGLE OF THE SUN (DEGREES) 

HOT 

1 

JRCOM 

I 

SPACECRAFT HEIGHT 

IT 

1 

TKPTC 

I 

STARTING MJD OF TC VALUES IN THIS COMMON 

KP 

21.8 

TKPTC 

I 

21 DAYS OF 3-HOUR GEOMAGNETIC INDICES 
STARTING WITH DAY IT - 1. 

PHJ 

1 

JRCOM 

0 

GEODETIC LATITUDE (RADIANS) 

PI 

1 

CONCOM 

I 

PI 

RH02 

1 

JRCOM 

I 

ATMOSPHERIC DENSITY AT HEIGHT = ZJO 

RTD 

1 

CONCOM 

I 

RADIANS TO DEGREES CONVERSION FACTOR 

TC 

20 

TKPTC 

I 

20 DAYS OF NIGHTTIME MINIMUM GLOBAL 
EXOSPHERIC TEMPERATURES START ON DAY IT 

TINF 

1 

JRCOM 

0 

EXOSPHERIC TEMPERATURE 

TWOPI 

1 

CONCOM 

I 

TWO TIMES PI 

TX 

1 

JRCOM 

0 

INFLECTION POINT TEMPEATURE 

ZJO 1 

SUBROUTINES USED 

JRCOM 

I 

MINIMUM HEIGHT 

HIALT 

- COMPUTES RHO 

AT 

ALTITUDES ABOVE 125 KM. 

JACCWF 

- FILLS 

COMMON 

TKPTC WITH DATA FROM PERMANENT FILE. 

LOWALT 

- COMPUTES RHO 

AT 

ALTITUDES BETWEEN 90 KM AND 125 KM 


PROGRAMMER 

J.P. MOLINEAUX. COMPUTER SCIENCES CORPORATION. 


Subroutine JACROB (2 of 2) 
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SUBROUTINE JR<N»POSf TPDUr lERR) 


PURPOSE : JR ROUTINE CONTROLS TEMPERATURE* PRESSURE* AND 

DENSITY CALCULATIONS FOR THE UPPER (JACCHIA- 
ROBERTS) SEGMENT OF THE ATMOSPHERE. 


METHOD 


1* COMPUTE THE JULIAN DATE BY CALLING DELTIM ROUTINE 
AND STORE THE UALUE IN COMMON. 

2. COMPUTE THE DECLINATION OF THE SUN. STORE IN 
COMMON. 

3. CALL JACROB SUBROUTINE TO COMPUTE THE TEMPERATURE* 
PRESSURE*AND DENSITY (T*P*D) VALUES. 

4. STORE VALUES IN ARRAY TPDU. 

5. RETURN. 


ARGUMENT LIST: 

ARGUMENT TYPE I/O 


DESCRIPTION 


N 

I 

I 

MODEL TYPE 
=1* DEFAULT 
=2* ADJUSTED 

POS 

R 

I 

FOUR WORD ! 

POS(l) 

P0S(2) 

P0S(3) 

P0S(4) 

ARRAY 

ALTITUDE(KM) 
LATITUDE(DEG) 
LOCAL SOLAR 
TIME (HOURS) 
LONGITUDE 

TPDW 

R 

0 

SEVEN WORD 
TPDH(l) -- 
(2) — 

(3) — 

(4) — 

(5) — 
(A) — 
(7) — 

ARRAY 

TEMP(DEG KELVIN) 
PRESSURE (N/M««2) 
DENSITY(KG/M««3) 

E-H WIND U(M/SEC) 

N-S WIND V(M/8EC) 
WIND SPEED US(M/SEC) 
WIND DIRECTION WD 
(DEO) 

lERR 

I 

0 

ERROR FLAG 



CALLING SUBROUTINES t MODELS* MODADJ 
SUBROUTINES CALLED : DELTIM* JACROB 
COMMON BLOCKS USED : JRCOM* gSECOM* FILCOM 
EXTERNAL DATA SETS USED t OUT (DEBUG ONLY) 

ANALYST/PROGRAMMER : A. K. KAPOOR (MATH AND COMPUTING APPLICATION 

SECTION) COMPUTER SCIENCES CORPORATION 


Subroutine JR 
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PROGRAM LAIRS(OUTPUT»OUTjPERM 1 TWRKrUSEf KPyNOIDSf USEMET» 
lTAPE4»TAPE5rTAPE21 ,TAPE22f TAPE23» 

2TAPE6=0UTfTAPE 1=PERM1 rTAPE2=WRK>TAPE3=USE» 

3TAPE75---=KPf TAPE20 = N0IDS»TAPE30 = USEMET) 


VERSION 1.2 SEPTEMBER 1981 

purpose: 

LAIRS IS THE MAIN DRIVER FOR THE LANGLEY ATMOSPHERIC 
INFORMATION RETRIEVAL SYSTEM. 

method: 

1. WRITE TO DATA SET NOIDS THE VALUES OF THE COMMON 
BLOCKS SO THAT FOR THE STACKED CASES THE COMMON BLOCKS 
MAY BE RESTORED. 

2. CALL READIN TO PROCESS USER INPUT. 

3. CALL OUTPT TO PRINT THE TITLE PAGE AND RUN INFORMATION. 

4. CALL WFILE TO BUILD THE WORKING FILE. 

5. CALL SETREF TO SET REFERENCE VALUES FROM FILES. 

6. CALL MODADJ IF POLYNOMIAL AND/OR JACCHIA MODEL ADJUSTMENT 
IS REQUIRED. 

7. CALL PARAMS TO CALCULATE THE ATMOSPHERIC PROFILE. 

8. CALL OUTPTl TO CONCLUDE THE PRINTOUT. 

9. IF REQUIRED, CALL PREUSE TO CREATE A METEOROLOGICAL FILE 
FROM THE OUTPUT USER FILE. 

10. REWIND FILES AND INITIALIZE COMMON FOR STACKED DECK CASE. 

11. STOP. 

ARGUMENT LIST*. NONE 
CALLING subroutines: NONE 

SUBROUTINES CALLED*. RE AD I N , OUTPT , SETREF , MODADJ » PARAMS , 

OUTPTl, PREUSE 


Program LAIRS (1 of 2) 
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COMMON BLOCK PARAMETERS USEDt 


COMMON NAME 

PARAMETER 

I/O DESCRIPTION 

/USECOM/ 

MODL 

I 

LOWER ATMOSPHERE MODEL TYPE 


MODU 

I 

UPPER ATMOSPHERE MODEL TYPE 


ICRT 

I 

METEOROLOGICAL FILE CREATION FLAG 

/FILCOM/ 

LINES 

I 

LINE COUNTER 


NPRM 

I 

PERMANENT FILE LUN 


NURK 

I 

WORKING FILE LUN 


MWFl 

I 

METEOROLOGICAL FILE 1 LUN 


MWF2 

I 

METEOROLOGICAL FILE 2 LUN 


MWF3 

I 

METEOROLOGICAL FILE 3 LUN 


NENT 

I 

ENTREE <BET) FILE LUN 


NOIDS 

I 

SAVE FILE LUN 


NUSE 

I 

USER FILE LUN 

EXTERNAL DATA 

SETS USED? 



NAME 

LUN 

I/O OPERATIONS PERFORMED 

PERM 

NPRM 

REWIND 


WORK 

NURK 

REWIND 


METl 

MWFl 

REWIND 


MET2 

MWF2 

REWIND 


MET3 

MWF3 

REWIND 


ENT 

NENT 

REWIND 


USE 

NUSE 

REWIND 


NOIDS 

NOIDS 

WRITEfI 

READ » REWIND 

PROGRAMMER? D» 

E» BOLAND » 

CSC» 

MARCH 1981 
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SUBROUTINE L AP ( NZ » ZNEXT ) 

purpose: 

LAP GENERATES ALTITUDE POINTS ON THE LOCAL ATMOSPHERIC PROFILE. 

method: 

1. COMPUTE LOCAL SOLAR TIME OF PROFILE. 

2. CALCULATE ALL Z VALUES NEEDED. 

3. RETURN. 


argument list: 

ARGUMENT TYPE 

I/O 

DESCRIPTION 


nz i 

I 

NUMBER OF POINTS IN 

BLOCK 

ZNEXT I 

I/O 

NEXT Z VALUE TO BE 

COMPUTED 


CALLING subroutine: PARAMS 


SUBROUTINES CALLED: RLST 


COMMON BLOCK 

PARAMETERS USED 

* 


COMMON NAME 

PARAMETER 

I/O 

DESCRIPTION 

/USECOM/ 

THETA 

I 

LATITUDE 


TAU 

I 

LONGITUDE 


Z1 

I 

LOWER ALTITUDE FOR LAP 


Z2 

I 

UPPER ALTITUDE FOR LAP 


DZ 

I 

ALTITUDE INTERVAL 

/CALCOM/ 

TRJ(5rlO) 

0 

BLOCK OF TRAJECTORY POINTS 
Z>LAT.LONGf GMTfLST 


EXTERNAL DATA SETS USED: NONE 

programmer: D. E. BOLAND. CSC» march 1981 


Subroutine LAP 


3-49 



SUBROUTINE L AT AD J ( NF f F , A 1 » PI » A2 » P2 » DFL AT » T AU » RLAT ) 

purpose: 

LATADJ COMPUTES NEU LATITUDE GRADIENTS FOR AMY INPUT 

parameter f (t»p» or D) from meteorological data, 
method: 

1. TRANSLATE EACH INPUT PARAMETER F (UP TO 3 POINTS — 

ONE FROM EACH MET FILE) FROM ITS LOCAL SOLAR TIME TAU 
TO THE MODEL REFERENCE LOCAL SOLAR TIME REFTAU TO 
REMOVE THE DIURNAL/SEMIDIURNAL COMPONENT OF THE 
VARIATION OF F. 

2. COMPUTE THE LATITUDE GRADIENT OF F FROM THE TRANSLATED 
values of F by THE METHOD OF LEAST SQUARES, LATITUDE 
UNITS ARE DEGREES. 

3. RETURN. 


ARGUMENT 

list: 



ARGUMENT 

TYPE 

I/O 

description 

NF 

R 

I 

NUMBER OF POINTS OF F 

F 

ft 

I 

parameter TO BE TRANSLATED 

A1 

R 

I 

DIURNAL AMPLITUDE OF F 

PI 

R 

I 

diurnal phase of F 

A2 

R 

X 

SEMIDIURNAL AMPLITUDE OF F 

P2 

R 

I 

SEMIDIURNAL PHASE OF F 

TAU 

R 

1 

LOr Al. SOLAR TIME OF MET PROFILE 

DFLAT 

R 

0 

LATITUDE GRADIENT OF F 

RLAT 

R 

1 

LATITUDE OF MET PROFILE 


CALLING bUBROUriNE: TRNADJ 

SUBROUTINES CALLED! TRNDSD 

COMMON BLOCK PARAMETERS USED! 

COMMON NAME parameter I/O DESCRIPTION 
/USECOM/ reflat I MODEL REFERENCE LATITUDE 

REFTAU 1 MODEL REFERENCE LOCAL SOLAR TIME 

lIiSDA 1 FLAG FOR D/SD ADJUSTMENT 

EXTERNAL DATA SETS USED: NONE 

programmer: D. E. B0LAND» CSCf march 1981 


Subroutine LATADJ 
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SUBROUTINE LOWALT ( RHO > T ■/ P ) 


PURPOSE 

LOUALT IS CALLED BY JACROE TO FURNISH ATMOSPHERIC DENSITY VALUES 
AT AND BELOW 1.25 KM. 

THIS SUBROUTINE SELECTS BETWEEN SUBROUTINE BARODE AND DIFFDE 
FOR COMPUTING THE LiNCORRECTED DENSITY^ DEPENDING ON WHETHER 
THE SATELLITE IS BELOW OR ABOVE lOOKM. 


METHOD ; 

1. EVALUATE CONSTANTS 

2. CALL SUBROUTINE ROOTS TO COMPUTE THE ROOTS R1»R2» 
X+IY, AND X-IY OF THE POLYNOMIAL 

3. COMPUTE THE CONSTANTS NEEDED BY BARODE AND DIFFDE 

4. COMPUTE THE TEMPERATURE AT THE SATELLITE 

5. IF 7-<100KMt call BARODE? 

IF Z>iOOKM» COMPUTE THE TEMPERATURE AT lOOKH AND 
CALL DIFFDE. 


ARGUMENT LIST? 

ARGUMENT TYPE I/O 
RHQ R 0 
T R 0 
P R D 


DESCRIPTION 

RHO IS A PARTIAL ATMOSPHERIC DENSITY 

TEMPERATURE 

PRESSURE 


CALLING subroutine; JACROB 


SUBROUTINES CALLED: 




BARODE 

COMPUTE RHO 

BE 

TWEEN 90 AND 100 KM 


DIFFDE 

COMPUTE RHO 

BETWEEN 100 AND 125 KM 


ROOTS 

TO COMPUTE 

THE 

ROOTS OF THE POLYNOMIAL 

COMMON BLOCK 

VARIABL ES 



NAME 

in MEN 

COMMON 

I/O 

DESCRIPTION 

CC 

5 

JRCOM 

T 

POWER SERIES COEFFICIENTS FOR T7 

DUHl 

1 

JRCOM 

I/O 

IMMAGINARY PART OF RLl 

DUM2 

1, 

JRCOM 

I/O 

IMMAGINARY PART OF RL2 

FKL 

1. 

JRCOM 

I 

FACTOR INVOLVED IN RHO COMPUTATION 

FLC4 

1 

JRCOM 

i: 

MODIFYING FACTOR- 

GLO. 

1 

JRCOM 

I 

MEAN SURFACE GRAVITY 

HGT 

1 

JRCOM 

I 

HEIGHT OF SPACECRAFT 

RC 

1 

JRCOM 

1 

UNIVERSAL GAS CONSTANT 

RCM 

1 

JRCOM 

1 

AVERAGE EARTH RADIUS 

RLl 

1. 

JRCOM 

1 

ROOT OF POLYNOMIAL IN INTEGRAND 

RL2 

;l 

JRCOM 

I 

ROOT OF POLFNOMIAL IN INTEGRAND 

TCIL 

1 

JRCOM 

I 

TEMPERATURE AT 100 KM. 

TX 

:i. 

JRCOM 

I 

INFLECTION POINT TEMPERATURE 

T7. 

1 

JRCOM 

I 

TEMPERATURE AT HEIGHT 7 

TO 

1 

JRCOM 

I 

TEMPERATURE AT MINIMUM HEIGHT 

UC 

9 

JRCOM 

I 

FUNCTIONAL VALUES AT RLl AN-D RL2 

VCD I 

.1 

JRCOM 

I 

FACTOR INCLUDEDIN RHO COMPUTATION 

WC 

2 

JRCOM 

I 

FUNCTIONAL VALUES AT RLl AND RL2 

XCDI 

1 

JRCOM 

I 

FACTOR INCLUDEDIN RHO COMPUTATION 

XLPS 

1 

JRCOM 

I 

ROOT OF POLYNOMIAL IN INTEGRAND 

YLPS 

1 

JRCOM 

I 

ROOT OF POLYNOMIAL IN INTEGRAND 

PROGRAMMER 




A, 

K. KAPOOR T COMPUTER SCIENCES CORPORATION 


Subroutine LOWALT 
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SUBROUTINE HETINT ( NH » POS r TPDUN » TPDUHU > lERR ) 


purpose: 

METINT IS THE EXECUTIVE SUBROUTINE FOR CONTROLING METEOROLO- 
GICAL INTERPOLATION FOR TEHPERATUREr PRESSURE* DENSITY* AND 
WIND COMPONENTS. 

method: 

1. CALL METINl TO OBTAIN INTERPOLATED VALUES FOR T* P* D* 

U* AND V IN ALTITUDE FOR EACH MET. PROFILE PROVIDED. 

2. PERFORM UNIVARIATE OR BIVARIATE INTERPOLATION FOR T* P* D* U* 
AND V BETWEEN THE TWO OR THREE MET. PROFILES PROVIDED. 

3. RETURN. 

ARGUMENT LIST: 


ARGUMENT 

TYPE 

I/O 

DESCRIPTION 

NH 

I 

I 

ATMOSPHERIC SEGMENT* 

«0* FLAG FOR CALL FROM SETREF 
=1* UPPER SEGMENT 
»2* LOWER SEGMENT 

POS 

R 

I 

4-WORD ARRAY FOR TRAJECTORY POINT: 
P08(1)-ALT1TUDE 
P0S(2)»LATITUDE 
P0S<3>»L0CAL SOLAR TIME 
P08<4)*L0NGITUDE 

TPDWM 

R 

0 

CALCULATED ATMOSPHERIC PARAMETERS: 
TPDWM<1>«TEMPERATURE 
TPDWH<2)«PRESSURE 
TPDWM(3)»DENSITY 
TP0WM(4>«E-W WIND * U 
TPDWN<5)»N-S WIND * V 
TPDWM<6)nWIND SPEED * W8 
TPDWM<7)«WIND DIRECTION * WD 

TPDWMU 

R 

0 

UNCERTAINTIES IN PARAMETERS (7) 

lERR 

I 

0 

ERROR flag: 

-0* NO ERROR 
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CALLING subroutines: HODELS 

SUBROUTINES CALLED: METIN1>DFLINT» IBIt TRNLATf TRNDSD 


COMMON BLOCK 

PARAMETERS USED 

• 

¥ 

COMMON NAME 

PARAMETER 

I/O 

/FILCOM/ 

NPRINT 

I 

/USECOM/ 

INTP 

I 


INT 

I 



NFILE 

I 


BOUND 

I 


RMGMT 

I 

/INTCOM/ 

II 

I 


ARR 

I 

/NETCOM/ 

RLAT 

I 


RLONO 

I 


RL8TM 

I 


DESCRIPTION 
PRINTER OUTPUT LUN 
ALTITUDE INTERPOLATION USED 
METEOROLOGICAL INTERPOLATION USED 
UNIVARIATE 

»lr 1ST ORDER BIVARIATE 
<^ 2 , 2ND ORDER BIVARIATE 
NUMBER OF MET* PROFILES 
BOUNDARY BETWEEN UPPER AND 
LOWER ATMOSPHERIC SEGMENTS 
GMT FOR EACH MET* PROFILE 
POINTER USED IN ARRAY ARR 
ARRAY OF COEFFICIENTSr INCLUDING 
LATITUDE GRADIENTS 
ARRAY FOR LATITUDES OF MET* FILES 
ARRAY FOR LONGITUDES OF MET FILES 
LST FOR EACH MET* PROFILE 


EXTERNAL DATA SETS USED: 

NAME LUN I/O OPERATIONS PERFORMED 

OUT NPRINT WRITE 

programmer: R*A* KUSESKIf computer SCIENCES CORP* 


Subroutine METINT (2 of 2) 
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SUBROUTINE METINl ( ALTf TPDW»TPDWhUf NFILRDf lERRl ) 

purpose: 

METINl READS THE GIVEN HETEOROLOGICAL PROFILES AND INTERPOLATES 
THEM IN ALTITUDE FOR Jr P» D, U» V* 

method: 

1. READS EACH MET. PROFILE IN SEQUENCE AND STORES VALUES IN 
ARRAYS. 

2. CALLS lUNI TO INTERPOLATE ARRAYS IN ALTITUDE. 

3. USES LOGARITHMIC INTERPOLATION FOR PRESSURE AND DENSITY. 

4. DETERMINES LARGEST UNCERTAINTY FOR EACH PARAMETER. 

5. RETURNS. 

ARGUMENT LIST’. 


ARGUMENT 

TYPE 

I/O 

DESCRIPTION 

ALT 

R 

I 

TRAJECTORY ATHTUDE 

TPDW 

R 

0 

5-WORD ARRAY WITH T»P»D»UfV 
FOR EACH PROFILE 

TPDWMU 

R 

0 

UNCERTAINTIES IN TrP»DfUfV 

NFILRD 

I 

0 

FLAG FOR READING MET. PROFILE 
=0» FILE READ 

■=1t EOF OR EXTRAPOLATION IN 

lERRl 

I 

0 

ERROR FLAG 

CALLING subroutines: 

METINTr TRNADJ 


SUBROUTINES CALLED: lUNI 


COMMON BLOCK PARAMETERS USED: 


COMMON NAME 

PARAMETER I/O 

DESCRIPTION 

/USECOM/ 

INPT I 

ALTITUDE INTERPOLATION FLAG 
1ST ORDER LAGRANGE 
=2f 2ND ORDER LAGRANGE 


IDBG I 

DEBUG FLAG 


NFILE I 

NUMBER OF MET. PROFILES 

/FILCOM/ 

NPRINT I 

PRINTED OUTPUT LUN 


METLUN I 

FIRST MET. PROFILE INPUT LUN 

/METCOM/ 

IPT I 

POINTER USED IN lUNI 


METFLG I 

FLAG FOR 1ST CALL TO METINl 


II I 

POINTER TO CLOSEST ALTITUDE 
IN DATA ARRAR (METARR) 


ALTM I 

ARRAY OF ALTITUDES FOR EACH 
PROFILE 


METARR I 

ARRAY OF T»P»D» AND WIND SPEED 
AND DIRECTION FOR EACH FILE 


FSTPT I 

LOWEST ALTITUDE FOR EACH FILE 


RLAT I 

LATITUDE OF EACH PROFILE 


RLONG I 

LONGITUDE OF EACH PROFILE 

EXTERNAL DATA 

SETS used: 


NAME 

LUN 

I/O OPERATIONS PERFORMED 

METl 

METLUN 

READ 

MET2 

METLUN+1 

READ 

MET3 

METLUN+2 

READ 

OUT 

NPRINT 

WRITE 

PROGRAMMER*. R 

.A. KUSESKIy COMPUTER SCIENCES CORP. 


Subroutine METINl 
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SUBROUTINE METREF (I SEG j NMET » Z » T » P » D » UT » UP » UD » I ERR ) 


purpose: 

METREF PRODUCES A SINGLE REFERENCE METEOROLOGICAL PROFILE 
<1 SEGMENT AT A TIME) FROM THE INPUT METEOROLOGICAL PROFILES* 

method; 

1, READ THE POINTS IN SEGMENT ISEG OF THE NEXT PROFILE, 

2, CALCULATE THE WEIGHTING FACTORS FOR EACH POINT, 

3, IF TRANSLATION IS REOUIREDf CALL DFLINT TO OBTAIN 
TRANSLATIONAL COEFFICIENTS AMD CALL TRNLAT AND TRNDSD 
TO PERFORM TRANSLATION, 

4, CALCULATE THE LOG BASE E OF PRESSUF:;E AND DENSITY, 

5, REWIND THIS METEOROLOGICAL PROFILE, 

6, REPEAT STEPS 1--5 FOR THE NEXT PROFILE UNLESS ALL PROFILES 
HAUE BEEN PROCESSED, 


7. RETURN 

♦ 




ARGUMENT LIST 

♦ 




ARGUMENT 

TYPE 

I/O 

DESCRIPTION 


ISEG 

I 

:l' 

ATMOSPHERIC SEGMENT 

NUMBER 

NMET 

I 

I 

NUMBER OF POINTS IN 

THIS SEGMENT 

Z 

R 

I 

ALTITUDE ARRAY 


T 

R 

I 

TEMPERATURE ARRAY 


P 

R 

I 

PRESSURE ARRAY 


D 

R 

I 

DENSITY ARRAY 


WT 

R 

I 

TEMPERATURE WEIGHTS 


UP 

R 

I 

PRESSURE WEIGHTS 


UD 

R 

I 

DENSITY WEIGHTS 


lERR 

I 

□ 

ERROR FLAG 



CALLING subroutine: POLADJ 
SUBROUTINES CALLED: DFLINT 


Subroutine METREF (1 of 2) 



COMMON BLOCK 

PARAMETERS USED 

COMMON NAME 

PARAMETER 

/USECOM/ 

ILATL 


IDSDL 


IWL 


NFIL 


BOUND 


RMERGE 


REFLAT 


REFTALI 


RM6MT 

/INTCOM/ 

12 

/FILCOM/ 

NPRINT 


WRK 


MWF 

/DCCOM/ 

WTDCT 


WTDCP 


WTDCD 

EXTERNAL DATA 

SETS used: 

NAME 

LUN 

OUT 

NPRINT 

WRK 

NWRK 

MET I 

MWF( 1 ) 

MET2 

MWF<2) 

MET3 

MWF< 3) 


I/Q DESCRIPTION 

I FLAG FOR LATITUDE VARIATIONS 

I FLAG FOR D/SD VARIATIONS 

I WIND FLAG 

I NUMBER OF MET FILES 

I 3 SEGMENT BOUNDARIES 

I REGION FOR SEGMENT MERGIMOi 

I REFERENCE LATITUDE 

I REFERENCE LOCAL SOLAR TIME 

I GMT OF MET FILES 

I FLAG TO REINITIALIZE DFLINT 

I PRINTED OUTPUT LUN 

I WORKING FILE LUN 

I METEOROLOGICAL FILES LUNS 

I A PRIORI TEMPERATURE WEIGHTING 

I A PRIORI PRESSURE WEIGHTING 

I A PRIORI DENSITY WEIGHTING 


I/O OPERATIONS PERFORMED 
WRITE (DEBUG ONLY) 
REWIND 

READ» REWIND 
READ* REWIND 
READ. REWIND 


programmer: 


D* E» BOLAND. CSC 


Subroutine METREF (2 of 2) 
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SUBROUTINE MODADJ(IERR) 


purpose: 

HODADJ IS THE MODEL ADJUSTMENT EXECUTIVE, 

method: 

1. CALL POLADJ IF POLYNOMIAL MODEL ADJUSTMENT IS REQUIRED. 

2. RESET JACCHIA-R0BERT8 MODEL DENSITY BOUNDARY VALUES AT 

90f lOOf AND 125 KILOMETERS. (NOTE THIS PROCEDURE REPLACES 
THE FIRST 10 WORDS OF THE ARRAY DD<35) IN JRCOM WITH THE 
ARRAYS FD125(5) AND D125(5J). 

3. RETURN. 

ARGUMENT LIST: 

ARGUMENT TYPE I/O DESCRIPTION 
lERR I 0 ERROR FLAG: 

=^ 0 * NO ERROR 

CALLING subroutine: LAIRS 
SUBROUTINES CALLED: POLADJf JR 


COMMON BLOCK 

PARAMETERS USED 

• 

• 


COMMON NAME 

PARAMETER 

I/O 

DESCRIPTION 

/USECOM/ 

MODL 

I 

TYPE OF LOWER ATMOSPHERIC MODEL 


MODU 

I 

TYPE OF UPPER ATMOSPHERIC MODEL 


IDB6 

I 

DEBUG FLAG 


BOUNDl 

I 

TOP SEGMENT BOUNDARY (KM) 

/JRCOM/ 

FD125 

I 

INDIVIDUAL NUMBER DENSITIES 




AT 125 KM 


D125 

0 

LOO INDIVIDUAL NUMBER DENSITIES 




AT 125 KM 


AVG 

I 

AVOGADRO'S NUMBER 


DBOUND 

I/O 

DENSITY VALUE AT BOUNDl 


ZD 

I/O 

COEFFICIENTS FOR DENSITY AT 100 

EXTERNAL DATA 

SETS used: 



NAME 

LUN 


I/O OPERATIONS PERFORMED 

OUT 

NPRINT 


WRITE (DEBUG ONLY) 

programmer: d 

. E. BOLAND r CSC 



Subroutine MODADJ 



SUBROUTINE MODELS < POS j UPOS y TPDU » UTPDW » lERR ) 


PURPOSE t 

MODELS CONTROLS THE COMPUTATION OF TEMPERATURE 1? PRESSURE Pr 
DENSITY D» AND WINDS U» V AND WSf WD FOR EACH POINT ON 
A SHUTTLE TRAJECTORY OR A LOCAL ATMOSPHERIC PROFILE* 

method: 

1* determine if the current point is in THE UPPER OR THE 
LOWER ATMOSPHERE* 

2* BRANCH TO THE APPROPRIATE CALCULATION ROUTINE ON THE. BASIS 
OF FLAGS BET BY THE USER* 

3* COMPUTE WINDS FOR NON- I NTERPOL AT I ON MODELS. 

4* RETURN* 


ARGUMENT LIST 

♦ 



ARGUMENT 

TYPE 

I/O 

DESCRIPTION 

POS 

R 

I 

TRAJECTORY POINT: 

POS< 1 )=ALTITUDE 
P0S<2)=LATITUDE 
P0S(3)=L0CAL SOLAR TIME 
P0S<4)=L0NGITUDE 

UPOS 

r 

I 

UNCERTAINTIES IN TRAJECTORY POINT (3) 

TPDU 

R 

0 

CALCULATED ATMOSPHERIC PARAMETERS: 
TPDW(1)=TEMPERATURE 


TPDW<2)=PRESSURE 
TPDW(3)=DENSITY 
TPDW<4)=E-W WIND U 
TPDW(5)=N-S WIND U 
TPDW(6)=WIND SPEED US 
TPDW<7)=WIND DIRECTION WD 

UTPDW R 0 UNCERTAINTIES IN PARAMETERS (7) 

lERR I 0 ERROR FLAGt 

=0» NO ERROR 

CALLING subroutine: PARAMS 

SUBROUTINES CALLED: DFLINT » METINT » JR y POL Y » WIND » WCNURT 

COMMON BLOCK PARAMETERS USED: 

COMMON NAME PARAMETER I/O DESCRIPTION 
/USECOM/ MODL I LOWER ATMOSPHERE MODEL TYPE 

MODU I UPPER ATMOSPHERE MODEL TYPE 

IDBG I DEBUG FLAG 

BOUND I UPPER-LOWER ATMOSPHERE BOUNDARY 

/FILCOM/ NPRINT I PRINTED OUTPUT LUN 

EXTERNAL DATA SETS USED: 

NAME LUN I/O OPERATIONS PERFORMED 

OUT NPRINT WRITE(DEBUG ONLY) 

programmer: D* E* B0LAND» esc. march 1981 


Subroutine MODELS 
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SUBROUTINE NLSPOL ( I SEG » KS » NMET » ZMET > TMET » DMET r PMET » TW2 » DW2 » PW2 » 
♦ZBOUNDrPBOUNDfSUMrE^Bf lERR) 

PURPOSE 

NLSPOL IS THE EXECUTIVE OF THE NON-LINEAR LEAST-SQUARE FITTING 
PROCEDURE 

METHOD? 

THE SUM OF THE WEIGHTED SQUARED RESIDUALS IS MINIMIZED ITERATIVELY, 

1, COMPUTE WEIGHTED RESIDUALS 

2, COMPUTE TEMPERATURE»PRESSURE» AND DENSITY PARTIAL DERIVATIVES 

3, CONSTRUCT NORMAL MATRIX 

4, INVERT THE NORMAL MATRIX AND COMPUTE THE CORRECTIONS TO -THE 
SOLVE-FOR PARAMETERS 

5, RETURN 

ARGUMENT LIST 

ARGUMENT TYPE I/O DESCRIPTION 

ISEG I I SEGMENT NUMBER 

KS I I CONTROLS THE NUMBER OF SOLVE-FOR PARAMETERS 

=1 SOLVES FOR ALL TEMP, COEFS AND PBOUND 
=2 SOLVES FOR TEMP. COEFS WITH BOUNDARY TEMP. 
AND PBOUND FIXED 

=3 SOLVES FOR ALL TEMP. COEFS WITH PBOUND 
FIXED 

=4 SOLVES FOR ALL PARAMETERS WITH BOUNDARY 
TEMP. FIXED 


NMET 

I 

I 

NUMBER OF MET. ALTITUDE POINTS 

ZMET 

R 

I 

ALTITUDE ARRAY OF MET. DATA 

TMET 

R 

I 

MET. TEMPERATURE DATA 

DMET 

R 

I 

LOG OF MET. DENSITY 

PMET 

R 

I 

LOG OF MET. PRESSURE 

TW2 

R 

I 

SQUARED WEIGHT FOR TEMPERATURE 

DW2 

R 

I 

SQUARED WEIGHT FOR LOG DENSITY 

PW2 

R 

I 

SQUARED WEIGHT FOR LOG PRESSURE 

ZBOUND 

R 

I 

BOUNDARY ALTITUDE 

PBOUND 

R 

I 

PRESSURE AT ZBOUND 

SUM 

R 

0 

SUM OF SQUARED WEIGHTED RESIDUALS 

E 

R 

0 

VARIANCE-COVARIANCE MATRIX 

B 

R 

I/O 

TEMPERATURE COEFS 

lERR 

I 

0 

ERROR FLAG 


CALLING subroutines: POLFIT 


Subroutine NLSPOL (1 of 2) 
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SUBROUTINES CALLED: TPAR» DPAR» PPAR> SOLVE? SIMP? OUTRES? OUTDC 


COMMON BLOCK 

PARAMETERS 

used: 


COMMON NAME 

PARAMETER 

I/O 

DESCRIPTION 

/USECOM/ 

ITRMAX 

I 

MAX* NUMBER OF ITERATIONS 


KTMET 

I 

=0 USE TEMPERATURE DATA 

=1 DO NOT USE TEMPERATURE DATA 


KPMET 

I 

=0 USE PRESSURE DATA 

=1 DO NOT USE PRESSURE DATA 


KDMET 

I 

=0 USE DENSITY DATA 

=1 DO NOT USE DENSITY DATA 


RSTEP 

I 

STEP SIZE FOR INTEGRATION 


EPSCON 

I 

CONVERGENCE CRITERIA 


RMSZRO 

I 

INITIAL RMS 


RMSADD 

I 

RMS ADDER 


FSIGMA 

I 

MULTIPLYING FACTOR TO RMS FOR 
EDITTING 

/SOLAND/ 

h i NT 

I 

NUMBER OF INTEGRALS 


KSOLA4 

I 

-0 NOT SOLVES FOR BOUNDARY TEMP. 
=1 SOLVES FOR BOUNDARY TEMP. 


KSOLAJ; 

I 

=0 NOT SOLVES FOR BOUNDARY PRESSURE 
-1 SOLVES FOR BOUNDARY PRESSURE 


NPOL 

I 

ORDER OF POLYNOMIALS FOR TEMP. 


MAXSOL 

I 

MAX. NUMBER OF SOLVE-FOR 
PARAMETERS 


NSOLV 

0 

NUMBER OF SOLVE-FOR PARAMETERS 


KEDTOT 

0 

TOTAL NUMBER OF MET. DATA USED 


K0NVR6 

0 

CONVERGENCE FLAG 


ZREF 

I 

SEGMENT BOUNDARY (ZBOUND) 


EPSLON 

I 

CONVERGENCE CRITERIA (EPSCON) 


SRSUM 

0 

TOTAL SUM OF SQUARED WEIGHTED 
RESIDUALS 


CURRMS 

0 

CURRENT RMS 


OLDRMS 

0 

PREVIOUS RMS 


RMSP 

0 

PREDICTED RMS 


A 

I/O 

SOI VE-FOR PARAMETER ARRAY 


6 

0 

RIGHT HAND SIDE OF NORMAL EQU. 


SH 

0 

LINEAR ARRAY OF NORMAL MATRIX 

/JRCOM/ 

RE 

I 

MEAN RADIUS OF THE EARTH 


RSTAR 

I 

UNIVERSAL GAS CONSTANT 


GS 

I 

GRAVITATIONAL ACC. AT 
THE SURFACE OF THE EARTH 


XM 

I 

MEAN MOLECULAR WEIGHT 

EXTERNAL DATA 

SETS used: 



NAME 

LUN 

I/O OPERATIONS PERFORMED 

OUT 

NPRINT 

WRITE 

(DEBUG ONLY) 

programmer: t» 

LEE? CSC 

MARCH 

1981 
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SUBROUTINE OUTDC< ISEG» ITER» NMETf KTT»KDT»KPT»K0NVR6f KB0LA4» 
*KS0LA5 » CURRMS » OLDRMS » RMSP » AP A » A » DEL » UCMAT ) 

purpose: 

OUTDO PRODUCES A DIFFERENTIAL CORRECTION RUN SUMMARY 


ARGUMENT 

ARGUMENT 

list: 

TYPE 

I/O 

DESCRIPTION 

ISEG 

I 

I 

=2 UPPER SEGMENT 

ITER 

I 

I 

=3 MIDDLE SEGMENT 
=4 LOWER SEGMENT 
ITERATION NUMBER 

NMET 

I 

I 

NUMBER OF METEOROLOGICAL DATA 

KTT 

I 

I 

TOTAL NUMBER OF TEMP» DATA USED 

KDT 

I 

I 

TOTAL NUMBER OF DENSITY DATA USED 

KPT 

I 

I 

TOTAL NUMBER OF PRESSURE DATA USED 

KONURG 

I 

I 

CONVERGENCE FLAG 

KS0LA4 

I 

I 

=1 CONVERGED 
=2 CONVERGING 
=3 DIVERGING 

=4 CONVERGED: LINEAR FITTING 
=0 DOES NOT SOLVE FOR BOUNDARY TEMP. 

KS0LA5 

1 

1 

=0 DOES NOT SOLVE FOR BOUNDARY PRESSURE 

CURRMS 

R 

I 

CURRENT WEIGHTED RMS 

OLDRMS 

R 

I 

PREVIOUS WEIGHTED RMS 

RMSP 

R 

I 

PREDICTED RMS 

APA 

R 

I 

A PRIORI SOLVE-FOR PARAMETERS 

A 

R 

I 

CURRENT VALUES OF SOLVE-FOR PARAMETERS 

DEL 

R 

I 

CORRECTION : CURRENT - PREVIOUS 

VCMAT 

R 

I 

VARIANCE-COVARIANCE MATRIX 


CALLING subroutine: NLSPOL 
SUBROUTINES CALLED: NONE 


COMMON BLOCK PARAMETERS USED: 
COMMON NAME PARAMETER I/O 
/FILCOM/ NPRINT I 

IPAGE I 


DESCRIPTION 
PRINTED OUTPUT LUN 
PAGE COUNTER 


EXTERNAL DATA SETS USED: 
NAME LUN 

OUT NPRINT 


I/O OPERATIONS PERFORMED 
WRITE 


programmer: 


A. K* KAPOOR r T* LEE» CSC 


MARCH 1981 


Subroutine OUTDC 
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SUBROUTINE OUTPRF< IPRFLrNZ) 


purpose: 

SUBROUTINE OUTPRF GENERATES A BLOCK PRINTOUT OF THE 
ATMOSPHERIC PARAMETERS* 

METHOD f 

1* PRINT HEADER AND COLUMN HEADINGS* 

2* PRINT BLOCK OF DATA <10 TRAJECTORY POINTS)* 

3* WRITE USER FILE* 

4* RETURN* 

ARGUMENT LIST : 

ARGUMENT TYPE I/O DESCRIPTION 

IPRFL I I TYPE OF DATA» 1=ENTREE » 2=LAP 

NZ II NUMBER OF POINTS IN BLOCK 

CALLING SUBROUTINES : PARAMS 

SUBROUTINES CALLED : NONE 

COMMON BLOCK PARAMETERS USED : 

COMMON NAME PARAMETER I/O DESCRIPTION 
/CALCOM/ TRJ I 10 TRAJECTORY POINTS: 

ALTiLATfLONGfGMTr AND LST 
UTRJ I UNCERTAINTIES IN ALT»LAT»LONO» 

OMTf AND LSTI 

TPDW I 7-UORD array: T»PrDt Ut MtUS»AND UD 

UTPDU I UNCERTAINTIES IN T>P»Di Uf VrWS» 

AND WD 

WTMOL 1 MEAN MOLECULAR WEIGHT (10) 

SCHT I PRESSURE SCALE HEIGHT (10) 

RMACH 1 MACH NUMBER (10) 

OBPRES I ON-BOARD PRESSURE (10) 

PCDEU I PERCENT DEVIATION FROM GAS LAW(IO) 

/FILCOM/ NPRINT I PRINTED OUTPUT LUN 

NUSE 1 USER FILE LUN 

LINES I/O LINE COUNTER 

IPAGE I/O NUMBER OF PAGES PRINTED 

EXTERNAL DATA SETS USED : 

NAME LUN I/O OPERATIONS PERFORMED 

OUT NPRINT WRITE 

USE NUSE WRITE 

PJ^OGRAMMER : R*A* KUSE8KI CSC> MARCK 1981 


Subroutine OUTPRF 
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SUBROUTINE OUTPT(IERR) 


purpose: 

OUTPT WRITES THE INITIAL CONDITIONS REPORT 

method: 

1» PRINTS HEADER AND OPTIONS CHOSEN FOR THE LOWER AND UPPER 
ATMOSPHERE. 

2. PRINTS DESCRIPTION OF OUTPUT PARAMETERS WITH UNITS. 

3. RETURN. 


ARGUMENT LIST: 

ARGUMENT TYPE 

I/O 

DESCRIPTION 

lERR I 

I 

ERROR FLAG 

CALLING subroutine: 

LAIRS 


SUBROUTINES CALLED: 

NONE 



COMMON BLOCK 

parameters used: 


COMMON NAME 

parameter 

I/O 

DESCRIPTION 

/USECOM/ 

IPRFL 

I 

TYPE OF PROFILE 


modl 

I 

TYPE OF LOWER ATMOSPHERIC MODEL 


IWL 

I 

TYPE OF L.A. WIND MODEL 


ILATL 

I 

LATITUDE VARIATIONS IN L.A. 


idsdl 

I 

DIURNAL/SEMIDIURNAL VARIATIONSr L.A 


MODU 

I 

TYPE OF UPPER ATMOSHPHERIC MODEL 


IWU 

I 

TYPE OF U.A. WIND MODEL 


ILATU 

I 

LATITUDE VARIATIONS IN U.A. 


IDSDU 

I 

diurnal/semidiurnal VARIATIONS^ U.A 


IPHYS 

I 

DERIVE PRESSURE AND DENSITY FROM 
TEMPERATURE VIA PHYSICAL LAWS 
(OVERRIDES ALL OTHER SELECTIONS) 


INTP 

I 

INTERPOLATOR TYPE - ALTITUDE ONLY 


IDBG 

I 

PRINT DEBUG 


INT 

I 

INTERPOLATER TYPE - BETWEEN 
MET. PROFILES 


IDNUM 

I 

ENTREE-LAP ID NUMBER 


IGAS 

I 

GAS LAW CALCULATION FLAG 


IDC 

I 

DIFFERENTIAL CORRECTION FLAG 


NPTSE 

I 

NUMBER OF POINTS IN ENTREE FILE 


NSAMP 

I 

SAMPLE RATE FOR ENTREE FILE 


YMD 

I 

YEARf MONTH» DAY OF LAP OR ENTREE 
TRAJECTORY 

/FILCOM/ 

NPRINT 

I 

PRINTED OUTPUT LUN 


NUSE 

I 

USER FILE LUN 

EXTERNAL DATA 

IPAGE 
SETS used: 

I/O 

NUMBER OF PAGES PRINTED 

NAME 

LUN 


I/O OPERATIONS PERFORMED 

OUT 

NPRINT 


WRITE 

USE 

NUSE 


WRITE 

programmer: 

R.A. KUSESKIf 

CSC» 

MARCH 1981 


Subroutine OUTPT 
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SUBROUTINE OUTPTJ. (lERR) 


PURPOSE? 

OUTPTl WRITES THE SECOND HALF OF THE OUTPUT PARAMETERS REPORT 
AND THE ERROR SUMMARY INFORMATION* 

METHOD? 

1. READS USER FILE AND PRINTS DESIRED INFORMATION 

2, WRITES ERROR RETURN SUMMARY INFORMATION 
3* RETURNS 

ARGUMENT LIST? 

ARGUMENT LIST? 

ARGUMENT TYPE I/O DESCRIPTION 
lERR I I ERROR FLAG 

CALLING SUBROUTINE? LAIRS 

SUBROUTINES CALLED? NONE 

COMMON BLOCK PARAMETERS USED? 


COMMON NAME 

PARAMETER 

I/O 

DESCRIPTION 

/FILCOM/ 

NPRINT 

I 

PRINTER OUTPUT LUN 


NUSE 

I 

USER FILE LUN 


LINES 

I/O 

LINE COUNTER 


IPAGE 

I/O 

PAGE COUNTER 

EXTERNAL DATA 

SETS USED? 



NAME 

LUN 

I/O OPERATIONS PERFORMED 

OUT 

NPRINT 

WRITE 


USE 

NUSE 

READ 



PROGRAMMER? R*A. KUSESKI» COMPUTER SCIENCES CORPORATION^ JULY 1981 


Subroutine OUTPTl 



SUBROUTINE OUTRES ( I SEG » I 1 » I 2 f NMET » Z » Y f WT » RES I D » SUM » lEUI T ) 

purpose: 

TO OUTPUT A RESIDUALS REPORT FOR A LAIRS LINEAR OR NONLINEAR 
PROGRAM RUN 

method: 

ALL THE INFORMATION THAT HAS TO BE REPORTED IS PASSED THROUGH 
THE ARGUMENT LIST* THE INFORMATION INVOLVED IS THE SEGMENT 
NUMBERt what type of data is BEING OUTPUT ( TEMPERATURE » 
PRESSURE^ OR DENSITY) THE ITERATION NUMBER^ NUMBER OF POINTS^ 
ALTITUDE^ OBSERVED VALUES^ WEIGHT OF RESIDUALS» RMSf AND 
THE EDITING CRITERIA FLAGS* 

THIS INFORMATION IS PARTITIONED IN ELEVEN COLUMNS AND OUTPUT 
IN A UNIFORM FORMAT* 

ARGUMENT LIST: 


VARIABLE 

TYPE 

DIMENSION 

DESCRIPTION 

ISEG 

I 

1 

SEGMENT NUMBER FLAG I*E*» 
UPPER» MIDDLE. OR LOWER 

11 

I 

1 

TEMPERATURE. PRESSURE. OR DENSITY 
FLAG 

12 

I 

1 

RUN TYPE AND ITERATION NUMBER FLAG 
=0. LINEAR RUN AND ONLY ONE REPORT 
IS PRINTED 

= ITER * . NONLINEAR RUN AND THIS 
FLAG IS THE ITERATION NUMBER 

NMET 

I 

NMET 

NUMBER OF POINTS 

Z 

R 

Z(NMET) 

ALTITUDE RELATED TO N TH POINT 

Y 

R 

Y(NMET) 

OBSERVED VALUE OF TEMPERATURE. 
PRESSURE. OR DENSITY 

WT 

R 

WT(NMET) 

WEIGHTS ASSOCIATED WITH N TH POINT 

RES ID 

R 

(NMET) 

RESIDUALS FOR THE N TH POINT 

SUM 

R 

1 

THE RMS 

lEDIT 

I 

(NMET) 

EDITING FLAG FOR THE N TH POINT 


CALLING routines: NLSPOL » POLF I T 


SUBROUTINES CALLED: NONE 


COMMON BLOCK PARAMETERS : 

COMMON NAME PARAMETER I/O 

/FILCOM/ NOUT I 

IPREV I 

NUMLIN I 

IPA6E I 


DESCRIPTION 
PRINTED OUTPUT LUN 
ITERATION NO. OF LAST ITERATION 
LINE COUNTER 
PAGE COUNTER 


EXTERNAL DATA SETS USED 
NAME LUN 

OUT NOUT 


I/O OPERATIONS PERFORMED 
WRITE 


PROGRAMMER/ANALYST : A* K* KAPOOR COMPUTER SCIENCES CORPORATION 


Subroutine OUTRES 
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SUBROUTINE OUTSUM(IA) 

purpose: 

OUTSUM PRODUCES A SUMMARY REPORT OF THE SOLVED-FOR 
COEFFICIENTS FOR A LINEAR POLYNOMIAL FITTING RUN AND FOR 
A NON-LINEAR LEAST SQUARES ESTIMATION RUN« 

method: 

THE VALUES OF TEMPERATUREr PRESSURE» AND DENSITY 
COEFFICIENTS RESIDE IN THE COMMON BLOCK COFCOM» FOR 
THE LINEAR POLYNOMIAL REPORT^ ALL THE COEFFICIENTS ARE 
PRINTED OUT. FOR THE NON-LINEAR REPORT » ONLY THE 
TEMPERATURE COEFFICIENTS AND THE FIRST PRESSURE COEFFICIENT 
ARE PRINTED OUT. 

ARGUMENT LIST: 

ARGUMENT TYPE I/O DESCRIPTION 

lA II FLAG FOR REPORT TYPE: 

“Or LINEAR POLYNOMIAL REPORT 
=lr NON-LINEAR REPORT 

CALLING subroutine: POLADJ 

SUBROUTINES CALLED: NONE 

COMMON BLOCK PARAMETERS USED: 

COMMON NAME PARAMETER I/O DESCRIPTION 

/COFCOM/ CT(3rA) I TEMPERATURE COEFFICIENTS 

CP(3r4) 1 PRESSURE COEFFICIENTS 

CD(3r4) I DENSITY COEFFICIENTS 

/FILCOM/ NPRINT I PRINTED OUTPUT LUN 

EXTERNAL DATA SETS USED: 

NAME LUN I/O OPERATIONS PERFORMED 

OUT NPRINT WRITE 

programmer/analyst: a. K. KAPOORr CSC MARCH 1981 


Subroutine OUTSUM 
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SUBROUTINE PARAMS(IERR) 
purpose; 

PARAMS IS THE PARAMETER CALCULATION EXECUTIVE. 

method; 

1. CHECK TO SEE IF AN ENTREE OR A LAP TRAJECTORY IS 

TO BE PROCESSED. CALL OUTPUT TO WRITE APPROPRIATE HEADERS. 

2. PRODUCE THE FIRST BLOCK OF TRAJECTORY POINTS BY CALLING 
ENTRD FOR AN ENTREE TRAJECTORY OR LAP FOR A LOCAL 
ATMOSPHERIC PROFILE. 

3. CALL MODELS TO RETURN TEMPERATUREr PRESSUREt DENSITY» AND 
WIND FOR EACH POINT IN THE BLOCK. 

A* CALL GASLAW TO CALCULATE PARAMETERS ASSOCIATED WITH THE 
IDEAL GAS LAW AND TO RECOMPUTE TEMPERATURE IN 
ACCORDANCE WITH THE GAS LAW. 

5. CALL PNCAL TO CALCULATE THE ONBOARD PRESSURE MEASUREMENT. 

6. CALL OUTPRF TO WRITE THE CALCULATED VALUES TO THE 
PRINTED REPORT AND TO THE USER OUTPUT FILE. 

7. CHECK TO SEE IF THE LAST POINT ON THE TRAJECTORY HAS 

BEEN PROCESSED. IF NOT, RETRIEVE THE NEXT BLOCK OF POINTS 
AND REPEAT STEPS 3 THROUGH A. 

8. RETURN. 

ARGUMENT LIST: 

ARGUMENT TYPE I/O DESCRIPTION 
lERR I 0 ERROR FLAG; 

0=N0 ERROR 

CALLING subroutine; LAIRS 

SUBROUTINES CALLED; LAP, ENTRD, MODELS, GASLAW, PNCAL, OUTPRF 


Subroutine PARAMS (1 of 2) 
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COMMON PLOCK PARAMETERS USED? 


COMMON NAME 

PARAMETER 

I/O 

DESCRIPTION 

/FILCOM/ 

NPRINT 

I 

PRINTED OUTPUT LUN 

/USECOM/ 

IPRFL 

I 

PROFILE TYPE 




=1» ENTREE 
=2» LAP 


IDBG 

I 

DEBUG FLAG 


Z1 

I 

LOWER LAP ALTITUDE 


Z2 

I 

UPPER LAP ALTITUDE 

/CALCOM/ 

TRJ(5» 10) 

0 

BLOCK OF TRAJECTORY POINTS Z» 
LATfLONGfGMTfLST 


UTRJ<5»10) 

0 

TRAJECTORY POINT UNCERTAINTIES 


PRM(7f 10) 

0 

BLOCK OF PARAMETERS CORRESPONDING 
TO EACH TRAJECTORY POINT: 

Tf Pr Df Uf USr UD 


UPRM<7» 10) 

0 

UNCERTAINTIES IN PARAMETERS 


MTMOLdO) 

0 

MEAN MOLECULAR WEIGHTS 


SCHT(IO) 

0 

PRESSURE SCALE HEIGHTS 


V(10) 

I 

SHUTTLE VELOCITIES (M/SEC) 


ALPHA(IO) 

I 

SHUTTLE ANGLES OF ATTACK 


RMC< 10) 

0 

CALCULATED MACH NUMBERS 


PNC(IO) 

0 

CALCULATED ONBOARD PRESSURES 


DEO< 10) 

0 

PERCENT DEVIATION FROM GAS LAW 

EXTERNAL DATA 

SETS used: 



NAME 

LUN 

I/O OPERATIONS PERFORMED 

OUT 

NPRNT 

WRITE 

(DEBUG ONLY) 


programmer: D* E» BOLAND* CSCr march 1981 


Subroutine PARAMS (2 of 2) 
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SUBROUTINE PHYSCS ( Z » T r P t D t lERR ) 

purpose: 

PHYSCS CALCULATES PRESSURE AND DENSITY FROM TEMPERATURE VIA 
PHYSICAL RELATIONSHIPS. 


method: 

1. CALL SIMP TO INTEGRATE THE INVERSE OF THE PRESSURE SCALE 
HEIGH-T TO THE CURRENT ALTITUDE ^SlMP IS A FORTRAN LIBRARY 
ROUTINE WHICH REQUIRES THE EXTERNAL ROUTINE PSCALE TO 
COMPUTE THE INVERSE SCALE HEIGHT). 

2. COMPUTE THE PRESSURE FROM THE INTEGRATED VALUE. 

3. COMPUTE THE DENSITY FROM THE PRESSURE. 

4 . RETURN. 


ARGUMENT LIST! 

ARGUMENT TYPE I/O 

Z R I 

T R I 

P R 0 

D R 0 

lERR I 0 


CALLING subroutines: poly 

SUBROUTINES CALLED: SIMP 


DESCRIPTION 
ALTITUDE (KM) 

TEMPERATURE AT Z (DEGREES K) 
PRESSURE AT 2 (N/M*»2) 
DENSITY AT Z (KG/M**3) 

ERROR FLAG 

NO ERROR 


COMMON BLOCK 

PARametFRS USED 

i 

COMMON NAME 

PARAMETER 

I/O 

/USECOM/ 

2REF 

X 


PREF 

I 


STEP 

1 


IDBG 

1 

/JRCOM/ 

CM 

I 


RC 

1 

/f ILCOM/ 

NPRINT 

I 

EXTERNAL DATA 

SETS used: 



NAME LUN 

OUT NPRINT 


DESCRIPTION 

REFERENCE ALTITUDE <KM) FOR 
INTEGRATION 
PRESSURE AT 2REF 
integration STEPSIZEi KM 
DEBUG flag: 

MEAN MOLECULAR MASS AT SEA LEVEL 
IN GM/MOLE 

UNIVERSAL GAS CONSTANT (MKS UNITS) 
PRINTED OUTPUT LUN 


I/O OPERATIONS PERFORMED 
WRITE (DEBUG ONLY) 


programmer: D. E. BOLANDr CSC 


Subroutine PHYSCS 
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SUBROUTINE PNCAL ( V » ALPHA » TPDW » EPHSEC t POS » RH f PN ) 


purpose; 

PNCAL CALCULATES THE ESTIMATED ONBOARD PRESSURE MEASUREMENT* 

method; 

1. CALCULATE THE SPEED OF SOUND AND MACH NUMBER. 

2. CALCULATE THE COEFFICIENT OF PRESSURE FOR THE SHUTTLE. 

3. CALCULATE THE DYNAMIC PRESSURE. 

4. CALCULATE THE ONBOARD PRESSURE MEASUREMENT. 

5. RETURN. 

ARGUMENT LIST; 


ARGUMENT 

TYPE 

I/O 

DESCRIPTION 

V 

R 

I 

SHUTTLE 

VELOCITY (M/SEC) 

ALPHA 

R 

1 

SHUTTLE 

ANGLE OF ATTACK (DEG) 

TPDW(7) 

R 

I 

CALCULATED ATMOSPHERIC PARAMETERS AT 




THIS TRAJECTORY POINT; T»P»D»U»V» 




USf WD 

(ONLY T»P»D -- THE FIRST THREE 




UORDS 

OF THE ARRAY -- ARE USED) 

EPHSEC 

R 

I 

SECONDS 

FROM EPOCH FOR THIS TRAJECTORY 




POINT 


P0S(4) 

R 

I 

SHUTTLE 

TRAJECTORY POINT; ALT, LAT, 




LOCAL 

SOLAR TIME, AND LONGITUDE 

RM 

R 

0 

CALCULATED MACH NUMBER 

PN 

R 

0 

CALCULATED ONBOARD PRESSURE MEASUREMENT 

CALLING SUBROUTINE 

; PARAMS 


SUBROUTINES CALLED 

; NONE 



COMMON BLOCK 

PARAMETERS 

used; 


COMMON NAME 

PARAMETER I/O 

DESCRIPTION 

/FILCOM/ 


NPRINT 

I 

PRINTED OUTPUT LUN 

/USECOM/ 


idbg 

I 

DEBUG FLAG 



HMSSEC 

I 

EPOCH (GMT SECONDS) 

/JRCOM/ 


RC 

I 

GAS CONSTANT (MKS UNITS) 



CM 

I 

MEAN MOLECULAR WEIGHT (G/MOLE) 

EXTERNAL DATA 

SETS 

used; 

NONE 


programmer; d 

. E. 

BOLAND 

» CSC» 

MAY 1981 


Subroutine PNCAL 



SUBROUTINE POLADJ(IERR) 

purpose: 

POLADJ IS THE POLYNOMIAL MODEL ADJUSTMENT EXECUTIVE. 

method: 

1. IF ADJUSTED POLYNOMIAL MODEL IS REOUESTEDt CALL METREF 
TO PRODUCE REFERENCE PROFILE. 

2. IF DEFAULT POLYNOMIAL MODEL IS REQUESTED. OR IF DEFAULT 
POINTS ARE TO BE ADDED TO REFERENCE METEOROLOGICAL 
PROFILE. CALL DFLINT. 

3. CALL POLFIT TO FIT THE POLYNOMIAL COEFFICIENTS TO THE 
REFERENCE PROFILE. 

A, RETURN. 


ARGUMENT list: 
ARGUMENT TYPE 

I/O 

DESCRIPTION 

I ERR I 

0 

ERROR flag: 



-O. NO ERROR 


CALLING subroutines: modadj 

SUBROUTINES CALLED: METREF. DFLINT. POLFIT 

COMMON BLOCK parameters used: 

COMMON NAME PARAMETER I/O 

/USECOM/ ILATL I 


IDSDL 1 


IDBG I 

IDFLPT X 


IDC X 

BOUND X 

RMERGE X 

/FILCOM/ NFRIMT X 

/XNTCOM/ IPTFLG 0 

external data sets used: 

NAME LUN I/O OPERATIONS PERFORMED 

OUT NPRZNT URITEtDEBUG ONLY) 

WRK NURK REWIND 

programmer: d.e. boland. csc 


Subroutine POLADJ 


description: 

INCLUDE LATITUDE VARIATIONS: 
«0 . NO 
-i. YES 

INCLUDE d/sd variations: 

>0 « NO 
-1. YES 
DEBUG FLAG 

INCLUDE DEFAULT POINTS: 

«0. NO 
«1. YES 
DC FLAG 

3 SEGMENT BOUNDARIES (KM) 
REGION FOR HERGING 
PRINTED OUTPUT LUN 
FLAG TO REINITIALIZE DFLINT 


3-71 



SUBROUTINE POLCAL ( Z » T » P f D » lERR > 


PURPOSE : 

POLCAL CALCULATES THE TEMPERATURE# PRESSURE# AND DENSITY 
FROM THE POLYNOMIAL COEFFICIENTS. 

method: 

1. COMPUTE TEMPERATURE# PRESSURE# AMD DENSITY FROM THE 
POLYNOMIAL COEFFICIENTS FOR THIS ALTITUDE (COMPUTE 
temperature only if the ‘DERIVE* OR *DC* OPTIONS HAVE 
BEEN SELECTED. 

2. RETURN. 

ARGUMENT list; 

ARGUMENT TYPE 

Z R 

T R 

P R 

0 R 

lERR I 


calling subroutines: poly# pscale 
subroutines called: none 


COMMON BLOCK 

PARAMETERS USED 

: 


COMMON NAME 

parameter 

I/O 

DESCRIPTION 

/USE COM/ 

I DTG 

I 

DEI'UC FLAG 


ICEG 

I 

ATMOSPHERIC SEGMENT NUMBER 


BOUNDO) 

1 

3 SEGMENT boundaries: 
DERAULT=90 KM# 6S KM# 25 KM 


RMERGE 

I 

REGION FOR MERGING ON EACH SIDE 
OF EACH SEGMENT BOUNDARY: 

default ■ 5 km 

/F ILCOM/ 

NPRINT 

I 

PRINTED OUTPUT LUN 

/COFCOM/ 

CT 

I 

A TEMFERATLfRE COEFFICIENTS FOR 
EACH SEGMENT 


CP 

1 

A PRESSURE COEFFICIENTS FOR 
EACH SEGMENT 

external DATA 

CD 

SETS USED! 

I 

A DENSITY COEFFICIENTS FOR 
EACH SEGMENT 

NAME 

LUN 


I/O OPERATIONS PERFORMED 

OUT 

NPRINT 


WRITE (DEBUG ONLY) 


programmer: d. e. boland# csc 


Subroutine POLCAL 


I/O description 
I altitude (KM) 

I temperature at Z (DEGREES K) 

0 PRESSURE AT Z (N/M**2) 

O DENSITY AT Z (KG/M*#3) 

0 ERROR FLAG 

»0» NO ERROR 
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SUBROUTINE POLERR < DZ f T r P » D » I SEG » UT f UP » I ERR ) 

purpose: 

POLERR ESTIMATES THE ERRORS (UNCERTAINTIES) ASSOCIATED WITH THE 
POLYNOMIAL CALCULATION OF TEMPERATUREf PRESSUREf AND DENSITY. 


method: 

1. SKIP THE NON-DC CALCULATION SINCE THE COUARINACE MATRIX 
IS UNAVAILABLE FROM THE LIBRARY ROUTINE LSOPOL . 

2. FOR A DC RUN» CALL DCERR TO CALCULATE THE SQUARE ROOT OF 
THE VARIANCE AT ALTITUDE DZ FOR EACH PARAMETER. 

3. RETURN. 


ARGUMENT LIST: 


ARGUMENT TYPE 

DZ R 

T R 

P R 

D R 

UT R 

UP R 

UD R 

lERR I 


I/O 

I 

I 

I 

I 

0 

D 

0 

0 


DESCRIPTION 

ALTITUDE (SEGMENT BOUNDARY - Z » KM) 

TEMPERATURE» DEGREES KELVIN 

PRESSURE^ N/M**2 

DENSITY » KG/M**3 

TEMPERATURE UNCERTAINTY 

PRESSURE UNCERTAINTY 

DENSITY UNCERTAINTY 

ERROR FLAG 


CALLING subroutine: POLY 


SUBROUTINES CALLED: LINERR» DCERR 


COMMON BLOCK 

PARAMETERS 

used: 



COMMON NAME 

PARAMETER I/O 

DESCRIPTION 


/COFCOM/ 

COVAT 

I 

TEMPERATURE COVARIANCE MATRIX 


COVAP 

I 

PRESSURE COVARIANCE 

MATRIX 


COVAD 

I 

DENSITY COVARIANCE 

MATRIX 


BOUND 

I 

POLYNOMIAL SEGMENT 

BOUNDARIES 

/USECOM/ 

IDC 

I 

DC FLAG 


EXTERNAL DATA 

SETS used: 

NONE 



programmer: d 

. E. BOLAND 

» CSC 

AUGUST 1981 



Subroutine POLERR 
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SUBROUTINE POLFIT ( ISEOf NHET » Z >T »P» D »UT r UPf UDf lERR ) 

PURPOSE • 

POLFIT CONTROLS THE FITTING OF POLYNOMIAL COEFFICIENTS, 

METHOD} 

1. IF THIS NOT A DC RUNf CALL LSQPOL 3 TIMES TO PERFORM LINEAR 
LEAST-SQUARES FITTING OF THE TEMPERATUREf PRESSUREf AND 
DENSITY COEFFICIENTS. CALL OUTRES TO OUTPUT THE RESIDUAL 
REPORTS. 

2. IF THIS IS A DC RUNt CALL NLSPOL TO ADJUST THE TEMPERATURE 
COEFFICIENTS VIA THE NONLINEAR LEAST SQUARES PROCESS. 

CALL OUTRES TO OUTPUT THE RESIDUAL REPORTS. 


3. RETURN 

ARGUMENT LIST 
ARGUMENT 

» 

t 

« 

TYPE 

I/O 

DESCRIPTION 

ISEG 

I 

I 

ATMOSPHERIC SEGMENT NUMBER 

NMET 

I 

I 

NUMBER OF POINTS IN THIS SEGMENT 

Z 

R 

I 

ALTITUDE ARRAY 

T 

R 

I 

TEMPERATURE ARRAY 

P 

R 

I 

PRESSURE ARRAY 

D 

R 

I 

DENSITY ARRAY 

WT 

R 

I 

TEMPERATURE HEIGHTS 

UP 

R 

I 

PRESSURE HEIGHTS 

UD 

R 

I 

DENSITY HEIGHTS 

lERR 

1 

0 

ERROR FLAG 


CALLING subroutine; POLADJ 


SUBROUTINES CALLED*. LSQPOLt NLSPOLf OUTRESf XALFAf UCMTRX 


COMMON BLOCK PARAMETERS USED} 


COMMON NAME PARAMETER I/O 

/COFCOM/ CT 0 

CP 0 

CD 0 

UT 0 

UP 0 

UD Q 

COUAT 0 


DESCRIPTION 

A TEMPERATURE COEFFICIENTS FOR 
EACH SEGMENT 

4 PRESSURE COEFFICIENTS FOR 
EACH SEGMENT 

4 DENSITY COEFFICIENTS FOR 
EACH SEGMENT 

TEMPERATURE STANDARD DEVIATION FOR 
EACH SEGMENT 

PRESSURE STANDARD DEVIATION FOR 
EACH SEGMENT 

DENSITY STANDARD DEVIATION FOR 
EACH SEGMENT 

TEMPERATURE COVARIANCE MATRIX FOR 
EACH SEGMENT 


Subroutine POLFIT (1 of 2) 
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/USECOM/ 

IDBO 

I 

DEBUG FLAG 


I DC 

I 

DC FLAG 


ITDAT 

1 

DC TEMPERATURE SELECTION FLAG 


IPDAT 

I 

DC PRESSURE SELECTION FLAG 


IDDAT 

I 

DC DENSITY SELECTION FLAG 


KWRITE 

1 

OUTPUT SUPPRESSION FLAG 


BOUND 

I 

SEGMENT BOUNDARIES 


ZREF 

I 

REFERENCE ALTITUDE FOR INTEGRATION 


PREF 

I 

PRESSURE AT ZREF 

/FILCOM/ 

NPRINT 

I 

PRINTED OUTPUT LUN 

/ JRCOM/ 

TJREF 

I 

TEMPERATURE AT TOP BOUNDARY 


RC 

I 

GAS CONSTANT (MKS UNITS) 


RMAV 

I 

MEAN MOLECULAR WEIGHT (CGS UNITS) 


DJREF 

I 

DENSITY AT TOP BOUNDARY 

/DCCOM/ 

MLDATA 

I 

DATA SELECTION FLAGS 


KSSEG 

I 

CONTINUITY FLAG 


ISGH 

I 

SEGMENTATION FLAG 


IDCSEG 

I 

DC FLAG FOR EACH SEGMENT 


ALFA 

I/O 

INTEGRANDS FOR P AND D 

EXTERNAL DATA 

SETS used; 



NAME 

LUN 


I/O OPERATIONS PERFORMED 

OUT 

NPRINT 


WRITE (DEBUG ONLY) 


PROGRAMhERJ D. £♦ BOLAND » CSC 


Subroutine POLFIT (2 of 2) 
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SUBROUTINE POLYCPOSt TPDU»UTPDW» lERR) 

purpose: 

POLY CONTROLS THE CALCULATION OF ATMOSPHERIC PARAMETERS 
VIA THE POLYNOMIAL MODEL* 


method: 

1. DETERMINE WHICH POLYNOMIAL SEGMENT THIS POINT IS IN* 

2* CALL POLCAL OR PHYSCS TO CALCULATE TEMPERATURE » PRESSUREf 
AND DENSITY* 

3* CALL DFLINT TO READ VALUES OF TRANSLATIONAL COEFFICIENTS 
FROM WORKING FILE* CALL TRNLAT AND TRNOSD TO TRANSLATE 
PARAMETERS IN LATITUDE AND LOCAL SOLAR TIME* 

4* CALCULATE UNCERTAINTIES 
5* RETURN* 


ARGUMENT LIST: 

ARGUMENT TYPE 
POS R 


TPDW R 


UTPDW R 

lERR I 


I/O DESCRIPTION 
I TRAJECTORY POINT: 

POS(l)=ALTITUDE 
POS(2)=LATITUDE 
P0S<3)=L0CAL SOLAR TIME 
POS(4)=LONGITUDE 

0 CALCULATED ATMOSPHERIC PARAMETERS: 

TPDW( I ):*TEMPERATURE 
TPDW<2)*PRESSURE 
TPDW<3)=DENSITY 
TPDU(4)»E>W WIND U 
TPDW(5)*N-S WIND V 
TPDW(6):bUIND SPEED WS 
TPDW(7)«WIND DIRECTION WD 
0 UNCERTAINTIES IN PARAMETERS (7) 

0 ERROR flag: 

-Of NO ERROR 


CALLING subroutine: MODELS 


SUBROUTINES CALLED: POLCALr PHYSCS* DFLINT* TRNLAT* TRNDSD* POLERR 


Subroutine POLY (1 of 2) 
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COMMON BLOCK 

parameters used: 



COMMON NAME 

parameter I/O 

DESCRIPTION 

/USECOM/ 

ILATL 

I 

INCLUDE LATITUDE VARIATIONS! 
= Of NO 
= lf YES 


IDSDL 

I 

INCLUDE D/SD VARIATIONS! 
= 0» NO 
=1 »YES 


IDRV 

I 

DERIVE P AND D FROM T! 
-Op no 
>=1p yes 


ZDBG 

I 

DEBUG FLAG 


iseg 

0 

ATMOSPHERIC SEGMENT NUMBER 


B0UND(3) 

I 

3 SEGMENT BOUNDARIES! 
DEFAULT=90 KMp 65 KMp 25 KM 


RMER6E 

I 

REGION FOR MERGING ON EACH SIDE 
OF EACH SEGMENT BOUNDARY! 
DEFAULT = 5 KM 


REFLAT 

I 

MODEL REFERENCE LATITUDE 


REFTAU 

I 

MODEL REFERENCE LOCAL SOLAR TIME 

/FILCOM/ 

NPRINT 

I 

PRINTED OUTPUT LUN 

/COFCOM/ 

CP 

I 

PRESSURE COEFFICIENTS ARRAY 


RMST(3) 

1 

TEMPERATURE RMS FOR EACH SEGMENT 


RMSP(3) 

I 

PRESSURE RMS FOR EACH SEGMENT 

EXTERNAL DATA 

RHSD(3> 
SETS used: 

I 

DENSITY RMS FOR EACH SEGMENT 

NAME 

LUN 


I/O OPERATIONS PERFORMED 

OUT 

programmer: d 

NPRINT 

. E* BOLANDf CSC 

WRITE (DEBUG ONLY) 


Subroutine POLY (2 of 2) 
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SUBROUTINE PPAR ( RINT » DPDS > 
purpose: 

PPAR COMPUTES PARTIAL DERIVATIVES OF PRESSURE WITH RESPECT TO 
SOLVE-FOR PARAMETERS 


method: 

USING PREVIOUSLY COMPUTED INTEGRALS AND CURRENT SOLVE-FOR 
PARAMETERS^ COMPUTES PARTIAL DERIVATIVES OF LOG PRESSURE 
WITH RESPECT TO TEMPERATURE COEFFIENTS AND LOG BOUNDARY 
PRESSURE 

CALLING subroutine: NLSPOL 
SUBROUTINES CALLED: NONE 


COMMON BLOCK PARAMETERS USED: 
COMMON NAME PARAMETER I/O 


/SOLAND/ IMAX I 

KSOLA4 I 

NPOL I 

A I 


DESCRIPTION 

MAX, NUMBER OF INTEGRALS 

=0 DOES NOT SOLVE FOR BOUNDARY TEMP. 

=1 SOLVES FOR BOUNDARY TEMP, 

ORDER OF TEMP. POLYNOMIAL 
SOLVE-FOR PARAMETER ARRAY 


EXTERNAL DATA SETS USED: NONE 


programmer; t, lee» esc 


MARCH 1981 


Subroutine PPAR 
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SUBROUTINE PSC ALE < Z » HP I ) 


PURPOSE : 

PSCALE COMPUTES THE INVERSE OF THE PRESSURE SCALE HEIGHT. 

method: 

1. COMPUTE THE ACCELERATION OF GRAVITY AT THIS ALTITUDE. 

2. COMPUTE THE INVERSE OF THE SCALE HEIGHT. 

3* RETURN. 

ARGUMENT list: 

ARGUMENT TYPE 
2 R 

HPI R 


I/O DESCRIPTION 
I ALTITUDE (KM) 

0 INVERSE PRESSURE SCALE HEIGHT (KM*»-1) 


CALLING subroutines: SIMP 

subroutines called: polcal 


COMMON BLOCK PARAMETERS USED: 

common name parameter I/O DESCRIPTION 


/ JRCOM/ 

GLO 

I 

SEA LEVEL value OF G <M/SEC**2> 


RAV 

I 

AVERAGE EARTH RADIUS (KM) 


CM 

I 

MEAN molecular MASS AT SEA 




LEVEL (GM/MOLE) 


TZ 

1 

TEMPERATURE AT ALTITUDE Z (KELVIN) 


RC 

X 

UNIVERSAL GAS CONSTANT (MKS UNITS) 

external data 

SETS used: 

NONE 


programmer: d. 

E. BOLAND 

» CSC 



Subroutine PSCALE 
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SUBROUTINE READIN<NFINf lERR) 

purpose: 

READIN READS THE KEYWORD CARD DECK AND» IF REQUIREDf THE 
ENTREE FILE HEADER^ AND STORES THE INPUT PARAMETERS IN COMMON, 

method: 

1, WRITE THE TITLE FOR THE KEYWORD CARD IMAGES, 

2, READ A CARD AND WRITE ITS IMAGE. 

3, BRANCH TO THE APPROPRIATE SECTION OF CODE TO SET THE COMMON 
BLOCK PARAMETERS INPUT ON THIS CARD, IF THIS AN ENTREE 
PROFILE RUN» READ THE HEADER ON THE ENTREE FILE ALSO. 

4, REPEAT STEPS 2 AND 3 UNTIL THE 'END' CARD IS READ, 

5, RETURN, 


ARGUMENT LIST 




ARGUMENT 

TYPE 

I/O 

DESCRIPTION 

NFIN 

I 

0 

FLAG TO END STACKED DECK PROCESSING 
=1» LAST DECK 

lERR 

I 

0 

ERROR flag: 

NO ERROR 


CALLING subroutine: LAIRS 
SUBROUTINES CALLED*. HMSCON 


COMMON BLOCK PARAMETERS USED: 
COMMON NAME PARAMETER I/O 


/USECOM/ IPRFL 0 

MODL 0 

IWL 0 

ILATL D 

IDSDL 0 

MO DU 0 

IWU 0 

ILATU 0 

IDSDU 0 

I DRV 0 


INTP 0 

IDBG 0 

INTPB 0 

IDNUM 0 

IGAS 0 

NSAMP 0 

ITER 0 

IWT 0 

IDSDA 0 

ILATA 0 

IDC 0 

ITBAT 0 

IPDAT 0 


DESCRIPTION 
TYPE OF PROFILE 

TYPE OF LOWER ATMOSPHERIC MODEL 
TYPE OF L.A. WIND MODEL 
LATITUDE VARIATIONS IN L.A. 
DIURNAL/SEMIDIURNAL VARIATIONS^ L.A, 
TYPE OF UPPER ATMOSPHERIC MODEL 
TYPE OF U.A. WIND MODEL 
LATITUDE VARIATIONS IN U.A. 
DIURNAL/SEMIDIURNAL VARIATIONS^ U.A, 
DERIVE PRESSURE AND DENSITY FROM 
TEMPERATURE VIA PHYSICAL LAWS 
(OVERRIDES ALL OTHER SELECTIONS) 
INTERPOLATOR TYPE 
PRINT DEBUG 

BIVARIATE INTERPOLATOR TYPE 
ENTREE-LAP ID NUMBER 
GAS LAW CALCULATION OVERRIDE 
ENTREE FILE SAMPLING INTERVAL 
(EACH NSAMPTH POINT READ) 

MAXIMUM NUMBER OF DC ITERATIONS 
WEIGHTING FLAG 
D/SD COEFFICIENT ADJUSTMENT 
LATITUDE GRADIENT ADJUSTMENT 
DC FLAG 

SOURCE OF T OBSERVATIONS FOR DC 
SOURCE OF P OBSERVATIONS FOR DC 


Subroutine READIN (1 of 2) 
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I 



IDDAT 

0 

SOURCE OF D OBSERVATIONS FOR DC 


IDMY 

0 

DUMMY 


NFIL 

0 

NUMBER OF MET FILES 


KWRITE 

0 

OUTPUT SUPPRESSION FLAG 


YMD 

0 

YEAR» MONTHf DAY OF LAP OR ENTREE 




TRAJECTORY 


HNS 

0 

HOURS* MINUTES* SECONDS OF FIRST 




POINT ON ENTREE TRAJECTORY 


THETA 

0 

LATITUDE OF LAP 


TAU 

0 

LONGITUDE OF LAP 


Z1 

0 

LOWER LAP ALTITUDE 


Z2 

0 

UPPER LAP ALTITUDE 


DZ 

0 

LAP ALTITUDE INTERVAL 


BOUND 

0 

BOUNDARY BETWEEN SEGMENTS 


HMSSEC 

0 

ENTREE FILE HMS START TIME IN 




DECIMAL SECONDS 


ZREF 

0 

REFERENCE ALTITUDE <KM) 


PREF 

0 

PRESSURE AT ZREF <N/M2) 


STEP 

0 

STEPSIZE 


CONV 

0 

CONVERGENCE CRITERION FOR DC 


RNERGE 

0 

REGION FOR SEGMENT MERGING 


REFLAT 

0 

REFERENCE LATITUDE 


REFTAU 

0 

REFERENCE LOCAL SOLAR TIME 


RMGNT 

0 

GMT OF MET FILES 

/FILCOh/ 

NPRINT 

I 

PRINTED OUTPUT LUN 


NREAD 

I 

CARD INPUT LUN 


NENT 

I 

ENTREE FILE LUN 

/DCCOM/ 

NLDATA 

0 

DC DATA SELECTION FLAGS 


K8SEG 

0 

DC CONTINUITY FLAGS 


I8GM 

0 

DC SEGMENTATION FLAG 


ISE6DC 

0 

DC XNCLUSON FLAG FOR SEGMENT 


MTDCT 

0 

WEIGHTING FOR TEMPERATURE 


WTDCP 

0 

WEIGHTING FOR PRESSURE 


WTDCD 

0 

WEIGHTING FOR DENSITY 

EXTERNAL DATA 

SETS USED 



NAHE 

LUN 

I/O OPERATIONS PERFORMED 

OUT 

NPRINT 

WRITE 


CARDS 

NREAD 

READ 


ENT 

NENT 

READ 


prograhner; d* 

E* BOLAND r 

CSC^ 

MARCH 1981 


Subroutine READIN (2 of 2) 
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FUNCTION RLSTtLONGf SECHHSfEPHSEC) 


purpose: 

rest converts longitude and GrtTdN SECONDS) INTO LOCAL SOLAR 
TIHE(SOLAR HOUR ANGLE). 

method: 

1. CONVERT LONGITUDE<DEGREES) TO SECONDS OF TIME EAST OF 
GREENUICH. 

2. ADD GMT IN SECONDS TO TIME EAST OF GREENUICH TO PRODUCE 
LOCAL SOLAR TIME IN SECONDS. 

3. CONVERT LST IN SECONDS TO LST IN HOURS. 


ARGUMENT LIST: 


ARGUMENT TYPE 
LONG R 

SECHMS R 

EPHSEC R 

REST R 


I/O DESCRIPTION 

I LONGITUDE EAST OF GREENUICH (DEGREES) 

I EPOCH HMS CONVERTED TO SECONDS 

I SECONDS FORM EPOCH (DECIMAL SEC. OF TIME) 

0 LOCAL SOLAR TIME (DECIMAL HOURS) 


CALLING SUBROUTINES: PARAM8 


SUBROUTINES CALLED: NONE 

COMMON BLOCK PARAMETERS USED: NONE 

EXTERNAL DATA SETS USED: NONE 

programmer: r.a. kusesm» esc 


Function RLST 



SUPROUT I NE ROOTS ( A » NA » CROOTS r I RL) 


PURPOSE t 

THIS SUBROUTINE COHPUTES ANY DESIRED COMPLEX ROOTS (ZEROES) OF A 
GIOEN POLYNOMIAL y IF APPROXIMATIONS TO THE DESIRED ROOTS ARE KNOWN* 

METHOD ? 

NEWTON'S METHOD? AS DESCRIBED IN THE REFERENCE* 

REFERENCE ? 

HENRI'CI? P,? 'ELEMENTS OF NUMERICAL ANALYSIS'? NEW YORK. WILEY. 
.1.965. PAGE 8^* 

CALI..;i;NG SEQUENCE 1 

CALL ROOTS ■ A ? N A . CROOTS . I RL) 

WHERE : 

A IS AN INPUT DIMENSION NA ARRAY OF REAL ( NON I MMAG I N ARY ) POLY- 

NOMIAL COEFFICIENTS. STARTING WITH THE LOWEST POWER. 

NA IS THE NUMBER OF IMPUT COEFFICIENTS (ONE MORE THAN THE 
D E G R E E 0 F T H E P 0 1... Y N 0 M I A 

CROOTS iS A CuMPLEX:i<i6 ARRAY OF' D I MENS TON NA CONTAINING THE 

STARTING APPROXIMATIONS TO THE DESIRED ROOTS ON INPUT AND THE 
SOLUTIONS ON OUTPUT... 

IRL IS THE NUMBER OF ROOTS (ZEROES) DESIRED TO BE SOLOED FOR. 

COMMON BLOCK PARAMETERS? 

NONE 

SUBROUTINES CALLED S 
NONE 

CALLED BY? 

LOUAL r 

programmer: 

•J. P, MOl. INEAUX? COMPUTER SCIENCES CORPORATION. 


Subroutine ROOTS 
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SUBROUTINE SETREF(IERR) 


purpose: 

SETREF SETS THE REFERENCE VALUES AT THE TOP SEGMENT BOUNDARY 
AND AT ZREF FOR THE JR MODEL* THE DC PROCESS* AND THE 
DERIVATION OF P AND D FROM T BY INTEGRATION* 

method: 

1. WRITE REPORT HEADER* 

2. IF METEOROLOGICAL FILES ARE BEING USED* READ HEADERS AND 
CALL METINT TO INTERPOLATE THE NEEDED VALUES* IF DEFAULT 
MODELS ARE BEING USED* CALL DFLINT TO INTERPOLATE NEEDED 
VALUES FROM THE DEFAULT FILES* 

3. REPORT THESE VALUES AND HEADER INFORMATION* 

4* RETURN* 

ARGUMENT LIST: 

ARROUMENT TYPE I/O DESCRIPTION 
lERR I 0 ERROR FLAG 

CALLING SUBROUTINE: LAIRS 

SUBROUTINES CALLED: METINT* DFLINT 

COMMON BLOCK PARAMETERS USED: 


COMMON NAME 

PARAMETER 

I/O 

DESCRIPTION 

/USECOM/ 

MODL 

1 

TYPE OF LOWER ATMOSPHERIC MODEL 


I DRV 

I 

FLAG FOR P AND D DERIVATION 


NFIL 

I 

NO* OF METEOROLOGICAL FILES 


BOUND<3) 

I 

ATMOSPHERIC SEGMENT BOUNDARIES 


ZREF 

I 

INTEGRATION REFERENCE HEIGHT 


PREF 

0 

PRESSURE AT ZREF 


REFLAT 

I 

MODEL REFERENCE LATITUDE 


REFTAU 

I 

MODEL REFERENCE LOCAL SOLAR TIME 

/JRCOM/ 

ZBND 

0 

TOP SEGMENT BOUNDARY 


TBND 

0 

TEMPERATURE AT ZBND 


DBND 

0 

DENSITY AT ZBND 

/FILCOM/ 

NPRINT 

I 

PRINTED OUTPUT LUN 


NURK 

I 

DEFAULT WORKING FILE LUN 


MWF<3) 

I 

METEOROLOGICAL FILES LUNS 

/INTCOM/ 

12 

0 

INTERPOLATION INITIALIZATION FLAG 

/METCOM/ 

M12 

0 

INTERPOLATION INITIALIZATION FLAG 

EXTERNAL DATA 

SETS used: 



NAME 

LUN 


I/O OPERATIONS PERFORMED 

OUT 

NPRINT 


WRITE 

WRK 

NURK 


READ* REWIND 

MET(3) 

MWF<3) 


READ* REWIND 

programmer: d» 

E* BOLAND* CSC* 

MARCH 1961 


Subroutine SETREF 



SUBROUTINE SOLVE ( KSOL > KDPMET » DELTAS f VCHAT i I ERR ) 


purpose: 

SOLVE SOLVES THE NORMAL EQUATION BY INVERTING THE NORMAL MATRIX* 

method: 

USE THE SUBROUTINE SYMINV TO INVERT THE NORMAL MATRIX* 

ARGUMENT LIST: 

ARGUMENT TYPE I/O DESCRIPTION 

KSOL I I =1 SOLVES FOR ALL FIVE PARAMETERS 

-2 SOLVES FOR THREE TEMP* COEFFS 
=3 SOLVES FOR ALL TEMP* COEFFS 
=4 SOLVES FOR THREE TEMP* AND ONE PRESS* 

COEFFS 

KDPMET I I =0 NO PRES* OR DEN* DATA USED 

=1 PRES* OR DEN* OR BOTH TYPES OF DATA USED 
DELTAS R 0 CORRECTION SOLVE-FOR VECTOR 

VCMAT R 0 VARIANCE-COVARIANCE MATRIX 

lERR I 0 ERROR FLAG 

CALLING subroutines: NLSPOL 

SUBROUTINES CALLED: SYMINV 

COMMON BLOCK PARAMETERS USED** 

COMMON NAME PARAMETER I/O DESCRIPTION 

/USECOM/ KPRNT I UNIT NUMBER FOR PRINTING 

/SOLAND/ KS0LA4 I =0 DOES NOT SOLVE FOR TEMP* BOUND* 

VALUE 

=1 SOLVES FOR TEMP* BOUND*VALUE 
KS0LA5 I =0 DOES NOT SOLVE FOR PRES* BOUND* 

VALUE 

=1 SOLVES FOR PRES* BOUND* VALUE 
NPOL I ORDER OF POLYNOMIALS FOR TEMP* 

MAXSOL I MAX* NUMBER OF SOLVE-FOR PARAMETERS 

NSOLV I NUMBER OF SOLVE-FOR PARAMETERS 

KEDTOT I TOTAL NUMBER OF MET* DATA INCLUDED 

K0NVR6 0 CONVERGENCE FLAG 

=1 CONVERGED 
=2 CONVERGING 
=3 DIVERGING 

=4 converged: linear FITTING 
EPSLON I CONVERGENCE CRITERIA 

SRSUM I TOTAL SUM OF SQUARED UEIGHTED RESIDUALS 

CURRMS 0 CURRENT WEIGHTED RMS 

OLDRMS 0 PREVIOUS UEIGHTED RMS 

RMSP 0 PREDICTED WEIGHTED RMS 

A I/O SOLVE-FOR PARAMETER ARRAY 

G I RIGHT-HAND SIDE OF NORMAL EQUATION 

SH I LINEAR ARRAY OF NORMAL MATRIX 

EXTERNAL DATA SETS USED: NONE 

programmer: T* LEE» esc march 19S1 


Subroutine SOLVE 
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SUBROUTINE SYMINV < SUMl , NDIM » NLIM i DELTA f lERR > 

purpose: 

SYMINV INVERTS A SYMMETRIC COEFFICIENT MATRIXf REQUIRING ONLY 
THE UPPER TRIANGULAR PORTION OF THE MATRIX* 

method: 

1. RECURSIVELY FIND THE INVERSE OF N X N MATRIX KNOWING THE 
INVERSE OF (N-1) X (N-1) MATRIX UNTIL THE INVERSE OF AN 
NLIM X NLIM SQUARE PARTITION MATRIX IS FOUND (SCHUR 
IDENTITY METHOD OF PARTITIONING)* 

2* RETURN* 

ARGUMENT LIST: 

ARGUMENT TYPE I/O DESCRIPTION 

SUMl R I/O i: UPPER TRIANGULAR FORM OF SYMMETRIC 

MATRIX STORED BY ROW 
01 UPPER TRIANGULAR FORM OF SYMMETRIC 
MATRIX WITH SUBMATRIX INVERTED 
NDIM I I DIMENSION OF MATRIX 'SUM!' 

NLIM I I DIMENSION OF SUBMATRIX TO BE INVERTED 

DELTA R I WORKING ARRAY 

lERR I 0 ERROR FLAG 

CALLING subroutine: SOLVE 

SUBROUTINES CALLED: COVUP 

COMMON BLOCK PARAMETERS USED: NONE 

EXTERNAL DATA SETS USED: NONE 

programmer: DOLORES P. M0RGAN» GSFC» NASA (FOR GTDS) 

MODIFIED FOR USE IN LAIRS BY D* E, BOLAND» CSC 


Subroutine SYMINV 
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SUBROUTINE TPAR ( ZR » DTDS ) 

purpose: 

TPAR COMPUTES PARTIAL DERIVATIVES OF TEMPERATURE WITH RESPECT TO 
SOLVE-FOR PARAMETERS 

method: 

USING THE POLYNOMIAL EXPRESSION OF TEMPERATURE» COMPUTE PARTIAL 
DERIVATIVES WITH RESPECT TO THE TEMPERATURE COEFFIENTS 

ARGUMENT LIST: 

ARGUMENT TYPE I/O DESCRIPTION 

ZR R I ALTITUDE WITH RESPECT TO THE REFERENCE 

ALTITUDE 

DTDS R 0 COMPUTED PARTIAL DERIVATIVE ARRAY 

CALLING subroutine: NLSPOL 

SUBROUTINES CALLED: NONE 

COMMON BLOCK PARAMETERS USED; 

COMMON NAME PARAMETER I/O DESCRIPTION 

/SOLAND/ NPOL I ORDER OF TEMP, POLYNOMIAL 

KS0LA4 I =0 DOES NOT SOLVE FOR BOUNDARY TEMP, 
=1 SOLVES FOR BOUNDARY TEMP. 

EXTERNAL DATA SETS USED: NONE 
programmer: T. LEEt CSC march 1981 


Subroutine TPAR 
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SUBROUTINE TRNADJ(IERR) 

purpose: 

TR/MDJ CONTROLS THE ADJUSTMENT OF LATITUDE GRADIENTS AND 

diurnal/semidiurnal coefficients, 
method: 

1. RETRIEVE DEFAULT COEFFICIENTS FROM UORK AREA /SCRATCH/. 

2. INTERPOLATE It Ft AND D FROM THE METEOROLOGICAL FILES 

AT EACH ALTITUDE REPRESENTED BY THE DEFAULT COEFFICIENTS. 

3. CALCULATE THE NEW LATITUDE GRADIENTS 

4. CALCULATE THE NEW D I URNAL/SEM I D I URNAL COEFFICIENTS. 

3. WRITE the new COEFFICIENTS TO THE WORK AREA SCRATCH. 

6. RETURN. 

ARGUMENT list: 

ARGUMENT TYPE I/O DESCRIPTION 
lERR I 0 ERROR FLAG 

CALLING subroutine: UFILE 

subroutines called: METINI »RLSTf LATAOJrDSDADJ 

COMMON BLOCK parameters USED.* 

COMMON NAME PARAMETER I/O 
/USECOM/ NFIL 1 

RGMT<3) I 

ILATA I 

IDSDA I 

/FILCDM/ NPRINT 1 

MURK I 

MWF<3> X 

/SCRATCH/ TAD<23»30) I 

/INTCOM/ 12 X 

/METCOM/ MI2 I 

EXTERNAL DATA SETS USEOI 

NAME LUN I/O OPERATIONS PERFORMED 

OUT NPRINT WRITE (DEBUG ONLY) 

MUF MWF(3) REWIND 

programmer: D. E. BOLANDp CSC» march 1981 


Subroutine TRNADJ 


DESCRIPTION 

NO. OF meteorological FILES 
TIMES OF MET, FILES <GMT) 

FLAG FOR LAT. GRAD. ADJUSTMENT 

FLAG FOR D/SD COEF . ADJUSTMENT 

PRINTED OUTPUT LUN 

WORKING file LUN 

MET. files LUNS 

CONTENTS OF WORKING FILE 

(WORK AREA) 

FLAG TO reinitialize DFLINT 
FLAG TO REINITIALIZE METINT 
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FUNCTION TRNDSD(F»DAF»DPF, SDAF,SriPFTTAU»R.EFTAU) 
purpose: 

TRNDSD TRANSLATES THE PARAMETER F FROM A REFERENCE SOLAR 
TIME TO LOCAL SOLAR TIME* 


method: 

EXTRAPOLATE F IN TIME UITH THE DIURNAL AND SEMIDIURNAL 
COEFFICIENTS OF F (PHASE IS LOCAL SOLAR TIME OF MAXIMUM)* 


ARGUMENT LIST 


ARGUMENT TYPE 
F R 

DAF R 

DPF R 

SDAF R 

SDPF R 

TAU R 

REFTAU R 

TRNDSD R 


I/O 

I 

I 

I 

I 

I 

I 

I 

0 


DESCRIPT ION 

PARAMETER TO BE TRANSLATED 
DIURNAL AMPLITUDE OF F 
DIURNAL PHASE OF F 
SEMIDIURNAL AMPLITUDE OF F 
SEMIDIURNAL PHASE OF F 
LOCAL SOLAR TIME 
REFERENCE SOLAR TIME 
TRANSLATED MALUE OF F 


CALLING subroutines: DFL I NT y MET I NT r LOWATM 


SUBROUTINES CALLED: NONE 

COMMON BLOCK PARAMETERS USED: 

COMMON NAME PARAMETER I/O DESCRIPTION 
/CONCOM/ TWOPI I TWO TIMES PI 

EXTERNAL DATA SETS USED: 

programmer: d* e* boland» esc 


Function TRNDSD 
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FUNCTION TRNLATCF, DF, RL AT » REFLAT) 

purpose; 

TRNLAT TRANSLATES THE PARAMETER F FROM A REFERENCE LATITUDE 
TO ANOTHER LATITUDE 


method; 

EXTRAPOLATE F IN LATITUDE WITH THE LATITUDE GRADIENT OF F 


ARGUMENT LIST 


ARGUMENT TYPE 

F R 

DF R 

RLAT R 

RE FLAT R 

TRNLAT R 


I/O DESCRIPTION 

I PARAMETER TO BE TRANSLATED 

I LATITUDE GRADIENT OF F 

I LATITUDE TO WHICH F IS TO BE TRANSLATED 

I REFERENCE LATITUDE 

0 TRANSLATED VALUE OF F 


CALLING subroutines; DFL I NT r MET INT ? LOWATM 


SUBROUTINES CALLED; NONE 


COMMON BLOCK PARAMETERS USED; NONE 


EXTERNAL DATA SETS USED; NONE 

programmer; d* e* boland» esc 


Function TRNLAT 
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SUBROUTINE UCNTRX ( ISEGHf BOUND »PBOUND t ALFA » COUAT » lERR ) 
PURPOSE} 

VCHTRX COMPUTES THE CONTRIBUTION TO THE 5X5 VARIANCE- 
COVARIANCE MATRIX FROM THE FIXED BOUNDARY TEMPERATURE 
AND PRESSURE. 

method: 

1. COMPUTE A 2 X 2 VARIANCE-COVARIANCE MATRIX FOR FIXED 
BOUNDARY TEMPERATURE AND PRESSURE USING THE VARIANCE- 
COVARIANCE MATRIX OF THE PREVIOUS SEGMENT* 

2* CONSTRUCT A DIRECT SUM OF THE 3X3 VARIANCE-COVARIANCE 
MATRIX OF THE CURRENT SEGMENT AND THE 2X2 MATRIX 
COMPUTED IN STEP 1. 


ARGUMENT LIST 
ARGUMENT 

• 

TYPE 

I/O 

DESCRIPTION 

1SE6M 

I 

I 

SEGMENT NUMBER 

BOUND 

R 

I 

SEGMENT BOUNDARIES 

PBOUND 

R 

I 

PRESSURE AT UPPER BOUNDARY 

ALFA 

R 

I 

INTEGRALS NEEDED FOR COMPUTING 

COVAT 

R 

I/O 

PARTIAL DERIVATIVES 
VARIANCE-COVARIANCE MATRIX 

lERR 

I 

0 

ERROR FLAG 


CALLING subroutines: POLFIT 


SUBROUTINES CALLED: NONE 

COMMON BLOCK PARAMETERS USED} NONE 

EXTERNAL DATA SETS USED: NONE 

programmer: t* lee» csc» September i9si 


Subroutine VCMTRX 
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SUBROUTINE UCNVRT <Ut WSr WD»UU»UVf UUSrUULD 

purpose: 

WCNURT CONVERTS UIND COMPONENTS U AND V TO WIND SPEED US AND 
WIND DIRECTION UD* IT ALSO CONVERTS THE UNCERTAINTIES* 


method: 

1* PERFORM COMPONENTS CONVERSION. 

2. PERFORM UNCERTAINTIES CONVERSION* 


3. RETURN 




ARGUMENT LIST 

♦ 



ARGUMENT 

TYPE 

I/O 

DESCRIPTION 

U 

R 

I 

EAST-WEST UIND COMPONENT (M/SEC) 

V 

R 

I 

NORTH-SOUTH UIND COMPONENT <M/8EC) 

US 

R 

0 

UIND SPEED (M/SEC) 

UD 

R 

0 

UIND DIRECTION (DEGREES) 

UU 

R 

I 

E-W WIND UNCERTAINTY <M/SEC) 

UV 

R 

I 

N-S WIND UNCERTAINTY (M/SEC) 

UUS 

R 

0 

UIND SPEED UNCERTAINTY (M/SEC) 

UUD 

R 

0 

WIND DIRECTION UNCERTAINTY (M/SEC) 


CALLING subroutines: MODELS 
SUBROUTINES CALLED: NONE 


COMMON BLOCK PARAMETERS USED** 
COMMON NAME PARAMETER I/O 
/CONCOM/ RTD I 

DTR I 


DESCRIPTION 

RADIANS TO DEGREES CONVERSION 
FACTOR 

DEGREES TO RADIANS CONVERSION 
FACTOR 


EXTERNAL DATA SETS USED: NONE 

prdorammer: d* e* boland* esc 


Subroutine WCNVRT 



SUBROUTINE UFILE<IERR) 


PURPOSE : 

ufile builds the working file from the permanent file. 


method: 

1. LOCATE THE PROPER MONTH ON THE PERMANENT FILE. 

2. READ THE PROFILE FOR THAT MONTH. 

3. REVERSE THE ALTITUDE ORDERING OF THE PROFILE AND WRITE IT 
TO THE UORKiNli TILE. 

4. ADJUST LATITUDE GRADIENTS AND 0 lURNAL/SEMI DIURNAL 
COEFFICIENTS IF REQUIRED. 

5. RETURN. 

ARGUMENT list: 

ARGUMENT TYPE I/O DESCRIPTION 
lERR 1 0 ERROR FLAG: 

aO. NO ERROR 


CALLING subroutine: lairs 


SUBROUTINES CALLED: TRNADJ 
COMMON BLOCK PARAMETERS USED: 


COMMON NAME 

PARAMETER 

I/O DESCRIPTION 

/USECOM/ 

IDBG 

I 

DEBUG FLAG 


RYMD 

1 

year»month»day of entree 

TRAJECTORY OR LAP 


ILATA 

X 

FLAG FOR LAT, GRAD. ADJUSTMENT 


IDSDA 

I 

FLAG r uR D/SD COEF . ADJUSTMENT 

/F ILCOM/ 

NPRINT 

1 

PRINTED OUTPUT LUN 


NF RM 

1 

permanent file LUN 


NURK 

I 

UORMNG file LUN 

/SCRATCH/ 

TAD(23» 

38) I 

WORKING file contents 
(WORK AREA) 

external DATA 

SETS USED: 



NAME 

LUN 


I/O OPERATIONS PERFORMED 

OUT 

NPRINT 


URITE(D£BUG ONLY) 

PRM 

NPRM 


READ 

URK 

NUI K 


WRITE 

programmer: D 

. e. BOLAND* 

CSC* 

MARCH 1981 


Subroutine WFILE 
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SUBROUTINE UIND(POS»TPDUf lERR) 

PURPaSET 

MIND CONTROLS THE CALCULATION OF THE WIND 

COMPONENTS U AND U AND THE MIND SPEED AND DIRECTION US AND UD* 

method; 

RETRIEVE D/8D COEFFICIENTS AND LATITUDE GRADIENTS FOR 
this POINT FROM COMMON /INTCOM/. 

2* CALL DIGRAD TO CALCULATE THE LONGITUDINAL PRESSURE GRADIENT* 
3* CALL OUIND TO CALCULATE THE 6E0STR0PHIC WIND* 

4* RETURN* 

ARGUMENT LIST: 


ARGUMENT 

TYPE 

I/O 

DESCRIPTION 

POS 

R 

I 

ARRAY CONTAINING TRAJECTORY POINT: 
POS(l)=ALTITUDE 
POS(2)=LATITUDE 
P0S(3)=L0CAL SOLAR TIME 
P0S(4)-L0NGITUDE 

TPDW 

R 

I/O 

ARRAY CONTAINING CALCULATED PARAMETERS 
TPDW(1)=TEMPERATURE <K DEGREESr IN) 
TPDW( 2) -PRESSURE (N/M**2f IN) 
TPDW<3)*DEN8ITY (KG/M««3f IN) 
TPD«(4)=E-W WIND U <M/SEC» OUT) 
TPDU(5)=N-S WIND V <M/SEC> OUT) 
TPDU(6)»UIND SPEED (UNUSED) 
TPDU(7)»WIND DIRECTION (UNUSED) 

I ERR 

I 

0 

ERROR flag: 
=^0? NO ERROR 


CALLING subroutines: MODELS 

SUBROUTINES CALLED: DIGRADrOUIND 

COMMON BLOCK PARAMETERS** 


COMMON NAME 

PARAMETER 

I/O 

DESCRIPTION 

/USECOM/ 

IDBG 

I 

FLAG FOR DEBUG 


BOUND 

R 

BOUNDARY BETWEEN UPPER AND LOWER 

/INTCOM/ 

JZ 

I 

ATMOSPHERIC MODEL SEGMENTS 
POINTER TO NEAREST ALTITUDE TO 


VALINT 

R 

POS(l> IN ARRAY RETURNED BY DFLINT 
6X38 ARRAY CONTAINING VALUES OF 

/FILCOM/ 

NPRINT 

1 

GRADIENTS AND D/SD COEFFICIENTS 
PRINTED OUTPUT LUN 

E)(TERNAL DATA 
NAME 

SETS used: 
LUN 

I/O 

OPERATIONS PERFORMED 

OUT 

NPRINT 

URITE(DEBU6 ONLY) 


programmer: d* e* Boland^ csc 


Subroutine WIND 



SUBROUTINE XALFA ( ISEGM f Z r ZREF » PBOUND » * lERR ) 

purpose; 

XALFA COHPUTES INTEGRALS NECESSARY FOR ERROR ESTIMATIONS 

method: 

1. COMPUTE INTEGRALS USING SUBROUTINE SIMP (INTEGRANDS 
PROVIDED BY SUBROUTINE GRAND)* 

2* RETURN 

ARGUMENT LIST: 

ARGUMENT TYPE I/O DESCRIPTION 

ISEGM I I SEGMENT NUMBER MINUS ONE 

Z R I ALTITUDE (KM) 

ZREF R I REFERENCE ALTITUDE (KM): ALTITUDE FOR 

THE TOP BOUNDARY OF CURRENT SEGMENT 
PBOUND R I PRESSURE AT ZREF (N/M*«2) 

ALFA R 0 AN ARRAY OF INTEGRALS 

lERR I 0 ERROR FLAG 

CALLING subroutines: POLFIT* DCERR 

SUBROUTINE CALLED** SIMP 

COMMON BLOCK PARAMETERS USED*. NONE 

EXTERNAL DATA SETS USEDt NONE 

PROGRAMMER*. T. LEE» CSC» SEPTEMBER 1981 


Subroutine XALFA 
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SECTION 4 - LAIRS MATHEMATICAL SPECIFICATIONS 

The mathematical details necessary for the LAIRS implementa- 
ion are presented in this section. Section 4.1 deals with 
the default atmospheric models. Section 4.2 discusses ad- 
justments to the models on the basis of meteorological meas- 
urements, including least squares fitting procedures for ver- 
tical profiles, followed by the computation of atmospheric 
parameters in Section 4.3. Section 4.4 discusses error cal- 
culations, Section 4.5 describes the interpolation methods used 
in the LAIRS Program, and Section 4.6 describes a method of 
adjusting translational coefficients using meteorological data. 

4 . 1 DEFAULT ATMOSPHERIC MODEL 

The LAIRS default atmospheric model consists of three parts: 
the Edwards Air Force Base Reference Atmosphere, extending 
from the ground to approximately 25 kilometers; the 1972 
Cospar International Reference Atmosphere (CIRA) , extending 
from 25 kilometers to 110 kilometers; and the 1971 Jacchia- 
Roberts atmospheric model above 90 kilometers. Between 90 
kilometers and 110 kilometers, either the CIRA or the Jacchia- 
Roberts model can be used. These models are described in the 
following sections and in References 2, 3, 5, and 6. 

4.1.1 THE JACCHIA-ROBERTS MODEL 

The Jacchia-Roberts atmospheric model determines the atmos- 
pheric mass density, , and the number density of individ- 
ual species, n^, using the following two equations: 

Iz ^ ' •L&\ooKm (4-1) 

i Z^IOoKm (4-2) 
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where 


0 !^ = thermal diffusion coefficient of the species i 
= molecular weight of the species i 
M = mean molecular weight 
g = gravitational acceleration 
R* = universal gas constant 

Equation (4-1) is the barometric equation and applies to 
the region of homogeneously mixed atmosphere below 100 
kilometers^ and Equation (4-2) is the diffusion equation 
for the number density of each individual species and applies 
to the atmosphere above approximately 100 kilometers, where 
the lighter species separate from the heavier species by 
the process of molecular diffusion. 

The temperature, T, in Equations (4-1) and (4-2) must be 
known as a function of Z , Temperature-height profiles were 
empirically obtained in the Jacchia 1971 model and will be 
slightly modified according to the formulation given by 
Roberts so that the density equations can be integrated 
analytically (References 6 and 7) . 

The temperature profile to be used in the Jacchia-Roberts 
atmospheric density model takes the form of a fourth-order 
polynomial for heights from 90 kilometers to 125 kilometers 
and an exponential function of Z for heights above 125 
kilometers . 


4 - 

Hz) =Tx+ (T.-TjUan fZ-zJ 
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(4-3) 
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This equation gives temperatures between Zq and Z^, and the 
coefficients a^ are determined from the following conditions; 


T (Z 

0^ 

= T * 
^0* 

T (Z 

) 

= T : 

X 

X 

/dT \ 

VdZ 

^Z=Zq 

= u 

/dT ^ 

1 dZ 

^z=z 

X 

= 

(- 

1 ) 

= 232 

\ dZ 

/ 

z=z 



= T^ : minimum temperature 


"x-"0 


(4-4) 


an empirical 
relation 


= 2^2 (VTq) = 0 


X 


The four coefficients a^ (i = 1~4) can be expressed in terms 
of Zq, Z^, Tq , and T^ as follows: 

a, = /.?/ (Zk-%c) 

dz 0 (4-5) 

Ui = -i.i /czf-t.f 
= -0-8/ (Zx- z.f’ 


The values of Zq and Z^ are fixed at 90 kilometers and 125 
kilometers, respectively. The value of Tq was 183 K in the 
Jacchia 1971 model and was changed to 188 K in the Jacchia 
1977 model. The quantity T is determined from the value 
of Too (Reference 5) . 


^ "t b I tfy 
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where 


a = 371.6678 


b = 0.0518806 
c = -294.3505 
k = -0.00216222 

For heights above 125 kilometers, Roberts used the following 
function to define the temperatures: 

T«) =T„ - (T.- T. ) H ) <.-7, 

where G = temperature gradient at Z 

R = 6356.766 kilometers = mean radius of the Earth 
a 

A temperature profile defined by Equations (4-3) and (4-7) 
is continuous at Z = Z regardless of the choice of G . 

The temperature gradient will be continuous at Z = Z if 
G = 1.9(T -T„)/(Z -Z,.) . However, in the Jacchia-Roberts 
density model implemented in the Goddard Trajectory Deter- 
mination System (GTDS) , G was determined as a function of 
T^ such that the resulting density profiles obtained using 
Roberts' temperature profiles gave the best least-squares 
fit to Jacchia's tabulated density data. In this case, 
the temperature gradient is not continuous at Z = Z . 

In addition to temperature profiles. Equations (4-1) and 
(4-2) require a knowledge of g(Z) and M(Z) . The following 
expressions were used for these quantities: 




(4-8) 
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2 

where = 9.80665 meters per second and 

S 

T\L%) = Y^ Cni'L-9o) (90<Z<loo) (4-9) 

n=o 


The numerical values of in this equation are given in 
Reference 5. Roberts analytically integrated Equations 
(4-1) and (4-2) using temperature profiles given by Equations 
(4-3) and (4-7), the gravitational acceleration, g(Z), and 
the mean molecular weight, M(Z) . The quantities g(Z) and 
M(Z) are defined by Equations (4-8) and (4-9), respectively. 
The solutions obtained by Roberts, however, contain rather 
lengthy algebraic expressions and will not be given in this 
report. (See References 6 and 7 for further details) . 

Each temperature profile is labelled by the so-called exos- 
pheric temperature, which sets the upper boundary value 

for the temperature. Too is computed in three steps in the 
Jacchia-Roberts model. First, the nighttime minimum exos- 
pheric temperature, T^ , is obtained from the solar 10.7- 
centimeter flux level, using the relationship 

Tc = 379' + 3’. 2^ ff.a ( F„.-j'F,o.-j ) (4-10) 

The second step involves computing the exospheric 

temperature that includes the diurnal variations, from T , 

c 

i .e . , 


~^NL ( ~T'D 


(4-11) 
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where 


Td = Tc ( It Rcos’”ite-Sl) 

■>N = Tc ( 1 +Rsin’"^le+S| ) (4-12) 

X =■ H T p T f> s!n ( H + D ) 


In these equations, H is the local solar time, 6 is the 
latitude of the subsatellite point, and 6 is the solar 
declination. The values of the parameters appearing in 
Equation (4-12) are: 

m = 2.2 

3 - -37° 

p =6° 

Y = 43^ 

The contribution to Too from the geomagnetic activity effect 
is given by 



28°Kp -t 0*03 (zyzooKm) 

f4-°Kp T o'.oze!^’’ c Z<ZoofCm.') 


(4-13) 


where is the 3-hour geomagnetic planetary index. An 
average time lag of 6^?7 was used for the time lag associated 
with the atmospheric variation resulting from the geomagnetic 
disturbance . 

The total exospheric temperature is given by 


^00 ~~ ( ^Tco ) 




(4-14) 
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The value of T«, uniquely determines the temperature profile 
to be used in the computation of the density, which will be 
further corrected for the semiannual and seasonal-latitudinal 
variations . 

The Jacchia-Roberts atmospheric density model described above 
gives density and temperature. The pressure values are com- 
puted using the following relations; 

Ptr fz) = IZ*T(Z) T rlji(Z) (Z>looKm) (4-15) 

P:yR(z) = R.*TCz)P^K('z)/i^(^z) (Z^looKm) <4-16) 

When the mean molecular weight is desired at various alti- 
tudes, the following relations can be used: 

= I?T(Z) 

In these equations, denotes the number density of an 

individual species given by Equation (4-2) , and denotes 
Avogadro's number. 

4.1.2 THE CIRA 1972 AND ERA MODELS 

The lower atmospheric part of the CIRA 1972 model (Reference 
3) gives tabulated atmospheric parameters, i.e, temperature, 
pressure, density, and wind velocities at intervals of 5 
kilometers in height, 10 degrees in latitude, and 1 month, in 
time. Furthermore, the meteorological data, used in the CIRA 
1972 model has a high concentration near a latitude of 


(4-17) 

(4-18) 
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30 degrees north, which is in the neighborhood of the 
Shuttle reentry path. The lowermost 25 kilometers of the 
atmosphere are not covered by the CIRA. 1972 model. For this 
region, the Edwards Air Force Base Reference Atmosphere 
(ERA) 1975 model (Reference 2) , which was designed with the 
Shuttle reentry in mind, is recommended as the default model. 

The combination of the CIRA 1972 and ERA 1975 models provides 
a general default lower atmospheric model for latitude 30 
degrees north during the major portion of the Shuttle re- 
entry and a more specific default model during the final 
descent phase into Edwards Air Force Base. These models 
are implemented in LAIRS in the form of tabulations on a 
permanent file. 

4.1.3. DIURNAL, SEMIDIURNAL, AND LATITUDINAL VARIATIONS 
IN THE LOWER ATMOSPHERE 

The diurnal and semidiurnal effects of CIRA 1972 can be 
incorporated in the following approximate manner. Let 
f(z,t) be an atmospheric parameter that is a function of 
altitude and local solar time, where f(z,t) is given approx- 
imately by 


-f rz.t) = -tAi(z) cosj^(t-t,u))^ 
-t kzCz) c.os^ Ct 


(4-19) 


where t is expressed in hours and 

= diurnal average of f(z,t) 

= diurnal amplitude at Z 

= semidiurnal amplitude at Z 

= time of the maximum diurnal effect at Z 

= time of the maximum semidiurnal effect at Z 


Aq (z) 
A^ (z) 
A 2 (z) 
t^ (z) 


t2(z) 
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The values of A^i and are supplied graph- 

ically in the CIRA 1972 publication (Reference 3) and are 
tabulated on the LAIRS permanent file at 5-kilometer inter- 
vals. (The temperature coefficients given by CIRA are in- 
complete and have been interpolated in LAIRS to smoothly 
join the Jacchia model.) Because the values of the param- 
eters tabulated by CIRA are not the average values Aq that 
are required by Equation (4-19) but, rather, are values 
biased toward noon, the calculation of diurnal and semi- 
diurnal effects in LAIRS takes the following form. First, 
the average A^(Z) is determined from the tabulated or 

interpolated f(Z, t ), where t represents the time of 

R R 

the reference tabulation (in this case, noon): 


Ao(Z) = -f CZ.,tR) -A,fZ)COS{X(tR-t,(Z))} 

(4-20) 


Once Aq(Z) has been determined, it is substituted into 
Equation (4-19) to yield f(Z, tg) , the value of the param- 
eter f at any input time, tg, corresponding to a Shuttle 
position. 

The latitude variations of the atmospheric parameters have 
been incorporated in LAIRS in the form of latitude gradients 
stored on the permanent file at 5-kilometer intervals. 

These gradients, g(f(Z, 6)), are simply the first deriva- 
tives in latitude of temperature, pressure, and density, as 
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determined from the tabulated CIRA values in the region of 
30 degrees north latitude: 


^Cf(z,eR)) = 


9f(Z,e) 

'96 


e = 0R =30° 


•f (Z.,40’)— 

Zo‘ 


(4-21) 


Below 25 kilometers, the tabulated values of f (temperature, 
pressure, and density) used in the above calculation were 
taken from the Handbook of Geophysics (Reference 4) , which 
gives tabulations at latitudes of 15 degrees and 45 degrees 
rather than at 20 degrees and 40 degrees. 

Using these tabulated gradients, extrapolation of the value 
of any atmospheric parameter f to the Shuttle latitude, 

6 , is easily achieved using the following relationship: 


f(z,6s) =f(z,0R>[i t 5(-f(z,eR))*(es-eR)] <4-22) 


This calculation is generally used only below 90 kilometers, 
where 6 is always relatively close to 30 degrees north 

O 

latitude, the default value of the reference latitude. 

Although the values of the diurnal and semidiurnal coeffi- 
cients and latitude gradients are tabulated in the LAIRS 
permanent file at 5-kilometer intervals (as are the values 
of the parameters themselves), they, like the parameters, 
are always interpolated to yield a smooth output. 
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4.2 ADJUSTED ATMOSPHERIC MODEL 


Both the upper and lower atmospheric models can be adjusted 
on the basis of meteorological data. For the lower model, 
which in its default form consists of simple tabular data, 
adjustment is achieved by fitting a set of polynomials to 
the meteorological data (in accordance with the gas law and 
hydrostatic equilibrium, unless the user selects otherwise) . 
Alternatively, the user may elect simply to interpolate from 
the meteorological data. For the Jacchia-Rober ts model, the 
boundary values of temperature and density at 90 kilometers 
are adjusted to agree with the data and to be continuous 
with the lower atmospheric model. Section 4.2.1 discusses 
the Jacchia-Roberts model adjustments, and Section 4.2.2 
describes the methods used in the polynomial fitting of the 
lower atmospheric parameters. 

4.2.1 JACCHIA-ROBERTS MODEL ADJUSTMENTS 

Atmospheric parameters computed using the Jacchia-Roberts 
model depend on the temperature and density values at the 
lower boundary (90 kilometers) of the Jacchia-Roberts 
model. In LAIRS, these boundary values are reset in such a 
way that the continuity of the computed atmospheric param- 
eters is satisfied at the boundary. A similar readjustment 
of the density value at 100 kilometers is also needed be- 
cause the Jacchia-Roberts formulation has an internal 
segment boundary at 100 kilometers that is affected by 
adjustments made at 90 kilometers. Boundary value read- 
justments at 90 kilometers are carried out by replacing the 
original Jacchia-Roberts values with the temperature and 
density values at 90 kilometers obtained using the polyno- 
mial fitting procedure for the lower part of the atmosphere 
(below 90 kilometers). In the case of temperature, the re- 
placement is simple; T(90), the Jacchia-Roberts default 

boundary temperature, is replaced by T (90), the polyno- 

P 

mial adjusted boundary temperature. 
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In the case of density, however, simply replacing the 

Jacchia-Roberts default boundary density p(90) with the 

polynomial adjusted density, p (90) , is insufficient, 

P 

because the Jacchia-Roberts model makes additional correc- 
tions (for geomagnetic, seasonal, and latitudinal effects) 
to the density. This means that the final density value 
output by the model at 90 kilometers is not equal to the 
input boundary value at that height, but rather is a cor- 
rected density of the form 


fj90) = fpC9o)o( 


(4-23) 


where a denotes a correction factor that has a value close 
to 1 and is assumed to be independent of pp(90). 

LAIRS uses a two-step procedure to ensure that density is 

continuous across the boundary. First, the Jacchia-Roberts 

model is run with the polynomial adjusted density p (90) 

P 

as the boundary value. From the output p^(90), the 
value of a is determined: 


0 ( 


fp(9o) 


( 4 - 24 ) 
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using the factor a, a new boundary value of density 
is defined as 


Second , 


P„(90) , 

It 



fp(9o) 

o< 


(4-25) 


The quantity Pp(90) is then used as the density 
boundary value in all future Jacchia-Rober ts calculations. 

It yields a Jacchia-Rober ts output density value at 90 kilo- 
meters that is equal to the polynomial density value 
Pp(90) . 

The 100-kilometer density value is linearly proportional to 
the 90-kilometer density value and is always computed using 
the Jacchia-Roberts density value at 90 kilometers. Thus, 
the 100-kilometer density is readjusted by scaling the density 
value at 100 kilometers obtained using the Jacchia-Roberts 

I 

model with the ratio of Pp(90) , which is the new boundary 

density at 90 kilometers, over p (90) , which is the poly- 
pe 

nomial adjusted density at 90 kilometers. 

4.2.2 POLYNOMIAL MODEL ADJUSTMENTS 

The idea behind the adjusted polynomial model is that the 
best approximation for temperature, pressure, and density in 
the entire lower atmospheric region (below 90 kilometers) 
can be achieved by combining the available meteorological 
data into a reference profile that can then be subjected to 
least squares polynomial fitting. The entire altitude range 
is divided into four segments, as shown in Figure 4-1. The 
Jacchia-Roberts model is used in the uppermost segment and 
the polynomial model is used in the lower three segments. 
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ALTITUDE (kilometers) 



Figure 4-1. Definition of Altitude Segments 
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In each of the three segments of the polynomial model, 
third-order polynomials are fitted to temperature and to the 
logarithms of pressure and density. The segment boundaries 
may be chosen by the user; their default values are 90 kilo- 
meters, 65 kilometers, and 25 kilometers. 

input meteorological data are received in the form of sep- 
arate vertical profiles taken at different times and dif- 
ferent places near the Shuttle reentry trajectory. Before 
these data can be combined for fitting, they must be trans- 
formed to a common latitude and local solar time. This is 
achieved by means of the diurnal and semidiurnal coeffi- 
cients and latitude gradients discussed in Section 4.1.3. 

The reference latitude and local solar time may be selected 
by the user, but default values of 30 degrees north latitude 
and local noon are supplied. Once the reference profile has 
undergone fitting, the diurnal and semidiurnal coefficients 
and latitude gradients are used to translate »from the refer- 
ence position to the Shuttle position. 

Two different procedures are used to fit the reference pro- 
files for the lower atmosphere. In the first scheme, tem- 
perature, pressure, and density data are directly fitted to 
a set of polynomials. There are three different polynomials 
for each segment: one for temperature, one for pressure, 

and one for density. In the second scheme, only the temper- 
ature polynomials are constructed, using temperature, pres- 
sure, and density data or using any combination of these 
three data types. In this procedure, the pressure and den- 
sity values are derived from the temperature polynomials 
using the ideal gas law and the barometric equation. De- 
tails of these two procedures are given in Sections 4. 2. 2.1 
and 4. 2. 2. 2. 
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4. 2. 2.1 Independent Polynomial Fittings for Temperature, 
Pressure, and Density Profiles 

In this procedure, temperature, pressure, and density poly- 
nomials are constructed independently using the corre- 
sponding meteorological profiles. The continuity condition 
is not imposed on the fitted atmospheric parameters at the 
boundary between two segments. 

Let F(Z) be the polynomial function of Z to be fitted to the 
meteorological profile f (Z^) . For pressure and density, 
f(Z^) will be the natural log value of the observed data. 


3 

FCZ) = />„ + X (4-26) 

•f(Zn) =r temperature 

= for pressure (4-27) 

for density 

Zq r appearing in the expression for F(Z), denotes a refer- 

I\ 

ence altitude. In the current version of LAIRS, Zj^ is 
defined as the upper boundary altitude of the segment. Then 
Aq becomes the upper boundary value of the computed atmos- 
pheric parameter. A set of coefficients (Aq, A^, A^) 

will be independently determined for temperature, pressure, 
and density. 
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In the least squares fitting procedure, the sum of the 
squares of the weighted residuals is minimized, i.e,. 


M 


Q(ft) 


n-1 


(4-28) 


where A denotes the array A. (i = O'V/S) and to is the 

1 n 

weighting factor for the residual, - F(Z^). 

n n 

The weighting factors used here are computed using the meas- 
urement uncertainty Af(Z_) 

j ' n 


“ y^-fCZn) 


(4-29) 


where Af (Z^) is the measurement uncertainty associated 

with the nth observation. For the case of temperature, 

Af(Z ) is taken to be a , the measurement standard 
n n 

deviation. (If no standard deviation is reported, a default 
value of 5 percent is assumed; however, the user may change 
this value via the DCEDIT keyword card.) For pressure, the 
uncertainty for 5,n(P^) is defined as 

AffCS.) = 1 + ] (4-30) 


where a and P denote the o and P values at Z 
n n 


Z 


n* 
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This expression was derived from the following simple con- 
sideration. If P is expressed as 


p = P + AP 


the corresponding expression in terms of £,n(P) will be 


inP = /nP f^rt[l + ^ ] 


Thus, iln(l +(AP/P)) can be interpreted as the uncer- 
tainty associated with the measurement of in(P) . 

The polynomial coefficients A^(i = 0~3) will be deter- 
mined when Q (A) is minimized with respect to A^ . This is 
done using subroutine LSQPOL, which is available in the CDC 
system library at LaRC. Mathematical details of this proce- 
dure will not be given here, since it is a special case of 
the procedure presented in Section 4. 2. 2. 2. 


4-18 



4. 2. 2, 2 Det-ermination of Temperature Polynomials Using a 

Differential Correction Procedure With Temperature, 
Pressure, and Density Data 

In this section, a differential correction procedure that 
determines a set of temperature coefficients and a boundary 
pressure value is described. This procedure involves mini- 
mizing the sum of the squares of the weighted residuals of 
temperature, pressure, and density with respect to temper- 
ature coefficients and the boundary pressure values. 

The sum of the squared weighted residuals to be minimized is 
given by 


Qis) = (4-32) 


fT) fP) ^D) 

where Q' , Q' , and Q' ' denote the contribution to Q from 
temperature, pressure, and density, respectively. The 
solve-for parameter vector S is defined as a column vector 
consisting of temperature coefficients and the boundary 
pressure value. The transpose of S is given by 



= r I , Aa , 3 , ^ Po ) 


(4-33) 
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The temperature polynomial is defined by Equation (4-26) . 
Q^^Ns), and are given by 


Q^^^CS) = Z TobsCZ- ) -T(z„) ) ] 

Y\ 


vn 


In these equations, the following notations were used. 


3 

TCz) = T(z,s)= 4, + Z AkCzc--z.) 



k =:1 


4-34a) 


4-34b) 


4-34c) 


4-35) 
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PoBS f 


(4-36a) 


where 
the g 
puted 


where 

R is 
e 


foBsfZ'n)= in CfoBS^Z-.)) 


P(Z„) = In P. - j 


2in 


1IBU± dz. 

T(Z.) 


f Cz«) 



P,M 

R'TCz„-) 



RT(z) 


dz 


/w 

g(Z), which appears in the integrands of P and p, 
avitational acceleration at the altitude Z and is 
using an approximate equation: 



g^ is the gravitational acceleration at Z = 0 and 
an effective radius of the Earth. 


-36b) 


-37a) 


-37b) 


is 

com- 


-38) 
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Below 90 kilometers, a constant mean molecular wieght, M, is 
used throughout. 

Using two column vectors, f and F(S), and a diagonal matrix 
W, the function Q(S^) can be expressed as 


Qcs) = (f -FCS))* w -F(S)) 


(4-39) 


where f, F, and W are defined as follows: 


f ^ 


fops ^ j PoBS ^^2 ^ f ■ ' / P 005 ^ 2; ^/ ) ; 

fops * h Jobs ’ y ?0BS f ^ N'' ^ } 


(4-40) 


F(£/={ T(z.),T(xO, •••.T(Zn') ; 

p(z/), p'czi) , ‘ . P ; 

f (zr), fczi), • •• . f 


(4-41) 
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1 4 ^ ^ rV 


Wii = 

■= C ^^i+^i' 

(4-42) 

*y N/f k/V t ^ t ^ aJ^ 

Wij — o C <.’ j ) 

where N^ denotes the total number of observations included, 
i.e., N, = N + N' + N", and {Z^;n = 1~N>, fz' :n = 1-N ' } , and 
{Z"^:n = 1~N"} denote arrays of altitude points associated 
with temperature, pressure, and density observations, re- 
spectively. These three arrays are taken to be identical in 
LAIRS. 

In general, F(S) is not a linear function of S. The differ- 
ential correction procedure presented here involves linear- 
izing F(S) using the Taylor expansion: 


pcs) = PCs.; + B^s 


(4-43) 


where 


B = 



X modrix 


(4-44) 
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and 


AS = 5 -So 


(4-45) 


N. and N denote the total number of observations and the 
t s 

total number of solve-for parameters, respectively. 

Using Equation (4-43) in Equation (4-39) , a Q-function that 
is quadratic in AS is obtained, i.e.. 


- E^sf V\/ (^f - B4S) 


(4-46) 


where 


= -P - F(So) 


(4-47) 


The function Qj^(AS) is readily minimized with respect to 
AS. The vector AS that minimizes Q (AS) is given by 

J_j 


W B) 


W 


(4-48) 
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S can then be determined as follows; 


-S ■= So •f' 


(4-49) 


This process is repeated with S as a new until AS 
becomes smaller than a prescribed tolerance limit. If F(S) 
is a linear function of S, i.e., F(S) = Fq + BS, then the 
S that minimizes Q(S) is given by 


s 


-I 


( B^W B) 3'*^ W 


(4-50) 


where 


Af = f - Fo 


In this case, the complete answer is obtained in one step 
and no iteration is required. The procedure used in Sec- 
tion 4. 2. 2.1 corresponds to the case where F(S) is a linear 
function of S. 
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In order to compute AS in Equation (4-48) , the matrix B 
must be evaluated. From the definition given in Equa- 
tion (4-44), component of B, is given by 


B 


nk 



3 F ^ ^ 'I 

/ S =?c 


(4-51) 


There are three different types of partial derivatives of 
this kind. The first type consists of the partial deriva- 
tives of the model temperature given in Equation (4-35) with 
respect to S. (See Equation (4-33) for the definition of 
S) , i .e. , 


B 


nK 


3T(Z„.s) 


V 'aSfe J 

1 

(Z.-Zn)^ 

: k = 1 ~ 

1 

"4- 

II 

0 

: « = 5 


(4-52) 


The second type of partial derivatives consists of the model 
pressure partial derivatives with respect to S. Using the 
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definition of P(Z) given in Equation (4-37a) , these partial 
derivatives are computed as follows; 


Bn'k. “ ( 

- s 


'dP(z'n',s) 


^r>' g(z) M 


/ S = S. 


S = S, 

.A 

2o 




/3S, 


/i'S. 




(4-53) 


= 1 


; k'S” 


Similarly, the third type of partial derivatives, the density 
partial derivatives, are given by 


B ^ 


'9P<:z;^3) n 

35^; ys=s. 

_i (2J(^hV] + 1 ' 

T(z>,s^l -3 Sk r *[ t(z 


^n" 


.=.>] Ws.!. 


(4-54) 


^ i 


j k= 1 

; 4c- 5 
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4.3 COMPUTATION OF ATMOSPHERIC PARAMETERS 


The final readjusted model vertical profile of an atmospheric 
parameter is given by Equation (4-26) in the case of inde- 
pendent polynomial fitting. The same equation will be used 
for temperature, pressure, and density, with an appropriate 
set of coefficients for each segment. In the case of com- 
bined fitting of temperature, pressure, and density data, 
the readjusted temperature profile is given by Equa- 
tion (4-35). Pressure and density at a particular altitude 
are computed from the temperature profile and the boundary 
pressure value using the gas law and the barometric equa- 
tion. Atmospheric parameters computed in this manner cor- 
respond to a reference latitude, 0^, and a reference 
local solar time, t^. Thus, an atmospheric parameter 
computed at (Z, 0 , t ) is denoted as F(Z, 0 , t ). The 

K. i\ R R 

computation of the same atmospheric parameter at a Shuttle 
position (Z , 0 , $ ) and a Greenwich mean time of t 
necessitates diurnal, semidiurnal, and latitudinal transla- 
tions as described in Section 4.1.3. Equations (4-19), 

(4-20), and (4-22) are used in this process (the values of 
0 and t in the equations are now those associated 
with the location of the reference profile) . Following 
these corrections, temperature, pressure, and density may no 
longer strictly obey the gas law because of uncertainties in 
the translational coefficients and gradients. The percent- 
age deviation from the gas law of these translated values is 
computed by deriving the value of the gas constant R* from 
these values and comparing it with the known value; 



PM 

PT 


(4-55) 
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X LoO 


( 4 - 56 ) 


= - R* 

R* 


where is the derived value of the gas constant and y 
is the percent deviation, which is reported in the LAIRS 
output. Temperature, pressure, density, and mean molecular 
weight are designated by T, P, p, and M, respectively. To 
ensure a set of parameters that are consistent with the gas 
law, LAIRS recalculates the temperature from the translated 
values of pressure and density as a final step (temperature 
was chosen because the diurnal coefficients for temperature 
are the least well known) simply by 




( 4 - 57 ) 


The user may set a keyword card (LOWMOD) to skip this final 
temperature correction if the correction is not desired. 
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In addition to the percent deviation, LAIRS prints out the 

values of the mean molecular weight, M, and the pressure 

scale height, H . The latter is calculated as 
P 


H 


? 


R*T 

M 3 


(4-58) 


where g is the acceleration of gravity and is given by Equa- 
tion (4-8) . The mean molecular weight reported above 
90 kilometers is the value given for the Jacchia-Rober ts 
model by Equation (4-9) , while below 90 kilometers it is 
simply the sea level value of 28.9644 grams/mole. 

Two other parameters are derived by LAIRS from the calcu- 
lated values of temperature, pressure, and density. These 
are the Shuttle Mach number m and the estimated onboard 

pressure measurement P . The Mach number is merely the 

n 

ratio of the Shuttle speed, V, to the speed of sound c: 


yn 



(4-59) 


For air, y = 1.4. The onboard pressure measurement is 
derived via a set of equations that are valid for a body 


4-30 



with the aerodynamic characteristics of the Shuttle, 
the coefficient of pressure, C , is calculated as 

p 


Cp = i.84 ( cosat’-oO)'^ 


where a is the angle of attack of the Shuttle and 


n = 2.7 — 


i. 5 


The dynamic pressure, q, is then given by 


f 


± 

z 



z 


and the onboard pressure measurement is given by 

Pjy = P ■+ ^Cp 


First , 


( 4 - 60 ) 


( 4 - 61 ) 


( 4 - 62 ) 


( 4 - 63 ) 
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where P and p are the free stream pressure and density 
(i.e., the values of pressure and density determined by 
LAIRS, respectively). As the Shuttle speed decreases during 
reentry, these equations break down between Mach 2 and 
Mach 1; however, they are reported down to Mach 1 by LAIRS. 

Finally, LAIRS must produce values for the wind components. 
These are not modeled in any way; they are merely interpo- 
lated from the meteorological data with no attempt at diur- 
nal or latitudinal translation. It should be recognized that 
while this approach cannot guarantee accurate values for the 
wind components, it was the only viable alternative available 
(see Reference 1). However, the option of calculating wind 
components using the geostrophic wind equations is provided 
in LAIRS (see Reference 1) . This technique requires the use 
of diurnal and semidiurnal coefficients and latitude gradients, 
which are, of course, default values. 

4 . 4 COMPUTATION OF UNCERTAINTIES ASSOCIATED WITH ATMOSPHERIC 
PARAMETERS 

An atmospheric parameter is computed using a function F of 
the altitude, Z, and the solve-for parameters, S: 

F = F (Z.,s) (4-64) 

The solve-for parameters, S, are obtained through a differ- 
intial correction procedure. The variance-covariance matrix 
of S is available from the differential correction 
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procedure and used to obtain an approximate variance of F 
V (F) , where 


vCf) 



( 4 - 65 ) 


In this equation, C(S) represents the variance-covariance 
matrix, and (8F/3S) represents a column vector of par- 
tial derivatives of F with respect to S computed using the 
converged differential correction solution. In a normal 
differential correction procedure, five parameters are de- 
termined for the top segment of the lower atmosphere, but 
only three parameters each are determined for the middle and 
bottom segments. The upper boundary temperature and pres- 
sure of the two lower segments are not solved for but are 
computed using the solutions for the upper segment. Thus, a 
normal differential correction process in LAIRS yields a 
5x5 variance-covariance matrix for the top segment and a 
3x3 variance-covariance matrix for each of the two lower 
segments. Contributions to the variance of F from uncer- 
tainties associated with the computed boundary temperature 
and pressure can be obtained in the following way. First, 
the variance-covariance matrix of the fixed boundary temper- 
ature and pressure is computed using the solutions of the 
upper segment. Then, the direct sum of the 3x3 variance- 
covariance matrix of the differential correction process and 
the computed 2x2 variance-covariance matrix is formed. The 
resultant 5x5 matrix will be used in Equation (4-65) to ob- 
tain V (F) . 
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The estimated uncertainty of F, AF, is obtained from 




( 4 - 66 ) 


and printed out in the LAIRS Output Report, Part I. The 
explicit forms of F(Z,S), 9F/3S, and C(S) that are used 
in LAIRS are listed below: 


S* = ( fl. , fis . fi. .jtnP.) 

fl = 

TCz) = A, + + A3 fz.-z) 
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where C{A*) = 3x3 variance-covariance matrix of 
A' = (Aj^, h^ 2 t and 

C ( 3) = 2x2 variance-covariance matrix of 
B = (Ag, to(Pg) and 



where s’ indicates the five-parameter array of the upper sec- 
tion. 

Somewhat different types of error estimations for the poly- 
nomial fitting scheme discussed in this section can be ob- 
tained using the average weighted residuals (0-C) and the 
weighted root-mean-square (WRMS) values for the temperature 
residuals, the log pressure residuals, and the log density 
residuals. First, the standard deviation of the weighted 

*See Section 4.2.2 for definitions of B and W. 
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residuals for temperature, log pressure, and log density is 
given by: 


Sfavidard Devia+ion o-f Residuals 

( 4 - 67 ) 

{Wei^nteJ RMS^ — (A\zer^e Residua]) 


Second, the WRMS values of the residuals can be used to 
estimate an average percentage deviation of the computed 
temperature, pressure, and density from the corresponding 
observed values, i.e.. 


St 

T 

iP 

p 

IL 

r 


. . (T(~r) 

^ WRMS(T)X — 

^ 6xp / w RM S ( ^np) x/n ( 1 + — i 

P ■' 

P) \ 'J 

<5; exp { WRMS (^mp)x ^ri { l + 


( 4 - 68 ) 


where a represents the uncertainty assumed for each meas- 
urement type. 

In a typical differential correction run, a(F)/F (where F 
denotes temperature, pressure, or density) is assumed to be 
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0.05 and the WRMS values are in the neighborhood of 0.5. 
Thus, for all three cases, the approximate percentage de- 
parture of the computed value from the observed value is 
given by 


F 


WRMS 



(4-69) 


If a(F)/F is not a constant in a given segment, an average 
of this quantity for the segment can be used in this error 
estimation . 

4 . 5 INTERPOLATION METHODS 

The LAIRS Program incorporates two different types of inter- 
polation, depending upon the type of input data used. This 
input data can be either the default atmospheric data or 
observed meteorological data. 

4.5.1 DEFAULT DATA 

The default interpolation mode uses data supplied on the 
atmospheric data working file. Either a first-order or 
second-order Lagrangian interpolation in altitude can be 
used by calling the LaRC CDC library subroutine lUNI. A 
matrix consisting of values of temperature, pressure, den- 
sity, and wind components for six different altitudes is 
passed to subroutine lUNI for interpolation. If , Y^ , 

..., Y^ are the values of the function evaluated at the 
tabulated values Xj^, •••/ of the independent 

variable (altitude), then the interpolation equations used 
are as follows. 
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1 . 


First order; 


Vo = Vi 


- ( tc . -y.i) 


rxc, - xj 

■%, ^ ^t'+f 


2. Second order: 


(4-70) 


V = — 






( -XJ 


(4-71) 




J 


where the index i is chosen so that IX^ - X . I + 

0 1 

IXq - is a minimum. 
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These equations are only used for the temperature and the 
wind components. The pressure and density are interpolated 
using the following logarithmic relationship 


f(k) = 

‘ (4-72) 

Ac ^ k. ^ Ac+i 


where h is the altitude. 

4.5.2 METEOROLOGICAL DATA 

When meteorological profiles based on observational data are 
available, the user has several options. If more than one 
meteorological profile is available, then each profile is 
interpolated in altitude using the technique described in 
Section 4.5.1. After each profile has been interpolated in 
altitude, the points from each profile corresponding to a 
given altitude can be interpolated again in latitude and 
longitude or in latitude and local solar time. 

If the user chooses to interpolate in latitude and longitude, 
a simple univariate interpolation scheme is used to obtain 
temperature, pressure, density, wind speed, and wind direc- 
tion between the meteorological profiles. This scheme is 
based on the great circle arc lengths from the Shuttle 
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position to each of two meteorological profiles. The equa- 
tion used to obtain the interpolated values is as follows 



(4-73) 


where C 12 

Cs 

fl 


great circle arc length between the two meteoro- 
logical profiles 

great circle arc length from the first profile 
to the Shuttle position 

parameter value at first profile 


f 2 = parameter value at second profile 


The user can also specify a bivariate interpolation in lati- 
tude and local solar time for temperature, pressure, and 
density. In this approach, the longitude and Greenwich Mean 
Time (GMT) for each profile are converted to local solar 
times. The bivariate library subroutine IBI is called to 
perform the interpolation, which can be either a first-order 
or second-order Lagrangian interpolation. However, before 
the interpolator can be called, a complete table (two- 
dimensional array) must be formed; the elements of this 
table are values of the function being interpolated at a 
given altitude. The diagonal elements are the values on the 
input meteorological profiles. The off-diagonal terms are 
obtained by translation from the diagonal terms using de- 
fault latitude gradients, since the latitude gradients are 
more accurate than the diurnal/semidiurnal gradients. If 
first-order interpolation is used, a 2x2 table must be con- 
structed. If second-order interpolation is used, then a 
3x3 table (requiring three meteorological profiles) must be 
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constructed. The governing equations for conventional two- 
dimensional Lagrangian interpolation are as follows: 


m n. 


Z. = z Z (xo) CYo) F 

t = 0 


( 4 - 74 ) 


wi 

= 7 [ 




Xo -x„ 

Xi -Xfc 


< - 0 , 1 , ■■■,/»■ 


( 4 - 75 ) 



Yo-Yic 

Yj-r. 


j=c, t, ■ ■■ , n. 


( 4 - 76 ) 


where m = order of interpolation in the X direction 
n = order of interpolation in the Y direction 
and F(I,J) is the value of the function at (X(I), Y(J)). 

The values of the X and Y arrays must be in algebraically 
increasing order. For interpolation of meteorological data, 
the X array corresponds to latitude and the Y array corre- 
sponds to local solar time. 

AS indicated above, values from the meteorological profiles 
will normally form the diagonal of the bivariate function 
table. However, in certain cases, such as when the latitude 
of the middle profile (in ascending longitude) is greater 


4-41 


than the latitude of the third profile, two rows or two 
columns must be interchanged in order to satisfy the re- 
quirement of algebraically increasing values of the two 
independent variables. The table can still be completed 
using the default latitude gradients, but the interpolated 
values may not be as accurate. 

4 . 6 READJUSTMENT OF DIURNAL AND SEMIDIURNAL COEFFICIEN TS AND 
LATITUDE GRADIENTS 

It is possible to construct a revised set of amplitudes and 
phases for the diurnal and semidiurnal coefficients and a 
new set of latitude gradients of the atmospheric parameters, 
given meteorological data with a sufficient spread in lati- 
tude and local solar time. The readjustment procedure is as 
follows : 

1. Interpolate the values of the observed parameters 
to the altitudes at which the default coefficients 
and gradients are tabulated. 

2. Using the default diurnal and semidiurnal coeffi- 
cients, translate the interpolated parameters to a 
reference local solar time. 

3. Construct a new set of latitude gradients from the 
translated parameters . 

4. Using the new latitude gradients, translate the 
original interpolated parameters to a reference 
latitude . 

5. Construct a new set of diurnal and semidiurnal co- 
efficients from the translated parameters. Steps 2 
through 5 may be repeated. 

In practice, the above procedure involves some difficulties. 
Some of the default translational coefficients are poorly 
known and may not be close enough to the actual values to 
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yield good results. Even more important is the fact that a 
sufficient spread of data may not be available. 

The currrent version of the LAIRS Program assumes that no 
more than three separate meteorological profiles are to be 
input to the system, and thus cannot solve for the entire 
set of five coefficients at each altitude level. Instead, 
only the diurnal and semidiurnal amplitudes and the latitude 
gradients are calculated, the phases remaining at their de- 
fault values. No iteration is attempted. The user should 
be aware of how severely this process is limited by the 
quality and amount of available data and should use it only 
with the utmost care. 
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SECTION 5 


OPERATIONAL USE OF LAIRS 


As stated in Section 1, the primary function of the LAIRS 
Program is to produce a set of atmospheric parameters for 
each point along the Shuttle reentry trajectory. In opera- 
tional use, as soon as both the BET file, containing Shuttle 
positions as a function of time, and the meteorological data 
files are available, the LAIRS Program must output a single 
"best estimated atmosphere" covering the entire reentry 
path. This section describes the development of an optimal 
procedure for producing the best estimated atmosphere and 
makes recommendations for use of the LAIRS Program in an R&D 
environment. 

5.1 LAIRS OPTIMAL PROCEDURE 

The production of a single, smooth, continuous atmosphere 
giving the best fit to all the available data and simulta- 
neously obeying the known laws of atmospheric physics over 
the entire range of the Shuttle reentry presents a number of 
challenges. The range covered by the Shuttle during reentry 
is quite large; over 200 kilometers in altitude, 50 degrees 
in latitude, and 120 degrees in longitude. The physical 
conditions and the quantity and quality of available meteor- 
ological data vary enormously over this range. Therefore, 
no single method can yield the best possible estimates of 
all atmospheric parameters along all portions of the reentry 
trajectory. However, it is possible to prescribe an "opti- 
mal" procedure that will produce a best overall atmosphere. 
This atmosphere can then serve as a baseline for investi- 
gators who may, if they wish, use LAIRS in different modes 
to study more fully particular regions of the atmosphere 
along the Shuttle reentry path. 
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5.1.1 SELECTION OF METEOROLOGICAL DATA 


During the development phase of the LAIRS Program, it was 
expected that the meteorological data for operational use 
would be provided by NOAA in the form and format discussed 
in Section 2.1.2, i.e., up to three complete sets of verti- 
cal meteorological profiles would be provided. In actu- 
ality, the data available for the STS-1 studies consisted of 
several vertical profiles that were supplied by NOAA, NASA, 
and the Air Force, some of which contained only a partial 
set of atmospheric parameters. These vertical profiles were 
made at Barking Sands, Hawaii, and at Point Mugu, Wheeler 
Ridge, Tehachapi, and Edwards Air Force Base, California. 

All profile observations were made near the time of the 
Shuttle reentry. Since the LAIRS Program was designed to 
use a maximum of three files and since the Wheeler Ridge and 
Tehachapi profiles lacked density data, the Barking Sands, 
Wheeler Ridge, and Edwards profiles were chosen as the pri- 
mary LAIRS input (subsequent investigations with the LAIRS 
differential correction process showed that the Wheeler 
Ridge and Tehachapi data yield results that are statisti- 
cally equivalent to those obtained using the Edwards data) . 
This data is discussed in more detail in Reference 8. 

5.1.2 PRELIMINARY MODELING APPROACHES 

The preliminary runs made with the LAIRS Program using this 
data used the straight interpolation model. Predictably, 
the results were quite rough, as shown in Figure 5-1, To 
produce a smoother fit, the Jacchia-Roberts model was used 
in the upper atmosphere and the polynomial (non-DC) fitting 
model was used in the lower atmosphere, with the temperature 
values calculated from the pressure and density values via 
the gas law. 

The polynomial model requires the establishment of a refer- 
ence profile as discussed in Section 4.2.2. It was 
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Figure 5-1. Interpolated Temperature Profile for STS-1 
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discovered that the location of the reference profile is of 
great importance in determining the smoothness of the esti- 
mated atmosphere because of the uncertainties in the diurnal 
and semidiurnal coefficients and in the latitude gradients. 
For example, when the reference profile is placed at Edwards 
Air Force Base at the time of Shuttle touchdown, the esti- 
mated atmosphere is somewhat rough in the vicinity of 
90 kilometers. Moving the reference profile to Barking 
Sands produces a smooth curve in the upper part of the lower 
atmosphere, but it becomes somewhat rougher in the lower 
part. A compromise position for the reference profile may 
be used, but the fact that the parameters must first be 
translated from the meteorological file location to the ref- 
erence location, then subjected to the fitting process, and 
then translated to the Shuttle position allows the uncer- 
tainties to propagate excessively when the total times and 
distances involved are large. 

5.1.3 TWO-STEP MODELING APPROACH 

To minimize the errors introduced in the translations and to 
ensure a smooth representation of the atmosphere, a revised 
procedure was developed which uses the LAIRS Program in a 
two-step mode. First, a LAIRS run is made using the inter- 
polation model to yield unsmoothed values of the atmospheric 
parameters translated to the Shuttle trajectory. A new 
LAIRS routine outputs these parameters in the same format as 
a standard LAIRS meteorological file. In the second step, 
this new file, called USEMET, is input to a LAIRS polynomial 
model run with the latitude, diurnal, and semidiurnal cor- 
rections turned off. 

This process, in effect, takes the Shuttle trajectory as the 
reference profile location. The result is a smooth set of 
parameters that are in agreement with the interpolated and 
translated data (see Figure 5-2) . The number of points on 
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the USEMET file can be controlled by the interval field on 
the ENTREE card and can be set approximately equal to the 
original number of data points. The creation of the USEMET 
file is invoked by the keyword card CREATE, 

5,1,4 MODEL REFINEMENTS 

The parameters obtained using the above procedure still 
cannot be guaranteed to be the "best" set that can be pro- 
duced, In the analysis of the STS-1 data, it was discovered 
that changing the atmospheric segment boundaries via the 
MODBOUND keyword card could substantially improve the poly- 
nomial fit. There is no formal procedure for setting these 
boundaries; the user must inspect the interpolated data 
(Figure 5-1) and experiment. For STS-1, boundaries of 
90 kilometers, 68 kilometers, and 14 kilometers produced 
good results. 

Another variable that can be adjusted is the merging inter- 
val at each segment boundary. This is the region around the 
boundary where the polynomial coefficients for the segments 
on both sides of the boundary are computed. For example, if 
the bottom segment boundary is 25 kilometers and the merging 
region is set at 5 kilometers, the meteorological data up to 
30 kilometers is used in determining the polynomial coeffi- 
cients for the bottom segment, and the data down to 20 kilo- 
meters is used in determining the coefficients for the 
middle segment. In the parameter calculations, however, the 
middle segment coefficients are used only down to 25 kilo- 
meters and the lower segment coefficients are used below 
25 kilometers. This scheme was devised to reduce possible 
discontinuities in the derivatives of the parameters at the 
boundaries in a differential correction (DC) run (in a 
non-DC run, this approach is used to reduce the possibility 
of discontinuities in the parameters themselves) . Thus, 
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making the merging region too small can have undesirable 
effects. Figure 5-3 shows the results of a differential 
correction run with a merging parameter of 1 kilometer. The 
curve mimics the data reasonably well in most places and it 
is continuous, but it looks decidedly nonphysical. Exper- 
ience with STS-1 suggests a merging factor of at least 3 or 
4 kilometers (see Figure 5-2) . The merging parameter is set 
on the MODBOUND keyword card. 

In a differential correction run, several other parameters 
may be changed to alter the fit. These include the editing 
criteria and, when no measurement uncertainties are avail- 
able, the relative a priori weights of the temperature, 
pressure, and density. These are controlled by the DCEDIT 
keyword card. 

5.2 R&D USE OF LAIRS 

LAIRS may be used in a variety of ways to improve the esti- 
mate of atmospheric parameters in any particular region of 
the atmosphere. If one is interested in only a small por- 
tion of the reentry trajectory, e.g., around 90 kilometers, 
the two-step approach described above becomes unnecessary, 
and the user may simply place the reference profile at the 
time and place of interest. Segment boundaries, editing 
criteria, weighting factors, etc., can be adjusted to give 
optimum results in the region under study as long as the fit 
in other regions is not critical. 

With regard to determining the reliability of the estimated 
parameters, the following procedure may prove useful. 

First, a LAIRS default model run for the day of reentry will 
give the user an approximate knowledge of the values to be 
expected from the actual data. When the actual data becomes 
available, the "best estimated atmosphere" can be produced 
via the optimal procedure. This atmosphere can then be com- 
pared with the original data (bearing in mind that diurnal 
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and latitudinal corrections have been made and that, in a 
differential correction run, all three parameters may be 
used to produce the values of any single parameter) . The 
best estimated atmosphere should also be judged against the 
intermediate interpolated file used to create it. Finally, 
if onboard pressure measurements are available from the 
Shuttle, they can be compared with onboard pressure measure- 
ments estimated by LAIRS. This provides an independent 
check on the reliability of the calculated parameters. 
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APPENDIX A 


COMMON BLOCK DESCRIPTIONS 


This appendix contains descriptions of all COMMON 
used in the LAIRS Program. The COMMON blocks are 
in alphabetical order. 


blocks 

listed 
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BLOCK DATA CALCOH 


COMMON BLOCK CALCOM 
PURPOSE t 

/CALCOM/ STORES ONE BLOCK OF SHUTTLE TRAJECTORY POINTS OR LAP 
POINTS AND THE ATMOSPHERIC PARAMETERS CALCULATED AT THOSE 
POINTS. ONE BLOCK « 10 POINTS. 


REFERENCED BY; LAIRSf PARAMS»OUTPRFrLAPf ENTRD 


COMMON BLOCK PARAMETERS: 


NAME 

TYPE 

DIMENSION 

TRJ 

R 

5rl0 

TPDW 

R 

7»10 


UTRJ 

R 

Sr 10 

UTPDW 

R 

7rl0 

UTMOL 

R 

10 

8CHT 

R 

10 

V(10> 

R 

10 

ALPHA 

R 

10 

RMC 

R 

10 

PNC 

R 

10 

DEV 

R 

10 


DESCRIPTION 
10 TRAJECTORY POINTS: 

ALTITUDE fLAT I TUBE » LONGITUDE » GMT » 
LOCAL SOLAR TIME 

10 PARAMETER SETS CORRESPONDI6 TO 
THE TRAJECTORY POINTS. EACH SET 

contains: 

TEMPERATURE > PRESSURE * DENSITY f 
E-U UIND Uf N-S WIND Ur WIND 
WIND SPEEDr AND UIND DIRECTION 
UNCERTAINTIES IN THE TRAJECTORY 
COMPONENTS 

UNCERTAINTIES IN THE CALCULATED 
PARAMETERS 

MEAN MOLECULAR WEIGHTS 
PRESSURE SCALE HEIGHT (KM) 

SHUTTLE VELOCITIES (M/SEC) 

SHUTTLE ANGLES OF ATTACK (DEG) 
MACH NUMBERS 

CALCULATED ONBOARD MEASUREMENTS 
PERCENT DEVIATIONS FROM GAS LAW 


TOTAL size: 310 WORDS 
programmer: D. E. BOLANDr CSC 


Block Data CALCOM 



BLOCK DATA COFCOM 

purpose: 

/COFCOM/ STORES POLYNOMIAL COEFFICIENTS AND MODEL STANDARD 
DEVIATIONS <UNCERTAINTIES> FOR TEMPERATURE, LOG 
PRESSURE, AND LOG DENSITY FOR EACH OF THE 3 
ATMOSPHERIC SEGMENTS, 

REFERENCED BY*, LAIRS , POLE IT , POLCAL , OUTSUM » POLY 
COMMON BLOCK PARAMETERS: 


NAME 

type 

dimension 

DESCRIPTION 

CT 

R 

3 , A 

4 TEMPERATURE COEFFICIENTS FOR 
EACH ATMOSPHERIC SEGMENT 

CP 

R 

3,4 

4 PRESSURE COEFFICIENTS FOR 
EACH ATMOSPHERIC SEGMENT 

CD 

R 

3,4 

4 DENSITY COEFFICIENTS FOR 
EACH ATMOSPHERIC SEGMENT 

UT 

R 

3 

TEMPERATURE STANDARD DEVIATION 
FOR EACH SEGMENT 

UP 

R 

3 

PRESSURE STANDARD DEVIATION 
FOR EACH SEGMENT 

UD 

R 

3 

DENSITY STANDARD DEVIATION 
FOR EACH SEGMENT 

COVAT 

R 

3,5,5 

COVARIANCE MATRIX FOR TEMPERATURE FOR 
EACH SEGMENT 

COVAP 

R 

3,5,5 

COVARIANCE MATRIX FOR PRESSURE FOR 
EACH SEGMENT (CURRENTLY DUMMY) 

COVAD 

R 

3,5,5 

COVARIANCE MATRIX FOR DENSITY FOR 
EACH SEGMENT (CURRENTLY DUMMY) 

TOTAL size: 

270 WORDS 


programmer: 

D t E! « 

BOLAND, CSC 



Block Data COFCOM 
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BLOCK DATA CONCOM 


purpose: 

/CONCOM/ STORES CONSTANTS NEEDED BY VARIOUS ROUTINES 
REFERENCED BY: 6WIND » DIGRAD f TRNDSD » DSDADJ f WCNVRT i HI ALT f JACROB » JR 


COMMON BLOCK PARAMETERS: 


NAME 

TYPE 

DIMENSION 

PI 

R 

1 

TWOPI 

R 

1 

DTR 

R 

1 

RTD 

R 

1 

HTD 

R 

1 

ANG 

R 

1 

R 

R 

1 

DUM 

R 

13 


DESCRIPTION 

PI 

TWO TIMES PI 

DEGREES TO RADIANS CONVERSION 
FACTOR 

RADIANS TO DEGREES CONVERSION 
FACTOR 

HOURS TO DEGREES CONVERSION FACTOR 
ANGULAR VELOCITY OF EARTH 
RADIUS OF EARTH 
SPARES 


TOTAL size: 20 WORDS 
PROGRAMMED BY: D. E* BOLAND» CSC 
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BLOCK DATA DCCOH 
COHMON BLOCK NAHEJ /DCCOH/ 

purpose: 


/DCCOH/ 

KEYWORD 

COHHOH BLOCK 

STORES DC PROCESS 
CARDS . 

PARAHETERS: 

CONTROL FLAGS READ IN THROUGH DC 

NAHE 

TYPE 

DIHENSION 

DESCRIPTION 

HLDATA<N»L> I 

(3»3) 

HETE0R0L06ICAL DATA SELECTION FLAGS: 
=0f DO NOT INCLUDE L TYPE HET. DATA 
IN DC FOR SE6HENT N 
=lf INCLUDE 
L=l» TEHPERATURE 
L=2» PRESSURE 
L=3f DENSITY 

KSSEG(N) 

I 

3 

CONTINUITY flags: 

“Ot CONTINUITY IS IHPOSED AT THE TOP 
BOUNDARY OF THE NTH SEGHENT 
=1» CONTINUITY IS NOT IHPOSED 

ISQH 

I 

1 

DC SEGHENTATION FLAG: 

»0* NORHAL DC (CONTINUITY IHPOSED AT 
ALL BOUNDARIES) 

=lr SEGHENTWISE INDEPENDENT DC 

IDCSEO(N) 

I 

3 

DC FLAG FOR EACH SEGHENT: 
aOf NO DC FOR THE NTH SEGHENT 
=tf DC FOR THE NTH SEGHENT 

WTDCT 

R 

3 

A PRIORI WEIGHTING FOR TEHPERATURE 
FOR EACH SEGHENT (5Z UNCERTAINTY) 

WTDCP 

R 

3 

A PRIORI WEIGHTING FOR PRESSURE 
FOR EACH SEGHENT (5X UNCERTAINTY) 

UTDCD 

R 

3 

A PRIORI WEIGHTING FOR DENSITY 
FOR EACH SEGHENT (5X UNCERTAINTY) 

ALFA 

TOTAL size: 

R 5 

30 WORDS 

INTEGRANDS FOR PARTIAL DERIAVATIVE3 
OF PRESSURE AND DENSITY 

prograhher: 

T. LEE 

r CSC 

JULY 1981 


Block Data DCCOM 
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BLOCK DATA FILCOH 


purpose: 

/FILCOH/ STORES THE LOGICAL UNIT NUMBERS OF VARIOUS FILES* 
AND STORES INFORMATION NEEDED BY THE OUTPUT ROUTINES* 
FILCOH IS ACCESSED BY VIRTUALLY ALL LAIRS ROUTINES* 


REFERENCED BY: 


COMMON BLOCK PARAMETERS: 


NAME 

TYPE 

DIMENSION 

NPRINT 

I 

1 

NREAD 

1 

1 

NUSE 

I 

1 

NPRM 

I 

1 

NURK 

I 

1 

NENT 

I 

1 

NKP 

I 

1 

NOIDS 

I 

1 

NUSMET 

I 

1 

NDUM 

I 

2 

MWRKl 

I 

1 

MHRK2 

I 

i 

MURK3 

I 

1 

NDUM2 

I 

8 

NITER 

I 

1 

LINES 

I 

1 

IPAGE 

TOTAL size: 

1 

25 WORDS 

1 


programmer: d* e* boland* esc 


DESCRIPTION 
PRINTED OUTPUT LUN 
CARD INPUT LUN 
USER FILE LUN 
PERMANENT FILE LUN 
WORKING FILE LUN 
ENTREE FILE LUN 
KP FILE LUN 

COMMON BLOCK VALUES LUN 
CREATED METEOROLOGICAL FILE LUN 
DUMMIES 

METEOROLOGICAL FILE LUN 
METEOROLOGICAL FILE LUN 
METEOROLOGICAL FILE LUN 
DUMMIES 

ITERATION COUNTER FOR DC OUTPUT 
LINE COUNTER 
PAGE COUNTER 


Block Data FILCOM 



BLOCK DATA INTCOM 

purpose: 

/INTCOM/ STORES VALUES RETRIEVED FROM THE WORKING FILE BY DFLINT» 

REFERENCED BY: LAIRS r DFL INT » METINT r POLY » SETREF f WIND » POLADJ » 

METREFf TRNADJ 

COMMON BLOCK PARAMETERS: 


NAME 

TYPE 

DIMENSION 

DESCRIPTION 

IPT 

1 

1 

INTERPOLATION POINTER (LAST POINT 
REFERENCED) 

IPTFLG 

I 

1 

FLAG FOR FIRST CALL TO DFLINT 

ALT 

R 

6 

BLOCK OF RETRIEVED ALTITUDES 

TPDWAR 

R 

6f 5 

BLOCK OF RETRIEVED PARAMETERS: 
TrP»D»U»V 

ARR 

R 

6»27 

ALL OTHER VALUES RETRIEVED FROM 
WORKING FILE (SAME ORDER AS ON FILE) 

TOTAL size: 

260 

WORDS 


programmer: 

R. 

A. KUSESKIf 

CSC. 


Block Data INTCOM 
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BLOCK DATA JRCOH 


COMMON/JRCOM 

/ PrtJ 

rHCJ 

fTINF 

fTX 

*DNHE * 

» 


CM(6) 

»CF 

fDD(5f7) 

*ADT(6) 

*ZJX * 

t 


AUG 

fRCM 

fTO 

*RC 

*GLO * 

« 


ZJO 

fFKL 

»TZ 

*TCIL 

*RL1 > 



DUHl 

»RL2 

rDUM2 

*XLPS 

*YLPS * 

t 


UC(2) 

fUC(2) 

fXCDI 

fUCDI 

»FLC4 * 

» 


CC(S> 

»ZD(7) 

»SD(5) 

*CMZ 

* 

t 


RHOZ 

f AD<6) 

f BD<6) 

*AC<7) 


« 


fCFL(5) 

»H6T 

fDPJDT 

*FITPRM<4) 


HASSES 

OF THE 

ATMOSPHERIC 

CONSTITUENTS IN GM/MOLE. 


DATA 

CM 

/28. 013400 » 39. 94800 

*4.002600 

*31.998600 

*15.999400* 

t 


1.0079700/ 





CORRECTION FACTOR IN KM 

DATA CF /U825.00 / 

POWER SERIES COEFFICIENTS FOR COMPUTING RHOS AT INFLECTION POINT 
<NO LONGER USED. SEE SUBROUTINE MODADJ DESCRIPTION.) 


DATA DD 

/0.1093155E2 

f0.8049405El 



.7646886E1 

» .9924237E1 



.1097083E2 

* .1186783E-2 



*2382822E~2 

*-.4383486E-3 



.1600311E-2 

* .6118742E-4 



~.1677341E-5 

*-.3391366E-5 



.4694319E-6 

»-.2274761E-5 



-.1165003E-6 

* .1420228E-8 



.2909714E-8 

*-.2B94886E-9 



.1938454E-8 

* .9239354E-10 



“♦7139785E-12 

*-.1481702E-ll 



.9451989E-13 

*-.9782183E-12 



-.3490739E-13 

* .1969715E-15 



.4127600E-15 

f-.1225838E-16 



.2698450E-15 

* .5116298E-17 



-.2296182E-19 

»“.4837461E“19 



~.3131808E-19 

*0.00 

/ 


THERMAL DIFFUSION COEFFICIENTS 


DATA ADT /2»0.00 f-.3800 »3»0.00 / 

INFLECTION POINT HEIGHT IN KM 

DATA ZJX /125.00/ 


Block Data JRCOM (1 of 3) 
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AVOGADRO'S NUMBER IN NUMBER / MOLE 


DATA AVG /6.02257E23/ 

AVERAGE EARTH RADIUS IN KM 
DATA ROM /6356. 76600/ 

TEMPERATURE AT MINIMUM HEIGHT IN DEGREES KELVIN 
DATA TO /183.00/ 

UNIVERSAL GAS CONSTANT IN JOULES/DEGREE KELVIN/MOLE 
DATA RC /8«3143200/ 

MEAN SURFACE GRAVITY IN M/SEC*!H2 
DATA GLO /9. 80665000/ 

MINIMUM HEIGHT IN KM 
DATA ZJO /90.00/ 

INITIAL APPROXIMATIONS TO RLl » RL2 i XLPS » YLPS 

DATA RLl »RL2rXLPSf YLPS /166.00 f61.00 i99.00 »28»00 / 

POWER SERIES COEFFICIENTS FOR TZ 

DATA CC /-,8928437500E8 f.35424E7 r-*526875E5 » 

» t3405E3 f-.800 / 

POWER SERIES COEFFICIENTS FOR RH0(100)/M0 

DATA ZD /.1985549E-10 » - « 1833490E-14 i 

* .1711735E-17 1021474E-20 i 

» »3727894E-24 f 77341 lOE-28 i 

» .7026942E-32 / 

COEFFICIENTS FOR COMPUTING NUMBER DENSITY AT Z^IOOKM 

DATA SD /♦7811000 » ,93432E-2f .61471E-5 f ♦ 16177800r .09554400/ 

MEAN MOLECULAR MASS AT SEA LEVEL IN 6M/M0LE 

DATA CMZ /28. 96000 / 
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ATMOSPHERIC DENSITY AT Z=90 KM IN GM/CM**3 


DATA RHOZ /3.46E-9/ 

CONSTANT COEFFICIENTS FOR THE B(N) SERIES 


DATA AD 

* 


/.3144902516672729E10 

12377488548329 17E9» * 1816141096520398E7 
11 40331 079489267E5» .243649861 21 05595E2 
,8957502869707995E~2/ 


LINEAR COEFFICIENTS FOR THE B(N) SERIES 


DATA BD 

* 

t 


/-♦5286448217910969E8» 

1663250847336828E5»-* 1308252378125E1 
3*0.00 / 


MOLECULAR POWER SERIES COEFFICIENTS 


DATA AC 

* 

« 

» 


/“.435093363387E6 

-.765334661080E3 

-.8958790995E-1 

-.697444E-6 


» .282755646391E5 
t .110433875450E2 
» .38737586E-3 
/ 


COEFFICIENTS FOR COMPUTING CF 


DATA CFL /.1031445E5 

* .1579202E-2 

* .2462708E-9 
..END 


f .2341230E1 
f-. 1252487E-5 
/ 


Block Data JRCOM (3 of 3) 



BLOCK DATA METCOM 

purpose: 

/METCOM/ STORES VALUES FROM METEOROLOGICAL PROFILES USED FOR 
METEOROLOGICAL INTERPOLATION. 

REFERENCED BYt LAIRS » MET INT » METINl i TRNADJ » SETREF 


COMMON BLOCK 
NAME 

parameters: 

TYPE DIMENSION 

DESCRIPTION 

IPT 

I 

1 

POINTER USED IN LIBRARY ROUTINE 

METFLG 

1 

1 

lUNl 

FLAG FOR FIRST CALL TO METINl 

11 

I 

1 

POINTER TO CLOSEST ALTITUDE IN 

ALTM 

R 

Zt6 

TRAJECTORY ARRAY 
ARRAY OF ALT, FOR EACH PROFILE 

METARR 

R 

3f 6r 10 

ARRAY OF T»PfD AND WIND SPEED AND 

FSTPT 

R 

3 

DIRECTION FOR EACH PROFILE 
LOWEST ALT. FOR EACH PROFILE 

RLAT 

R 

3 

LATITUDE OF EACH PROFILE 

RLONG 

R 

3 

LONGITUDE OF EACH PROFILE 

RLSTM 

R 

3 

LOCAL SOLAR TIME FOR EACH PROFILE 


TOTAL size: 240 WORDS 

programmer: r.a. kuseski» computer sciences corp. 


Block Data METCOM 
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BLOCK DATA SOLAND 
PURPOSE! 

/SOLAND/ STORES DIFFERENTIAL CORRECTION CONTROL PARAMETERS 
AND STATISTICS 

REFERENCED BY! NLSPOL » SOLVE » GRAND » DPAR » TPAR » PPAR 
COMMON BLOCK PARAMETERS! 


NAME 

TYPE 

DIMENSION 

DESCRIPTION 

IMAX 

I 

1 

MAX. NUMBER OF INTEGRALS FOR SIMP 

KS4 

I 

1 

FLAG FOR BOUNDARY TEMPERATURE SOLVE-FOR 

KS5 

I 

1 

FLAG FOR BOUNDARY PRESSURE SOLVE-FOR 

NPOL 

I 

1 

ORDER OF TEMP. POLYNOMIAL 

MAXSOL 

I 

1 

MAX. NUMBER OF SOLVE-FOR PARAMETERS 

NSOL 

1 

1 

ACTUAL NUMBER OF SOLVE-FOR PARAMETERS 

KEDTOT 

I 

1 

TOTAL NUMBER OF MET. DATA USED 

K0NVR6 

I 

1 

CONVERGENCE FLAG 

IDUM 

I 

2 

SPARES 

ZREF 

R 

1 

REFERENCE ALTITUDE 

EPSLON 

R 

1 

CONVERGENCE TOLERANCE CRITERIA 

SRSUM 

R 

1 

SUM OF SQUARED WEIGHTED RESIDUALS 

CURRMS 

R 

1 

CURRENT WEIGHTED RMS 

OLDRMS 

R 

1 

PREVIOUS WEIGHTED RMS 

RMSP 

R 

1 

PREDICTED RMS 

A 

R 

10 

SOLVE-FOR PARAMETER ARRAY 

6 

R 

10 

RIGHT HAND SIDE OF NORMAL EQUATION 

SH 

R 

55 

LINEAR ARRAY OF UPPER TRIANGULAR PART 
OF Nuh-nwL MATRIX 

RDUM 

R 

9 

SPARES 

PROGRAMMER! 

T. LEE 

CSC MARCH 

1981 


Block Data SOLAND 
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BLOCK DATA TKPTC 
COHMON BLOCK NAME: /TKPTC/ 

purpose; 

/TKPTC/ INITIALIZE VALUES FOR THE KP FILE 


REFERENCED BYt L A I RS » J ACCWF , JACROB 


COMMON 

BLOCK 

parameters; 




NAME 

TYPE 

dimension 

DESCRIPTION 



IT 

1 

1 

TIME OF FIRST DAY 
IN COMMON 

OF TC 

DATA 

KP 

1 

21f8 

MAGNETIC ACTIVITY 

3-HOUR 

INDICES 

TC 

R 

20 

EXOSPHERIC TEMPERATURE 



TOTAL size: 189 WORDS 

PROGRAMMER ANALYST; A. K. KAPOORr COMPUTER SCIENCES CORPORATION 


Block Data TKPTC 
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BLOCK DATA USECQM 


COMMON BLOCK NAME: /USECOM/ 


purpose: 

/USECOM/ STORES THE DATA INPUT BY THE USER ON KEYWORD CARDS» 
USECOM IS ACCESSED BY VIRTUALLY ALL LAIRS ROUTINES* 


COMMON BLOCK 

parameters: 


NAME 

TYPE 

DIMENSION 

DESCRIPTION 

IPRFL 

I 

i 

TYPE OF profile: 
=1» ENTREE 
-2r LAP 

MODL 

I 

1 

TYPE OF LOWER ATMOSPHERIC MODEL: 
=lf DEFAULT INTERPOLATION 
=2» METEOROLOGICAL INTERPOLATION 
=3r DEFAULT POLYNOMIAL 
-At ADJUSTED POLYNOMIAL 

lUL 

I 

1 

TYPE OF L«A. WIND MODEL: 

=Or NO WIND 

=lf EMPIRICAL (INTERPOLATED) WIND 
=2f ANALYTIC WIND 

ILATL 

1 

1 

LATITUDE VARIATIONS IN L*A.: 
=Of DO NOT INCLUDE 
=lf INCLUDE 

IDSDL 

I 

1 

DIURNAL/SEMIDIURNAL VARIATIONS! L.A 
=Of DO NOT INCLUDE 
= lf INCLUDE 

MODU 

I 

1 

TYPE OF UPPER ATMOSHPHERIC MODEL: 
=l! DEFAULT INTERPOLATION 
=2» METEOROLOGOCAL INTERPOLATION 
=3r DEFAULT JACCHI A-ROBERTS 

lUU 

I 

1 

TYPE OF U»A* WIND MODEL 
=Or NO WIND 

=1t empirical (INTERPOLATED) WIND 

ILATU 

1 

1 

LATITUDE VARIATIONS IN U*A* 

=Of do not INCLUDE! EXCEPT IN JR 
=l! INCLUDE 

IDSDU 

I 

i 

DIURNAL/SEMIDIURNAL VARIATIONS! U«A 
=0! DO NOT INCLUDE! EXCEPT IN JR 
=l! INCLUDE 

IDRV 

I 

1 

DERIVE PRESSURE AND DENSITY FROM 
TEMPERATURE VIA PHYSICAL LAWS 
(OVERRIDES ALL OTHER SELECTIONS) 

INTP 

1 

1 

ALTITUDE INTERPOLATOR TYPE: 
=l! FIRST ORDER LAGRANGE 
=2! SECOND ORDER LAGRANGE 


Block Data USECOM (1 of 3) 
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I DBG 

I 

1 

PRINT debug; 
= 0f NO 
=1» YES 

INT 

I 

1 

HORIZONTAL INTERPOLATOR TYPEt 
=Of UNIVARIATE 
=lr FIRST ORDER BIVARIATE 
=2» SECOND ORDER BIVARIATE 

IDNUM 

1 

1 

ENTREE-LAP ID NUMBER 

I6AS 

I 

1 

GAS LAW CALCULATION OVERRIDE 

NSAMP 

I 

1 

ENTREE FILE SAMPLING INTERVAL < POIN' 

ITER 

I 

1 

MAXIMUM NUMBER OF ITERATIONS 

IWT 

I 

1 

INCLUDE weighting; 

= 0» NO 
=1» YES 

IDSDA 

I 

1 

ADJUST D/SD coefficients; 
= Or NO 
=1» YES 

ILATA 

I 

1 

ADJUST LATITUDE GRADIENTS; 
= Of NO 
=1» YES 

I DC 

I 

1 

PERFORM DIFFERENTIAL CORRECTIONS; 
= NO 

-1? YES 

ITDAT 

I 

1 

INCLUDE TEMPERATURE DATA IN DC ; 
= Of NO 
=lr YES 

IPDAT 

I 

1 

INCLUDE PRESSURE DATA IN DC; 
= Or NO 
=1» YES 

IDDAT 

I 

1 

INCLUDE DENSITY DATA IN DC; 
= 0» NO 
=1» YES 

IDMY 

I 

1 

DUMMY 

NFIL 

I 

1 

NUMBER OF METEOROLOGICAL FILES 

KWRITE 

I 

1 

FULL OUTPUT PRINT FLAG 

ICRT 

I 

1 

FLAG FOR CREATION OF NEW 
METEOROLOGICAL FILE FROM USER FILE 

IDUM 

I 

2 

SPARES 

YMD 

R 

1 

YEAR» MONTH* DAY OF LAP OR ENTREE 
TRAJECTORY 

HhS 

R 

1 

HOURS* MINUTES* SECONDS OF FIRST 
POINT ON ENTREE TRAJECTORY 

THETA 

R 

1 

LATITUDE OF LAP 

TAU 

R 

1 

LOCAL SOLAR TIME OF LAP 

Z1 

R 

1 

LOWER PROFILE ALTITUDE 

Z2 

R 

1 

UPPER PROFILE ALTITUDE 


Block Data USECOM (2 of 3) 



ZDUM 

R 

6 

DUMMIES 

BOUND 

R 

3 

BOUNDARIES BETWEEN 4 SEGMENTS: 
90KM» 65KM» 25KM 

HMSSEC 

R 

1 

ENTREE FILE HMS START TIME IN 
DECIMAL SECONDS 

ZREF 

R 

1 

REFERENCE ALTITUDE (KM) FOR 
DERIVATION OF P AND D 

PREF 

R 

1 

PRESSURE AT ZREF 

STEP 

R 

1 

INTEGRATION STEPSIZE (KM) 

CONU 

R 

1 

CONVERGENCE CRITERION FOR DC 

RHERGE 

R 

1 

REGION FOR POLYNOMIAL MERGING 
ON EACH SIDE OF SEGMENT 
BOUNDARIES (KM) 

REFLAT 

R 

1 

MODEL REFERENCE LATITUDE 

REFTAU 

R 

1 

MODEL REFERENCE LOCAL SOLAR TIME 

RMGMT 

R 

3 

GMT OF METEOROLOGICAL FILES 

DUH 

TOTAL size: 

R 

70 WORDS 

16 

SPARES 


Block Data USECOM (3 of 3) 
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APPENDIX B 


FILE DESCRIPTIONS 


This appendix specifies the formats of the files that are 
created and used by the LAIRS Program. Detailed descriptions 
of the working file, the user file, and the permanent file 
are presented in this appendix. 



WORKING FILE 


The working file is an indexed sequential file that stores 
the default atmospheric profile for the day of the Shuttle 
reentry or the local atmospheric profile. Each record is 
38 words in length and there is one record per altitude 
level. The number of records depends on the altitude in- 
terval at which the default data has been tabulated,- tabu- 
lation at each 5 kilometers up to 110 kilometers results 
in 23 records. Up to 50 records are allowed. The record 
format is the same as that of the permanent file. 
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USER FILE 


The user file is a sequential file containing the calcu- 
lated atmospheric parameters for each point on a local at- 
mospheric profile or Shuttle trajectory. There is a title 
record describing the file, a header record, and one rec- 
ord for each altitude level at which the parameters were 
calculated. Each record contains 29 words. 


The format 
Variable 

of the 
Type 

title record is 
1 Dimension 

as follows: 

Description 

TITLE 

A 

8 

Description of the file 

The format 

of the 

header record is 

as follows: 

Variable 

Type 

Dimension 

Description 

ID 

I 

1 

ENTREE-LAIRS ID Niomber 

DATE 

R 

1 

YYMMDD. of LAP or 
ENTREE Trajectory 

TIMEO 

R 

1 

HHMMSS. SS (GMT) of first 
point on ENTREE file 

NSPR 

R 

1 

Spare 

The foirmat 

of the 

data record is as follows: 

Variable 

Type 

Dimension 

Description 

Z 

R 

1 

Altitude (km) 

THETA 

R 

1 

Latitude (decimal de- 
grees) 

PHI 

R 

1 

Longitude (decimal de- 
grees) 

GMT 

R 

1 

Time (HHMMSS. SSSS) 

TAU 

R 

1 

Local solar time 


(HHMMSS .SSSS) 
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Variable 

Type 

Dimension 

Description 

T 

R 

1 

Temperature ( K) 

P 

R 

1 

2 

Pressure (n/m ) 

D 

R 

1 

3 

Density (kg/m ) 

U 

R 

1 

E-W wind component 
(m/sec) 

V 

R 

1 

N-S wind component (m/sec) 

WS 

R 

1 

Wind speed (m/sec) 

WD 

R 

1 

Wind direction (degrees 
from N, positive clockwise) 

TU 

R 

1 

Temperature uncertainty 
{ K) 

PU 

R 

1 

Pressure uncertainty 
(n/m^) 

DU 

R 

1 

Density uncertainty 
(kg/ta^) 

UU 

R 

1 

U uncertainty (m/sec) 

VU 

R 

1 

V uncertainty (m/sec) 

WSU 

R 

1 

WS uncertainty (m/sec) 

V7DU 

R 

1 

WD uncertainty (degrees) 

ZU 

R 

1 

Altitude uncertainty (km) 

THETAU 

R 

1 

Latitude uncertainty (degrees) 

PHIU 

R 

1 

Longitude uncertainty (degrees) 

GMTU 

R 

1 

Time uncertainty (blank) 

LSTU 

R 

1 

Local solar time uncertainty 
(blank) 
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MWT 

R 

1 

Mean molecular weight 
(g/mole) 

SCHT 

R 

1 

Pressure scale height (km) 

M 

R 

1 

Mach number 

PN 

R 

1 

Estimated onboard pressure 

2 

measurement (n/m ) 

PCDEV 

R 

1 

Percent gas law deviation 

PERMANENT FILE 



The permanent 

file is 

a sequential 

file that stores the 


default atmospheric profile for each month of the year. 

There are 12 sets of records, one record per altitude level, 
containing atmospheric parameters for each month. There 
are 38 words per record, including spares. To uniquely 
distinguish each record by month and by altitude, the record 
identification key is a floating-point number, the whole 
part of which is the number of the month and the decimal 
part of which is the altitude, e.g., 

1.025 = record for altitude of 25 kilometers in January 
10.105 = record for altitude of 105 kilometers in October 


The format of the header record 
Variable Type Dimension 


KEY 


R 


1 


NREC I 


1 


NJAN I 


1 


is as follows : 

Description 

Header identification 
keyword 

Total number of records 
in the file 

Number of records for 
January 
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Variable 

Type 

Dimension 

Description 

T 

R 

1 

Temperature ( K) 

P 

R 

1 

2 

Pressure (n/m ) 

D 

R 

1 

3 

Density (kg/m ) 

U 

R 

1 

E-W wind component 
(m/sec, positive east) 

V 

R 

1 

N-S wind component 
(m/sec, positive north) 

DTLAT 

R 

1 

Latitude temperature gra- 
dient { K/degree latitude) 

DPLAT 

R 

1 

Latitude pressure gradient 
(n/m /degree latitude) 

DDLAT 

R 

1 

Latitude density gradient 

3 

(kg/m /degree latitude) 

DULAT* 

R 

1 

Latitude U gradient 
(m/sec/degree latitude) 

DVLAT* 

R 

1 

Latitude V gradient 
(m/sec/degree latitude) 

DAT 

R 

1 

Diurnal temperature am- 
plitude ( K) 

DPT 

R 

1 

Diurnal temperature phase 
(hours) 

SDAT 

R 

1 

Semidiurnal temperature 
amplitude ( K) 


•k 

Latitude and diurnal/semidiurnal coefficients for winds 
(U and V) are not presently included on the file. The value 
"-9999.” fills these locations. Below 25 kilometers, lati- 
tude and diurnal/semidiurnal coefficients are not presently 
included for any parameters. The value "-9999." fills 
these locations . 
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Variable 

Type 

Dimension 

Description 

SOFT 

R 

1 

Semidiurnal temperature 
phase (hours) 

DAP 

R 

1 

2 

Diurnal P amplitude (n/m ) 

DPP 

R 

1 

Diurnal P phase (hours) 

SDAP 

R 

1 

Semidiurnal P amplitude 
(n/m^) 

SDPP 

R 

1 

Semidiurnal P phase (hours) 

DAD 

R 

1 

3 

Diurnal D amplitude (kg/m ) 

DPD 

R 

1 

Diurnal D phase (hours) 

SDAD 

R 

1 

Semidiurnal D amplitude 
(kg/m^) 

SDPD 

R 

1 

Semidiurnal D phase (hours) 

DAU* 

R 

1 

Diurnal U amplitude (m/sec) 

DPU* 

R 

1 

Diurnal U phase (hours) 

SDAU* 

R 

1 

Semidiurnal U amplitude 
(m/sec) 

SDPU* 

R 

1 

Semidiurnal U phase (hours) 

DAV* 

R 

1 

Diurnal V amplitude (m/sec) 

DPV* 

R 

1 

Diurnal V phase (hours) 

SDAV* 

R 

1 

Semidiurnal V amplitude 
(m/sec) 

SDPV* 

R 

1 

Semidiurnal V phase (hours) 

SPARES 

R 

5 

Spares 

* 

See footnote 

on page 
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METEOROLOGICAL FILES (METl , MET2 , MET3) 

The meteorological files are sequential files containing 
observed atmospheric parameters versus altitude. As there 
will be a maximum of three meteorological profiles obtained 
per Shuttle flight, only three meteorological data files 
have been created, METl, MET2, and MET3. 

Each file has a header record that contains the following 
information; total number of points for all the profiles, 
the number of profiles, the number of points for this pro- 
file, the number of points between 95 kilometers and 65 
kilometers, the number of points between 65 kilometers and 
25 kilometers, the number of points between 25 kilometers 
and 0 kilometers, the lowest altitude, the highest altitude, 
the latitude corresponding to the first point, and the 
longitude corresponding to the first point. There is one 
data record for each altitude level at which the atmospheric 
parameters are observed. All the parameters are in mks units 
and the altitude is in kilometers. The METl, MET2, and 
MET3 files start at the highest altitude and go to their 
lowest altitude. The program PREMET sets up these files. 

File METl is written to unit 21, MET2 is written to unit 22, 
and MET3 is written to unit 23. 

The input meteorological tape has to be of a certain format; 
otherwise the routine PREMET may have to be modified. 
Currently, the input tape has to have a header record that 
contains ten card images of 80 characters each (20A4). . The 
data should then follow in units described in Reference 1. 

PREMET is designed to be executed as a preprocessor routine 
that is independent of the LAIRS Program. PREMET could be 
modified to execute a meteorological tape having fewer or 
more header records but the output from PREMET may not be 
altered without altering the LAIRS Program. 
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The format of the header record for the meteorological files 
is as follows: 


Variable 

IMET* 

NPRO* 

* 

NUMPTS 

ININTY* 

IMIDLE * 

ILOW* 

RFRST 

RLAST 

RLATl 

RLONl 


Type 


I 


I 


I 

I 


I 


I 


R 

R 

R 


R 


Dimension 


Description 


1 Total number of data 

points for n meteoro-' 
logical profiles; 
n = 1, 2, or 3. 

1 Total number of profiles 

n for this data; 
n = 1 , 2 , or 3 , 

1 Total number of points 

for this profile, 

1 Number of points between 

95 and 65 kilometers for 

this profile, 

1 Number of points between 

65 and 25 kilometers for 

this profile, 

1 Number of points between 

25 and 0 kilometers for 
this profile. 

1 The lowest altitude for 

this profile. 

1 The highest altitude for 

this profile, 

1 Latitude corresponding 

to point RFRST and used 
for this profile. 

1 Longitude corresponding 

to point RFRST and used 
for this profile. 


* 

This is not essential information. 


B-10 



The format of the data record is as follows: 


Variable 

Type 

Dimension 

Description 

RLAT 

R 

1 

Latitude (degrees north) 

RLON 

R 

1 

Longitude (degrees east) 

ALT 

R 

1 

Altitude (km) 

U 

R 

1 

E-W wind component 
(m/sec, positive east) 

V 

R 

1 

N-S wind component 
(m/sec, positive north). 

TE 

R 

1 

Ambient temperature ( K) 

PR 

R 

1 

2 

Ambient pressure (n/m ) 

D 

R 

1 

Ambient density (kg/m^) 

TEU 

R 

1 

Ambient temperature 
systematic uncertainty 
( K) 

PRU 

R 

1 

Ambient pressure 
systematic uncertainty 
(n/m2) 

DU 

R 

1 

Ambient density 
systematic uncertainty 
(kg/m^) 

UtJ 

R 

1 

U uncertainty 
systematic uncertainty 
(m/sec) 

VU 

R 

1 

V uncertainty 
systematic uncertainty 
(degrees) 
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A sample job setup that can be executed to generate the 
meteorological files is shown below. 


/JOB 

LAIRS rT200» CM 120000. 1251JMP 

USER,XXXXXXX . 

CHARGE, XXXXXX^ LRC . 

GET f PREMET. 

FTN j R = 3 r B”0BPRE f I=-=PREMET . 

LABELr TAPE 1 tUSN = ST 1412 fLB = KU,D = 800rF = S. 

FI LE f TAPEl f RT-F r MBL-BO . 

LDSETf FILES = TAPE1 T PRESETA=-NGINFvHAP = SBEX. 
LOADrOBPRE. 

OBPRE. 

REWINDf OUTl *MET1 T MET2rMET3. 

COPY f OUT 1 . 

COPY » MET 1 . 

C0PY^MET2. 

C0PYrMET3. 

REPLACE r OUT 1 y MET 1 »MET2fMET3. 

DAYFILErBAYFX. 

REPLACEt BAYFX. 

EXIT. 

REW I ND y OU i 1 y MET 1 y MET2 y MET3 . 

COPY y OUT 1 . 

COPY y MET 1 . 

C0PYyMET2. 

C0PYyMET3. 

R E P L. A C E y 0 U r 1 y M E T 1 y M Ef. T 2 y M E T 3 . 

DAYFILEyDAYFX. 

REPLACE yliAYFX. 

/ 


PRICE 
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GE0I4AGNETIC INDICES FILE (KP) 


A geomagnetic planetary indices file (KP) is required for 
the Jacchia-Roberts model in LAIRS. This file contains eight 
Kp values per day taken at 3-hour intervals . It also contains 
the daily solar flux number. Three LAIRS routines were 
developed to generate this file for a period of 1 year. 


Program JRKP generates the file (KP) . The required 

K and solar flux data is obtained from the technical journal 
P 

entitled, "Solar Geophysical Data Prompt Reports," which is 
issued monthly. The input data for Program JRKP must be in 
the following format: 

Card 1 

Columns Format Description 

1-6 16 Total number of input data cards. 


Card 2 through n 


Columns 


Format 


Description 


1-6 16 Date, YYMMDD . , corresponding to 

the solar flux and K values. 

P 

8-12 F5.1 Flux number 

14-38 8(I1,A2) Eight K^ values with their signs. 

Routine JRREAD was developed to read the KP file and dump 
it for inspection. It lists the K values in the processed 
format, not in their raw format. 


A sample deck set-up is shown below to execute either of 
these options. When a new K^ file is to be generated, the 
user has to execute Program JRKP with a full year's worth of 
data. This restriction exists because of dimension statements 
in the coding. 
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To create the K file (KP) 

P 

GET , JRKP , KPDATA . 

FTN,R=3,I=JRKP. 

LG0 

REWIND , KP . 

REPLACE, KP. 

To dump the contents of the K^ file 
GET , JRREAD , KP . 

FTN , K= JRREAD , L=0 , TS , B=0B JRAD . 

L0AD , 0B JRAD . 

KPDATA is a file consisting of card images in the described 
format, containing the raw K^ and solar flux data. KP 
is a direct access file containing this data in a form 
suitable for use in the LAIRS Program. 
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APPENDIX C 


LAIRS BENCHMARK DECKS 


Figure C-1 shows a set of LAIRS keyword card decks used as 
benchmark tests during LAIRS development. These decks are 
set up for a LAIRS stacked deck run and illustrate the cards 
that exercise the various LAIRS options. Each deck begins 
with an ENTREE card and ends with an END card. The decks 
are processed in sequence until the FIN card is reached. 

The first deck calls for each twentieth point on the input 
ENTREE file to be processed. The Jacchia-Roberts model is 
used in the upper atmosphere, and default interpolation 
(i.e., interpolation from the working file) is used in the 
lower atmosphere. Default interpolation is used for winds 
in both regions. Latitude variations and diurnal/semidiurnal 
effects are selected for the lower atmosphere (they are 
automatically included in the Jacchia-Roberts model) . The 
interpolator type is second-order Lagrangian univariate in 
altitude. 

The second deck is similar to the first, except that inter- 
polation from meteorological files has been selected as the 
model to be used in the lower atmosphere (and for wind in 
both atmospheric segments) . The two STATION cards indicate 
that two meteorological files are to be used. The first 
STATION card, which corresponds to the file input on tape 21, 
is labelled with a GMT of 8.0 hours. The second STATION card 
corresponding to the file on tape 22, has a GMT of 11.0 hours 
The altitude interpolator type is second-order Lagrangian 
univariate, and the horizontal interpolator type is first- 
order Lagrangian bivariate. The horizontal interpolation 
eliminates the need to specify latitude and diurnal/semi- 
diurnal corrections on the LOWMOD card. 

The third deck calls for the polynomial model to be used in 
the lower atmosphere. In the simple linear mode selected 
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in this deck, temperature, pressure, and density measurements 
will be independently fitted by separate polynomials. 

The fourth deck is a duplicate of the third, except that 
pressure and density calculations will be performed by an 
integration of temperature at every ENTREE point processed. 

The reference value of pressure needed for this integration 
will be interpolated from the meteorological files at 30 kilo- 
meters altitude. The integration step size will be 0.5 
kilometers . 

The final deck calls for a nonlinear differential correction 
fit. All three data types (temperature, pressure, and density) 
will be used in this fit, although only the temperature poly- 
nomial coefficients will be produced. Calculation of pressure 
and density at each ENTREE point will proceed via integration, 
as was the case for the fourth deck. However, no DERIVE card 
needs to be included, as the values needed for integration 
will be set by the polynomial fitting routines. 

Figure C-2 shows an example of a deck using the Local Atmos- 
pheric Profile (LAP) option, which will produce a vertical 
test profile of temperature, pressure, density, and winds. 

The date of the profile is April 12, 1980, at 10 hours, 

0 minutes, 0 seconds GMT, The profile will be produced at 
latitude 35 degrees north and longitude 220 degrees east 
of Greenwich, and will extend in altitude from 0.0 to 120.0 
kilometers. Values will be printed every 2.0 kilometers. 

Figure C-3 shows two card decks used in the optimal procedure 
described in Section 5. The first deck interpolates meteo- 
rological data from the location of the meteorological pro- 
files to the Shuttle trajectory. The keyword card CREATE 
causes the values resulting from this interpolation to be 
output on a new file in the same format as a LAIRS meteoro- 
logical file. 
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The second deck uses this new file as input to the LAIRS 
Program and performs a differential correction on the data. 
Since all diurnal, semidiurnal, and latitude effects have 
already been included in the first run, they are turned off 
in this second run. A STATION keyword card is needed so 
that the system will recognize that the source of the data 
is a meteorological file, but the GMT field on the card will 
be ignored. 

A sample job stream is shown in Figure C-4. The files used 
are ; 

• CARDS - the keyword card images (the deck shown in 

Figure C-1, for example) 

• PERMl - The permanent file 

• ENT - The ENTREE file (the Best Estimate Trajectory 
(BET) file) 

- The meteorological files produced by the 
preprocessor PREMET 

• KP - The Jacchia-Roberts geomagnetic index and 

solar flux file 

All files beginning with "OBJ" are object modules. 

The job streams required to run the decks in Figure C-3 are 
essentially the same as the job stream of Figure C-4. For 
the first deck of Figure C-3, however, an additional file 
(USEMET) is created during the run and must be saved after 
execution. This file is simply renamed and substituted for 
the meteorological files METl, MET2, and MET3 in the run 
controlled by the second deck. 


• METl 

• MET 2 

• MET3 
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ENTREE 

20. 





UPMOD 

3.0 

1.0 

0.0 

0.0 

0.0 

LOWHOD 

1.0 

1.0 

1.0 

1.0 

0,0 

INTERPOL 

2.0 

0.0 




E^;D 






entree 

20. 





IJPNOD 

3.0 

2.0 

0.0 

0.0 

0,0 

LOUNOD 

2.0 

2.0 

0.0 

0.0 

0,0 

INTERPOL 

2.0 

0.0 




STATION 

8.0 





STATION 

11.0 





END 






ENTREE 

20 . 





UPMQD 

3.0 

2.0 

0.0 

0.0 

0.0 

LONMOD 

2.0 

2.0 

0.0 

0.0 

0.0 

INTERPOL 

2.0 

2.0 




STATION 

8.0 





STATION 

11.0 





END 






entree 

20. 





UPMOD 

3.0 

2.0 

0.0 

0.0 

0,0 

LOWMOD 

4.0 

2.0 

1.0 

1.0 

0.0 

STATION 

8.0 





STATION 

11.0 





END 






entree 

20. 





UPMOD 

3.0 

2.0 

0.0 

0.0 

0.0 

LDU'MOD 

4.0 

2.0 

1.0 

1.0 

0.0 

STATION 

8.0 





station 

11.0 





DERIVE 

0.5 

2.0 

30.0 



END 






ENTREE 

20. 





UPMOD 

3.0 

2.0 

0.0 

0.0 

0.0 

LOUIMOD 

4.0 

2.0 

1.0 

1.0 

0.0 


DC 

STATION 8.0 

STATION 11.0 

END 
FIN 


Figure C-1. LAIRS Benchmark Keyword Card Decks 
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LAP 

800412. 

100000.0 

35.0 

220.0 

0.0 

12 ). ■> 

2.0 

UPMDD 


3.0 

1.0 

0.0 

0.0 

0.0 


LOUMOD 


1.0 

1.0 

1.0 

1.0 

0.0 


INTERPOL 


2.0 

0.0 





END 








LAP 

800412. 

100000.0 

35.0 

220.0 

0.0 

120.0 

O 

CM 

UPMOD 


3.0 

2.0 

0.0 

0.0 

0.0 


LOWMOD 


2.0 

2.0 

0.0 

0.0 

0.0 


INTERPOL 


2.0 

0.0 





STATION 


8.0 






STATION 


11.0 






END 








LAP 

800412. 

100000.0 

35.0 

220.0 

0.0 

120.0 

2.0 

UPMOD 


3.0 

2.0 

0.0 

0.0 

0.0 


LOWMOD 


2.0 

2.0 

0.0 

0.0 

0.0 


INTERPOL 


2.0 

2.0 





STATION 


8.0 






STATION 


11.0 






END 








LAP 

800412. 

100000.0 

35.0 

220.0 

0.0 

120.0 

2.0 

UPMOD 


3.0 

2.0 

0.0 

0.0 

0.0 


LOWMOD 


4.0 

2.0 

1.0 

1.0 

0,0 


STATION 


8.0 






STATION 


11.0 






END 








LAP 

800412. 

100000.0 

35.0 

220.0 

0.0 

12C .( 

2.0 

UPMOD 


3.0 

2.0 

0.0 

0.0 

O.C 


LOUMOD 


4.0 

2.0 

1.0 

1.0 

o.c 


STATION 


8.0 






STATION 


11.0 






DERIVE 


0.5 

2.0 

30.0 




END 








LAP 

800412. 

100000.0 

35.0 

220.0 

0.0 

120,0 

2.0 

U'^-rtOD 


3.0 

2.0 

0.0 

0.0 

0.0 


LOUMOD 


4.0 

2.0 

1.0 

1.0 

0.0 



nc 

STATION 0.0 

STATION 11.0 

END 
FIN 


Figure C-2 . 


Sample LAIRS Deck Using LAP Option 
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ENTREE 

810414 * 

12t0 

60,0 

95,00 



UPMOD 


2t0 

2,0 

1,0 

1,0 

0,0 

L.0WM0D 


2,0 

2,0 

1,0 

1,0 

0,0 

INTERPOL 


2,0 

1,0 




STATION 


17,083 





STATION 


17,2167 





STATION 

CREATE 

END 

FIN 


18,367 





ENTREE 

810414, 

00,00 

00,0 

000,0 



UPMOD 


3,0 

2,0 

0,0 

0,0 

0,0 

LOWMOD 


4,0 

2,0 

0,0 

0,0 

0,0 

MODREF 


90, 

68, 

14, 

4, 


INTERPOL 


2,0 

1,0 




DC 

STATION 

END 

FIN 


0,0 






Figure C-3, Sample LAIRS Decks for Two-Step Modeling 



/JOB 

LAIRSf T200 tCH1.40000 » 

USERjXXXXXXX* 

CHARGErXXXXXX? LRC* 

GET y CARDS y PERM 1 f ENT. 

GETy TAPE2:l.===:METl. . 

6ETrTAPE22===^ME'T2. 

GETyTAPE23==MET3. 

6ET»0BJJR>KP. 

GETyOBJUTLyOBJWND. 

G E T y (3 B J C 0 M y □ B J 1.. R y 0 B J I □ y 0 B J F- 1 L y 0 B J 1 N 'T . 

GET 7 OB JPOL y QB JABJ y QEJBC v 
MAP y OFF- . 

CTIME. 

ST I ME. 

ATTACF-ly FTNML I B/UN^L I BRAR Y . 

LDSETyL:tB==FTNMLIB. 

L 0 A B y G B J C G M y (] B .J I... R y 0 B J I G y G B .J J R y G B J F- 1 1... y 0 B J I N T y 0 B J U T L y G B J W N D . 
LGAD y OB JPOL y OB JAD J y OB JBC . 

EXECUTE. 

DAYFILEy DAYFLE . 

REPLACE yDAYFLE. 

REWIND y OUT . 

COPY y GUT. 

REPLACE y OUT . 

REWIND y USE . 

REPLACE y USE. 

EXIT . 

DAYFILEy DAYFLE. 

REPLACE yDAYFLE. 

REWIND y OUT. 

COPY y OUT. 

REPLACEyOUT. 

/ 


Figure C-4. Sample LAIRS Job Stream 



^^ENDIX D_- FILE M AINTENANCE ROUTINE DESCRIPTIONS 

This appendix contains descriptions of the preprocessor rou- 
tines used for LAIRS file maintenance. These descriptions 
are in alphabetical order. 
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FUNCTION FLUXAV (IlDAYf I81DAY) 

purpose: 

FLUXAU CALCULATES THE AVERAGE VALUE OF THE OBSERVED OTTAWA 2800 FLUX 
FOR AN 81 DAY PERIOD BEGINNING WITH IlDAY AND ENDING WITH I81DAY, 

method: 

1, ADD THE FLUX VALUES FOR 81 DAYS AND STORE IN VARIABLE 
FLUXSM 

2. DIVIDE FJ-UXSM BY 81 TO GET THE 81 DAY AVERAGE 
ARGUMENT LIST: 

VARIABLE TYPE DESCRIPTION 

IlDAY I 40 DAYS BEFORE THE DAY IN QUESTION 

I81DAY I 40 DAYS AFTER THE DAY IN QUESTION 

CALLING subroutine: JRKP 

SUBROUTINES CALLED! NONE 

COMMON BLOCK PARAMETERS: 

VARIABLE TYPE DESCRIPTION 

FLUX R OTTAWA 2800 FLUX NUMBER 

programmer/analyst: a* K. KAPOOR, computer SCIENCES CORPORATION 


Function FLUXAV 
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PROGRAM JRKP ( INPUT » OUTPUT r URIiATA » KPI LE r TAPE5 = JRDATA » T APE75 = KP I LE ) 

purpose: 

THIS PROGRAM BUILDS A JACCHI A-ROBERTS ATMOSPHERE FILE (FILE 75) 

IN GTDS DATA SET LAYOUT FORMAT, CURRENT DIMENSIONS ARE FOR A 
FILE LENGTH OF ONE YEAR* 

method: 

1* READ YEARr MONTH» DAY (IN PACKED FORM* I.E** YYMMDD)* DAILY 
OBSERVED OTTAWA SOLAR FLUX AND THE THREE-HOURLY KP GEOMAGNETIC 
ACTIVITY INlrEX* I IS THE NUMBER OF DATA CARDS READ* 

2* CALCULATE THE MODIFIED JULIAN DATE (MJD) OF PACKED DATE 
JUST READ* 

3* CHECK IF DATA CARDS ARE IN ORDER* I*E** THAT DATES ARE 
INCREASING BY ONE AS EACH CARD IS READ* 

4* SET UP THE INFORMATION FOR THE HEADER RECORD* IDRECl IS THE 
MJD OF THE FIRST DAY OF AVAILABLE DATA* N THE FILE* IDRECl 
MUST BE A MULTIPLE OF 20* 

5* LOOP TO FINISH READING ALL THE INPUT CARDS* 

6* DETERMINE THE DATA RECORDS IN THE FILE NREC . 

7* DETERMINE IF THE LAST RECORD IS A FULL 20 DAYS* 

8* DETERMINE THE MJD OF THE LAST DAY OF AVAILABLE DATA IN THE 
FILE (IDRECE) AND THE MJD OF THE LAST GOOD DATA POINT(LGD)* 

9* WRITE THE HEADER RECORD. 

10* CALCULATE TC ARRAY* WHICH IS COMPOSED OF TWENTY DAYS OF NIGHT- 
TIME MINIMUM GLOBAL EXOSPHERIC TEMPERATURES* STARTING WITH DAY 
IT* CALCULATE KP ARRAY* WHICH IS COMPOSED OF TWENTY-ONE DAYS OF 
PACKED KP3HR VALUES STARTING WITH DAY IT-1* 

11* ZERO THE KP AND TC ARRAYS SO THAT IF THE LAST DATA RECORD IS 
LESS THAN TWENTY DAYS* THE U^|IFILLED PORTION OF THE RECORD 
WILL BE PADDED WITH ZEROES* 

12* WRITE A DATA RECORD* NRECWR COUNTS THE NUMBER OF RECORDS WRITTEN* 
13* UPDATE RECORD NUMBER AND IF EQUAL TO NREC* STOP* IF NOT* THEN 
GO BACK TO POINT (10) AND CONTINUE PROCESSING. 

ARGUMENT LIST: 

CALLING subroutine: 


SUBROUTINES CALLED: 


FLUXAV 


COMMON BLOCK PARAMETERS: 

VARIABLE TYPE DESCRIPTION 

FLUX I OTTAWA 2800 FLUX NUMBER 

EXTERNAL DATA SETS USED: 

NAME LUN I/O OPERATION PERFORMED 

KP 75 WRITE THE FILE TO FILE 75 


PROGRAMMER/ANALYST 


A* K* KAPOOR* COMPUTER SCIENCES CORPORATION 


Program JRKP 
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PROGRAM JRREAD ( INPUT OUTPUT » KPILE r TAPE75 = KPILE ) 


purpose: 

TO LIST THE KP FILE DATA SET BY EACH RECORD 


method: 

1. read the header RECORD. NRECR COUNTS THE NUMBER OF RECORDS 
READ 

2. PRINT THE CONTENTS OF HEADER RECORD 

3. READ THE DATA RECORDS 
4* PRINT THE DATA RECORDS 

5* ADUANCE THE NRECR COUNTER AND GO BACK TO STEP 3. IF LAST 
RECORD» THEN STOP* 


ARGUMENT LIST: NONE 


CALLING subroutines: 


SUBROUTINES CALLED: OPENMS » RE ADMS t CLOSMS 

COMMON BLOCK PARAMETERS: NONE 

EXTERNAL DATA SETS USED: 

NAME LUN I/O OPERATIONS PERFORMED 

KP 75 READ ONLY 

PROGRAMMER/ANALYST : A. K* KAPOOR» COMPUTER SCIENCES CORPORATION 


Program JRREAD 
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PROGRAM PREMET ( I NPUT» OUTPUT » OUT 1 fMETl,MET2FMET3f 

1 TAPEl=/80f TAPE8=0UT1»TAPE25=MET1»TAPE26=MET2» 

2 TAPE27=MET3) 

PURPOSE} 

TO READ ANY STRANGE METEOROLOGICAL DATA TAPE AND 
OUTPUT PROCESSED FILE(S) FOR LAIRS PROGRAM TO USE, 
THESE FILES ARE IN A FORMAT COMPATIBLE WITH THE 

Cairs read routine, 
method: 

1. READ THE TAPE AND EXTRACT INFORMATION NEEDED BY 
LAIRSJ SUCH AS NUMBER OF FILES (PROFILES ON TAPE 

NUMBER OF POINTS/PROFILE 
STARTING POINT/PROFILE 
END POINT/PROFILE 

2. READ THE DATA AND OUTPUT AN INVERTED PROFILE FOR 
EACH PROFILE (FILE) 


3. CONVERT THE DATA TO METRIC UNITS 

4. WRITE A SINGLE HEADER RECORD WITH THE INFORMATION 
STATED IN POINT 1 

THEN WRITE THE NUMBER OF FILES DEPENDING ON THE 
NUMBER OF PROFILES PROVIDED ON THE ORIGINAL TAPE 

ARGUMENT LIST} NONE 

CALLING routines; NONE 

SUBROUTINES CALLED? NONE 

COMMON BLOCK PARAMETERS USED? NONE 

EXTERNAL DATA SETS USED? 

NAME FILE ♦ DESCRIPTION 

1 /80 ORIGINAL TAPE FROM LARC 

25 METl FIRST FILE WITH HEADER 

26 MET2 SECOND FILE 

27 MET3 THIRD FILE (PROFILE) 

PROGRAMMER/ANALYST J A, K. KAPOOR i COMPUTER SCIENCES CORPORATION 


Program PREMET 


D-5 



DIFFERENTIAL CORRECTION NOISE ANALYSIS 


APPENDIX E - 


A series of tests of the differential correction process 
using simulated data was conducted to determine the response 
of the program to various levels of random noise. A non-DC 
polynomial fitting run was used with default CIRA data to 
produce a simulated meteorological file via the CREATE key- 
word option. This file was used as input to a modified 
version of LAIRS which, prior to fitting, added a specified 
amount of random noise to the simulated data. The LARC 
library routine URAN, which returns a uniformly distributed 
random number between 0 and 1, was used in this process. 

The random number was shifted to lie in the interval between 
-0.5 and 0.5 and scaled by a factor supplied in each test 
to yield noise with average amplitudes ranging up to 40 per- 
cent of the values of the parameters on the file. 

It was found that, for relatively low noise levels, the dif- 
ferential correction process consistently yielded excellent 
results, with the estimated parameters often being within a 
few percent of the known values even when the noise amplitudes 
reached 20 or 25 percent. For larger noise amplitudes, the 
results were generally good although, of course, larger vari- 
ations in the values of the estimated parameters tended to 
accompany these larger noise amplitudes. In all test cases, 
the results were well within the range of acceptability for 
the noise levels and number of data points used (generally 
30 to 50 points per segment) . 

The LAIRS differential correction process proved to be very 
stable in all the tests. High levels of random noise did 
not result in excessive editing of data points or in runs 
that would not converge in a few iterations . 
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